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FOREWORD 


Sir Arthur Eddington once remarked that one of the most 
profound mysteries of the universe consists in the fact that 
everything in it rotates. The meteors, asteroids, the planets and 
satellites, the Sun, the stars, the clusters, and nebulae, and 
even the galaxies spin around their axes. As they do so their 
shapes are distorted—^imperceptibly in the case of the Sun with 
its equatorial rotational velocity of but 2 km./sec., but easily 
distinguishable, with a small telescope, in the case of Jupiter or 
Saturn. 

There are many hot stars whose equatorial rotational velocities 
approach values of the order of 500 km./sec. In a few white 
dwarfs even larger velocities may be found. Such stars must be 
greatly flattened, and in a few of them we actually observe the 
outflow of gases at the sharp edges of their equators, producing 
flat, extended gaseous rings in their equatorial planes. 

Professor Jardetzky's book deals with the figures of equilibrium 
and distortions of rotating bodies, and it comes at a very opportune 
time. This theory is usually associated with the names of Poincar6, 
Darwin, Liapounov, and Jeans, although many other distinguished 
astronomers and mathematicians—among them Dr. Jardetzky 
himself—have made contributions to it. In its early stages the 
aim of the theoreticians was to trace the consequences of an ever 
increasing angular rate of rotation. This problem has assumed 
new importance in stellar astronomy because it is now possible 
to trace, largely from the work of L. G. Henyey in Berkeley and 
of B. Strdmgren at the University of Chicago, how a diffuse 
gaseous mass condenses through the gravitational process known 
as Kelvin contraction, and becomes after a lapse of some hundreds 
of thousands of years, a main-sequence star. The angular mo¬ 
mentum of the original cloud is preserved in the contraction, and 
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the newly-formed star rotates around its axis. But even before 
the star reaches the main sequence nuclear processes start operat¬ 
ing in its interior, and they maintain the star for a long time— 
millions, billions, or even tens of billions of years—in equilibrium 
on the main sequence. Gradually, however, the hydrogen nuclear 
fuel begins to be nearly exhausted. The star then starts expanding 
again—at first slowly, by a process first investigated by S. Chan¬ 
drasekhar and M. Schonberg, and later more and more rapidly 
according to a new theory by M. Schwarzschild and A. Sandage. 
The star blows up and cools down: it becomes a diffuse giant. 
Its surface rotation must change again, and accordingly also the 
shape of its surface. The observations show that many of these 
'"old’" giants rotate much faster than had been expected; values 
of the order of 75-100 km./sec. are not unusual. It should be 
possible to infer from these velocities the distribution of the 
material inside such giant stars. 

Professor Jardetzky was the last of many pupils of the late 
A. M. Liapounov whose work he has treated in this volume. In 
1918, during the Russian revolution, Liapounov settled in Odessa 
where Jardetzky was then a young mathematician. In the last 
year of his life (he committed suicide in 1918) Liapounov gave 
a series of lectures at the university about his investigations, 
which have since been published by the Soviet Academy of 
Sciences, in 1925, under the title Sur certaines sSries de figures 
d'equilibre d’un liquide heterogdne en rotation. In this memoir he 
gave the most general results in the field in which his Master’s 
thesis, in 1884, on the stability of ellipsoidal forms of equilibrium 
of a rotating liquid won for him immediate international recog¬ 
nition. It should be emphasized that many of the later results of 
Liapounov have remained virtually unknown in the western world. 

Astronomers may be interested in a few personal remarks. A.M. 
Liapounov was the son of the astronomer M. Liapounov, who 
was in charge of the Kazan Observatory around the middle of 
the 19th century. There the father was engaged in various astro¬ 
nomical observations, including a detailed study of the positions 
of stars associated with the Orion nebula. This work was published 
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by him in collaboration with Otto Struve Sr. in 1862 . The son 
learned the fundamentals of astronomy from his father. His 
mathematical inclination was, however, the result of the inspi¬ 
ration he received from the great mathematician Tchebyshev at 
St. Petersburg University, The son Liapounov spent some of the 
most productive years of his life as professor of analytical me¬ 
chanics in the University of Kharkov, where he was associated 
with the father of the present writer. Unfortunately, the latter 
never had the opportunity of studying under Liapounov. His 
knowledge in this field derives from one of Liapounov’s earliest 
students, N. N. Saltykov. 

Otto Struve 

Professor of Astrophysics, University of California; 

Director of the Leuschner Observatory, Berkeley, California 
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PART I 


Figures of Equilibrium of a Rotating Fluid Mass 




CHAPTER I 


IntroJttction 

The theory of figures of celestial bodies has been developed 
by a nuruber of most outstanding mathematicians, primarily by 
Maclaurin, Jacobi, Poincar6, and Liapounov. Their investiga¬ 
tions, however, have been largely limited to the case of an isolated 
fluid mass in uniform rotation. This approach was justified for 
two reasons. First, the figures of celestial bodies actually differ 
but little from the figures given by the solutions of this problem 
(e.g., a sphere, ellipsoids of Maclaurin and Jacobi, figures of 
Poincar^ and Liapounov), and second, the existence of such 
figures could be demonstrated by the most rigorous methods of 
mathematical analysis. 

Some progress may still be possible in this area, but the develop¬ 
ment of an exact theory which would take into account conditions 
more closely approaching reality appears more interesting and 
important. Variations in figures and stratification of celestial 
bodies which result from internal movements should be considered 
as the next approximation. There are variations important 
primarily in the investigation of stellar structure. On the other 
hand, there are features observed in the Earth's crust and similar 
features undoubtedly existing in the crust of any other solidified 
planet which may be also explained by the departure from 
conditions of equilibrium. Thus, many problems arise in which 
the figures of celestial bodies have to be determined under con¬ 
ditions different from equilibrium. 

The most exact methods used in the theory of figures of equilib¬ 
rium and the most important results are discussed in the first 
six chapters. The second part of this book deals with some 
other'problems concerning the figures of celestial bodies. Many 
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of these problems still await their solution by equally rigorous 
methods as those developed in the theory of figures of equilibrium. 


LL Fundamental Problem 


The solution of the following problem of Hydrodynamics is 
being presented for the purpose of explaining the shape of the Sun 
and other celestial bodies. Let a fluid mass be rotating in space as 
a rigid body. Let us assume that it is isolated from other bodies 
and that its particles are subject to mutual attractions according 
to Newton’s law of gravitation. What is the figure of this mass? 

When the problem is given in this general form, not all condi¬ 
tions are taken into account. Nevertheless, we can make some 
important preliminary conclusions using the equations of motion 
of a fluid mass. We call ^-axis the axis of rotation of the mass and 
take the origin of a rectangular system of axes Oxyz at the mass 
center. Let to be the angular velocity. The velocity v of a point M 
is then given by the expression 

(1.1.1) V = to X r 

if the vector r[x, y,z) is equal to OAi. The components of ac¬ 
celeration at M are: — — co^y, 0. Let fU be the gravitational 

potential, / the constant of gravitation and and p the density 
and pressure of the fluid. We can write the equations of motion 
of a fluid mass either in the vector form 

(1-1-2) — = grad/[/--grad 

or in the scalar form 


(1.1.3) 


dU 

-co^x = f — 
ox 


1 dp 
X dx’ 


(o^y = / 


du 

dy 


1 dp 
% dy' 


0 = /-^ 


1 dp 
K dz 


If we assume that co is constant, we obtain by (1.1.3) 


(1.1.4) 


fdU ^ — d{x^ + y^) dp 
2 % 
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In two cases, namely, k = const, ox n ~ (p{p) we have an inte¬ 
gral of this equation in the form 

(1.1.5) fU H-— = — + const. 

2 J >c 

where s = Vx^ + is the distance of a particle from the axis 
of rotation. 

For Newton's law of gravitation we have^ 

,o.e, 

where D is the distance of the particle at M from an element dV'. 

If a fluid mass is isolated in space, the pressure at its boundary- 
surface vanishes: 


(1.1.7) ^ = 0 

Thus the boundary surface belongs to the set of surfaces of equal 
pressure. The general condition for this set may be immediately 
obtained by putting p = const, in (1.1.5) 


( 1 . 1 . 8 ) 


' >c'dV' 
V D 


+ 


co^s^ 

—— = const. 

2 / 


In the case x = <p(j>) the surfaces p — const, coincide with the 
surfaces of equal density. 

In general, the volume V of the fluid mass in the equation 
(1.1.8) is unknown and the shape of this mass must be chosen in 
such a way that this equation will be satisfied. Thus, the problem 
we have posed in its simplest form, leads to the solution of a 
functional equation (1.1.8), where the limits of integration in the 
first term are unknown.^ 

We obtain a particular case of the problem on considering figures 
of equilibrium of a fluid mass at rest. Since we assume now that 


^ The theory of figures of celestial bodies can be readily extended to the 
case of a law of gravitation which yields a potential differing but little from 
(1.1.6) as it was shown by Jardetzky [5]. 

^ This equation can be transformed in an integral equation by using the 
method of Dirichlet. The resulting equation is a nonlinear one with four-fold 
integrals having limits i oo and a very complicated integrand. 
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cw = 0, these figures must be determined by the solutions of the 
equation 

f H'dV' 

(1.1.9) —— =const. 

Jv D 

We have to keep in mind that on deriving the integral (1.1.5) 
we have considered fluids as having the property that either 
X = const., i.e., homogeneous and incompressible fluids, or 
^ i.e., compressible fluids which are homogeneous by 

composition but can have a density varying with pressure. It 
is evident, however, that the equations (1.1.3) hold without these 
restrictions, and we can investigate also the figures of equilibrium 
of a fluid mass that is composed of a finite or even of an infinite 
number of incompressible fluids having different densities or 
figures of a compressible and inhomogeneous mass. 

There are types of motion other than a uniform rotation of a 
fluid mass about a fixed axis which are also determined by 
equations (1.1.3). The problems concerning these more general 
movements will be discussed later. 

L2, Remarks Concerning Basic Assumptions 

In order to solve a problem in the theory of figures of equilibrium 
or, in general, to determine the figures of a rotating fluid mass, 
certain basic assumptions are made concerning the physical 
properties of the fluid, the character of motion, the mass distribu¬ 
tion, and others. Some of these assumptions will be discussed now. 

For the motion of a fluid the condition of continuity can be 
written either in the form 

^ + K div V = 0 
dt 


3 ^ 

— + V ■ grad x -j- ?€ div v = 0 
at 

In case of a permanent motion, we have dxjdt = 0, i.e., the 
density x does not depend on time explicitly. But, in some 


( 1 . 2 . 1 ) 

or 

( 1 . 2 . 1 ') 



INTRODUCTION 


7 


invariable figures of a fluid mass, the density of a particle will 
not vary along its path and, therefore, we have dKjdt = 0. By 

(1.2.1) and (1.2.1') we obtain div v = 0 and 

(1.2.2) V • grad k = 0 

Since the vector grad x is perpendicular to the surface of equal 
density = const., the velocity v must be in the tangential 
plane passing through the given particle or point in the fluid. 

If the motion of a fluid mass has the characteristics just 
mentioned, the equation of continuity has the form 

(1.2.3) divv^O 


for an incompressible fluid as well as for a compressible one. It 
is easy to see that, if the velocity v is given, by the expression 
(1.1.1), the equation (1.2.3) will be satisfied either by co = const, 
or by CO = F{x^ + z). It is, namely, 


div V = 


d{o)y) d(cox) 

dx dy 


in both cases. Thus, the assumption of a rotating figure of 
equilibrium is compatible with the condition (1.2.3). 

The second assumption about the character of motion concerns 
the existence of an axis of rotation. It has been proved by 
Appell [1] that a necessary condition of the existence of a figure 
of equilibrium of a fluid mass is that the rotation occurs with 
respect to an axis having a direction invariable in space and 
coinciding with one of the principal axes of inertia of this mass. 
Jardetzky [3] has shown that this condition holds also for a 
more general kind of rotation, namely, for the so-called zonal 
rotation and for mixed systems, such as a body composed of solid 
and fluid parts. 

It is not necessary to make special assumptions about the 
symmetry of figures of equilibrium. It is known (see, for example, 
Appell [1]) that the axis of rotation can be but is not necessarily 
an axis of symmetry of a rotating fluid mass. It can be shown, 
for example, in case of a homogeneous mass that there can be 
a finite number of planes of symmetry passing through the axis. 
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On the other hand, Lichtenstein [1] has proved that a figure of 
equilibrium must have a plane of symmetry perpendicular to 
the axis of rotation. 

As to the assumptions concerning the physical properties of 
the fluid, few remarks have to be added. We have seen that we 
have to discuss the figures of a homogeneous or heterogeneous, 
compressible or incompressible, fluid mass. It is usually assumed 
that the fluid is isotropic. A theory of figures of an anisotropic 
fluid has not yet been developed. 

In the most simple case, a fluid mass is considered incompressible 
and homogeneous in the sense that the density has equal values 
everywhere and that these values cannot be changed by a varying 
pressure. Then, the equation (1.1.8) takes the form 


(1.2.4) 


I, 


dV' 

v~D 


+ 


co^s^ 


= const. 


The left side of this equation is an unknown function of coor¬ 
dinates determined by the shape of the boundary surface S of 
the volume V. If the value of the constant right-side member 
is taken for p = 0, the problem of figures of equilibrium is to 
find a set of surfaces S such that for the corresponding volume V 
the condition (1.2.4) will be satisfied, i.e., that (1.2.4) will be the 
equation of any of these surfaces. 

The second problem in the theory of figures of equilibrium deals 
with the case where the density depends on the pressure only 
x = <i>{p). In classical Hydrod 3 niamics, these fluids are usually 
referred to as the inhomogeneous. To solve the problem in this 
case, we have to make use of equation (1.1.8) or (1.1.5), since 
the last equation yields the level surfaces. Under our assump¬ 
tion about the physical properties of the fluid, the surfaces 
p —■ const., X = const., and the level surfaces (the surfaces of 
equal gravity potential in case of the Earth) will coincide if 
(1.1.6) is satisfied. In this equation, the density jc is the variable 
which will determine the stratification. The outer layer of the 
fluid corresponds to x = q>{0). 

As to the equation of state x = <pip), it can determine a con- 
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tinuous function, or we can consider a system of a finite number 
of fluids differing in their physical properties superimposed on 
each other. For each layer, then, a particular equation of state 
will hold. A fluid for which no assumption such as = const. 
or H = (p{p) is made can be also considered. Then a continuous 
stratification is postulated with a density expressed in terms of 
coordinates. For such fluids the conditions (1.1.4) and (1.1.5) 
do not always hold and, therefore, the discussion must begin 
with the equations (1.1.3), 

In problems concerning the figures of equilibrium mentioned 
before, the stratification of an inhomogeneous fluid is not known 
in advance and is to be found. The figure of such a fluid is 
determined by the shape of the outer layer. Of course it is not 
always necessary to follow this direct way. One can make 
assumptions about the shape of the boundary surface and postulate 
certain stratification. Then, it must be verified that all other 
conditions are not violated by the assumptions made. 

Volterra [2] has proved that in case of equilibrium the layers 
of equal density cannot have the form of a surface of the second 
degree, and this proof was generalized for a zonal rotation (see, 
for example, Jardetzky, [4]). However, when an approximate 
solution is considered, for example, for a slow rotation for which 
we know that the surfaces of equal density will differ but little 
from a spherical shape, an ellipsoidal stratification presents a 
certain degree of approximation. 

It is usually assumed that in a figure of equilibrium the density 
is increasing towards its centrum. For example, the mass distribu¬ 
tion in the Earth has been represented by the well-known laws 
of Levy, Lipschitz, Roche, and some others of a similar type. 

We shall mention one more assumption, namely, that in the 
case of an inhomogeneous fluid the surfaces of equal density form 
a set of closed surfaces each enclosing the preceding one. 

1.3. Density Distribution 

In the investigations concerning the figures of equilibrium, it 
is usually assumed that the density is an increasing function of 
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the distance from the surface of the body. However, for one 
case, at least, a proof of this fact can be given. Having demon¬ 
strated that a sphere is a unique figure of equilibrium of an isolated 
fluid mass at rest, Liapounov [10] had also shown that the density 
will increase from the surface toward the center because of the 
conditions of equilibrium. His proof is as follows. 

Let A be the radius of a sphere which is the boundary surface 
of a fluid mass at rest and fc the density. It is obvious that if 
an inhomogeneous fluid mass is in equilibrium, it has to be 
composed of concentric spherical layers. In general, the density 
can be a function of a parameter a. However, if the fluid is 
compressible and a = (p{p) it may be shown that « = y (a) cannot 
be an arbitrary function. By equation (1.1.9) we can put on 
a level surface 


(1.3.1) 


' K'dV' 
V D 


const. = a 


It follows then by (1.2.4) and (1.1.5) that 

1 r # 

1 . 3.2 ^ = 

fU(p{p) 

if CO = 0 and ao corresponds to /> = 0, i.e., to the free surface 
of the fluid. Thus, we obtain from this equation p ^ p{oi) and, 
on inserting this function in the equation of state, it is ^ = -^(oc). 
In order to prove that the density is a function increasing toward 
the center, we shall transform the integral (1.3.1). Let be 
the potential at a point M[x, y, z) at a distance r from the center 
0 of the fluid mass and the element >c' dV' taken at a point 
M'[x', y', z') at a distance r\ Taking the polar axis along OM, 
we assume that the angle 6 corresponds to the colatitude and X 
to the longitude. Let 0' be the angle between r and r', do' the 
element of the sphere having unit radius. Then da' = sin 6' dO' dX 
and Z) = 2rr' cos 0'. Now we can write 




(1.3.3) 
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since the volume element dV' at the distance r' is r'^^da'dr'. 
The integral taken over the unit sphere is 


Cda' r 

\-~=27t ■ 

J D Jo 


sin 6' dd' 


271 


_j_ y '2 — 2rr' cos 0' yy 
271 


^ — [y'/S /2 — 2r/ cos 0']J 




+ Ir —/|] 


If y' < f this expression is equal to 

^^da' 4:71 


rda' __ - 


D y 

and if < y' we have 

r da' 4:7t 

For each point M in fluid's interior it is, then, by (1.3.3) 
>€'dV' 47t 


r x'dv' 4:77 , ,, , , , 

—-— = — H r^dr + 4^7 

jv D y Jo J 


x'y'dy' 


IV D y 

Now, if we take all points of a sphere r = a, it is 

xada 


x'dV' 

■ 1 -J— 4-'7 t 

]v D 

1 ✓tfCv {JUiAf j~ jCJlf 

a «/ 0 •' 


where x must be expressed in terms of a. 
Hence 

, d ^x'dV' 

(1.3.4) 


dCx'dV' 4.71 r 

- - =:- \ da 

daj D 


The integral taken with respect to the variable a is obviously 
a positive number because of the condition x> 0, Since 

d C df d C dp dp ^ dp 
daJ (p{p) dpJ cp{p) da x da 
we obtain by (1.1.5) and (1.3.4) that 


(1.3.6) 


dp 


da da 


d Cx'dV' 


J 


D 


< 0 


Thus, the pressure in a spherical figure of equilibrium is 
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increasing toward the center. If we consider the case where the 
density of a fluid is an increasing function of the pressure, the 
density in a figure of equilibrium must also increase from the 

free surface toward the center. 

In case of a heterogeneous and incompressible liquid we usually 
assume that the density is an increasing function of the depth 
and that otherwise the equilibrium would be unstable. 

Another proof of the theorem that a sphere is a unique solution 
of the problem of equilibrium of a fluid mass at rest has been 
given by Carleman [1]. 

1,4. Dev&loptnctit of Methods of Solution 

We have seen that the problem of figures of equihbrium of a 
fluid requires even in its simplest form, the solution of a 

functional equation (1.1.8). This problem is not yet solved in 
its general form, i.e., all the possible figures of equihbrium are not 
yet known despite the fact that the most competent investigators 
have suggested highly ingenious methods. Using these methods, 
solutions were obtained for a set of particular problems. The 
most complete review of the theories concerning the problem of 
equihbrium of a fluid mass is given in Volume IV of the Traitd 
de Mecanique Rationelle by P. Appell. Only a short survey 
necessary for understanding the development of these theories 
will be presented here. 

As is weU known, Newton has seen the importance of the law 
of gravitation for the explanation of figures of celestial bodies 
and, on investigating the conditions of equihbrium in canals 
directed from the Earth’s center toward the equator and the 
poles, reached the conclusion that the figure of the Earth is an 
eUipsoid. Maclaurin was the first who gave (in 1742) a rigorous 
mathematical proof of the fact that an elhpsoid of rotation can 
be a figure of equihbrium of an isolated, rotating, homogeneous 
fluid mass. One year later the investigation of Clairaut dealing 
with an inhomogenous fluid was published. There were no exact 
solutions known in the 18th Century other than Maclaurin’s 
elhpsoid, but attempts to obtain such solutions were made by 
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Maupertius, Simpson, Legendre, and Laplace. It was Legendre's 
suggestion to identify the equation of the free surface of fluid 
with the condition that the sum of potentials due to gravitation 
and to the centrifugal force is a constant on this surface. 

The easiest case in the theory of figures of equilibrium was 
obviously that of a homogeneous fluid (% = const.), and we shall 
mention first the progress made in this case. The potential of 
a homogeneous ellipsoid E corresponding to Newton's law of 
gravitation at a point M{x, y, z) inside the mass can be written 
in the form 

^2 qj2i ^ 2 

(1.4.1) fU = const. - 


Inserting this expression in (1.3.5) of (1.1.8) and taking the 
value of the constant term which the latter has at the surface 
of equal pressure passing through the point M, we obtain 


(1.4.2) (L, - co^)x^ + {L, ^ co^)y^ + = C 


This condition holds for all points of such a surface, and a par¬ 
ticular value of the constant will determine the surface of the 
ellipsoid E. We shall now assume that this surface is given by 
the equation 


(1.4.3) 


^2 4. ^2 ^2 



1 


On putting = Ly = L in (1.4.2) we have for an ellipsoid 
of revolution the equation 
(1.4.2') (L 0)2) (x2 + ^ ^ c 

and, since the equations (1.4.3) and (1.4.2') must represent the 
same surface, we obtain two conditions 
(1.4.4) a^L - 0)2) = = C 

The expressions and Ly are known in terms of the axes of 
the ellipsoid as well as of its density.^ Therefore, we can find the 
corresponding values of o) and C from the conditions (1.4.4). 
Thus, the surface of the fluid mass will be completely determined. 


See (2.2.7). 
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It may be easily proved (as shown on p. 30) that, if an ellipsoid 
will be a figure of equilibrium, it must be an oblate one. It may 
be seen from (1.4.4) and. the expressions for and Ly that in 
different conditions the ratios of axes play an important part. 
If we put I = Va^ — c^jc a condition may be found which w, 
I, and have to satisfy. This condition is 

(3 + l^) tan-i I — 31 
=^- 

The discussion of this equation shows that there is a certain 
limi t value of co and that, in general, three cases must be con¬ 
sidered corresponding to the condition. 


In the first case there are two ellipsoids of Maclaurin which 
correspond to the same value of co. If h has exactly the numerical 
value given in (1.4.6), or is larger, only one such ellipsoid is 
possible. Finally, when o) > coq where coq is the limit value 
mentioned above, no ellipsoid of revolution can be a figure of 
equilibrium of a rotating fluid mass. 

In 1834, Jacobi had shown that an ellipsoid with three different 
axes can also be a figure of equilibrium, if the axes satisfy certain 
conditions. Now, if 


(1.4.7) 


— + —+ - 


is the equation of such an ellipsoid, it must be identical to (L4.2). 
Therefore, the conditions 

(1.4.8) a^{L^ - «2) = b^{Ly - co2) = 

are required. Further analysis of this case was given also by 
Otto Meyer, LiouviUe, Smith, and Plana. We will only mention 
here that L^,Ly, and can be expressed in terms of a, b, c and, 
therefore, if w is given the ratios s = c^ja^ and t — c^jh^ are 
determined by (1.4.8). On making use, then, of the expression 
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for the given mass of a fluid m = ^oxaboc, we can compute the 
axes themselves. A new limit value for the angular velocity is 
determined for these so-called Jacobi’s ellipsoids. This value 
CO == coi, is given by the condition h = 0.1871 .... When 

(1.4.9) h < 0.1871. . . 

there is only one ellipsoid of Jacobi possible which corresponds 
to a given value of co and represents a figure of equilibrium. 
For the limit value of the angular velocity, Jacobi’s ellipsoid 
becomes an ellipsoid of revolution. 

For the discussion of results of the theory the angular momentum 

(1.4.10) ikf = 7co 

(where I is the moment of inertia with respect to the axis of 
rotation) can be used instead of the angular velocity. The 
initial state of a given mass can also be characterized by it. Then, 
the results mentioned above are as follows, (a) If il^ = 0, the 
angular velocity co is equal to zero too, and the corresponding 
figure of equilibrium is a sphere, (b) If the initial angular mo¬ 
mentum M has values between 0 and (or co varies from 0 
to coq — the limit value), the corresponding figure of equilibrium 
is an ellipsoid of revolution, (c) When M belongs to the interval 
(Mq, oo), the corresponding value of co decreases from coq, but 
the corresponding ellipsoids of revolution become more com¬ 
pressed. (d) For M = M-i < or co = co^ the ellipsoid of 
revolution coincides with the limiting ellipsoid of Jacobi. Now 
an ellipsoid having three different axes can be taken into account. 
These ellipsoids can be figures of equilibrium, when M varies 
between and oo, while the angular velocity decreases to 0, 
and the ratios of axes vary between certain limits. 

Whether some figures of equilibrium can exist or not if co surpasses 
the value coq, was an important problem posed by Tchebysheff. 
Liapounov was able to prove in 1884, that there are no new 
figures of equilibrium in the neighborhood of the limit ellipsoid 
(co = coo) but that there exist figures of equilibrium differing 
but little from some well-determined ellipsoids of Maclaurin and 
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Jacobi. As long as M is smaller than M-^, the ellipsoids of revo¬ 
lution are stable. When M is larger than Maclaurin's 
ellipsoids become unstable, but the ellipsoids of Jacobi remain 
stable. On investigating this property, Liapounov found that 
it holds only until a certain limit is reached. From the 
first approximation, the conclusion could be drawn that, if the 
value of M were increased, Jacobi's ellipsoids would change into 
figures of equilibrium of a new kind. These would be algebraic 
surfaces of the third order. One year later than Liapounov, 
Poincar6 also discovered this new kind of figures of equilibrium.^ 
Liapounov expressed his doubt about the possibility of proving the 
existence of new figures of equilibrium by using the method of 
successive approximations, and only after twenty years did he make 
clear the reason for the difficulties met by the theory. In order to 
find a new figure of equilibrium which differs but little from an 
ellipsoid, one has to compare it, not to the given ellipsoid, but to 
a variable ellipsoid which is confocal to the former and passes 
through the point for which the value of potential is considered. 
This was the essential point of the new method developed by 
Liapounov- On using this method, the approximation of an 
arbitrary high order can be evaluated. Then, on applying the 
method of Cauchy, the convergence of successive approximations 
can be shown. Therefore, the existence of these figures of 
equilibrium can be proved with all the required precision. We 
shall later discuss in more detail the figures of a homogeneous 
fluid mass differing but little from the ellipsoides of Maclaurin 
and Jacobi. 

Not connected with this group if investigations was the proof 
given by Sophie Kovalevsky concerning the existence of a 
ring-shaped figure of equilibrium. Despite the fact that a torus 
may be a figure of equilibrium of a homogeneous liquid mass, 
there are no applications of this result, since the figure is unstable. 
To explain the ring of Saturn, one will rather admit the Cassini 
hypothesis according to which this ring is composed of a swarm 


* C. R. Acad. Sci., Paris 1885, and Memoir in Acta Mathematica. 
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of solid particles. (See, for example, Lichtenstein [1], [2].) 

The problem of figures of equilibrium becomes much more 
difficult if an inhomogeneous mass is considered. The first exact 
solution in a case where the density is a function of the pressure 
only was also given by Liapounov. It was published in the 
posthumous memoir of this investigator (1925-1927) and represents 
the final form of his method. Although it is a masterpiece of 
mathematical thought, it is, however, very difficult to study, 
since many details of calculations were omitted. As to the 
classical investigations, the first of them is represented by the 
memoir of Clairaut. The question it deals with is known as 
ClairauLs problem. To determine the compression of level 
surfaces and the distribution of gravity acceleration at the 
surface of a slow-rotating inhomogeneous fluid mass is its main 
purpose, a concentrical stratification being assumed. In the 
same direction the problem was developed by Hamy, V^ronnet, 
Wavre, and Dive, and certain particular cases were considered 
by Wiechert, Klussman, and Haalck. 

The generalization of another kind was given by Liapounov 
[1], [2]. He assumes that a planet, say the Earth, is formed 
by a set of thin layers, each having a constant density and dif¬ 
fering but little from a sphere. Since the density of these con¬ 
centric layers might not be a continuous function of the radius 
of the layer, Liapounov considers an equation more general than 
that of Clairaut. To take into account the discontinuities of 
density, he made use for the first time, in this theory, of integrals 
representing a generalization of Riemann integrals in the sense 
that was determined also by Stiltjes. Liapounov [10] gave the 
following characteristics of methods used in the Clairaut problem. 
It has been known that in the second approximation, the level 
surfaces are determined as certain surfaces of revolution of the 
fourth order. Thus, it is not correct to look again for the elements 
of ellipsoids in the first approximation, since they are elements 
of those ellipsoids which themselves represent in the first approxi¬ 
mation the unknown surfaces. This was the reason why Clairaut 
could not go further than the first approximation. The theories 



18 


FIGURES OF CELESTIAL BODIES 


of Legendre and Laplace were not subject to this objection, but 
Laplace made use of assumptions which can raise certain doubts. 
He assumes, namely, a priori that the unknown function can be 
developed in a series of spherical functions and makes use of suck 
a series the convergence of which is questionnable. If we put, 
in the equation (1.3.3) 

(1.4.11) r=a(l + t), r' = ^'(l + r) 


the potential takes the form 


(1.4.12) 


r K^dV' r 

Jv D Jo J 




D 


Now, Laplace takes 


(1.4.13) 



P„ (cos y') 
Pn (cos y') 


r 

for — < 1 
r 

r 

for ” < 1 
r 


where y' is the angle between r and r\ 

If Id <.li the limits by means of which the convergence of these 
series is determined are 


\—I 1 +/ 

(1.4.14) a'< a - a'> a -- 

1 ~j“ I 1 V 

Poisson pointed out that there is a layer 

(1.4.15) a - <a < a - 

for which the expressions (1.4.13) are not valid. This question 
was also investigated by Callandreau and later by Hopfner [1]. 
In connection with Liapounov's remarks the so-called desideratum 
of Tisserand [1] (Vol. II, p. 317) should also be mentioned. This 
desideratum was the starting point for investigations of Wavre, 
which will be discussed later. 

In one particular case, namely, that of the Earth's figure, the 
theory has been presented in a form differing from those used 
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by other investigators. The problem posed by Bruns [1] is, 
namely, to investigate the level surfaces as those surfaces which 
are determined by general characteristics of the gravity potential 
and its derivatives. The deformations which actually occurred 
during the history of a planet are, of course, not taken into 
account, since Bruns considers the present distribution of masses 
in the Earth's interior and makes no assumptions about the con¬ 
nection between this distribution and the laws of Hydrostatics 
or some properties of the density. In this respect, the approach 
given by Bruns differs from that formulated by Clairaut. The 
latter identified the Earth's figure with one of figures of equilib¬ 
rium. The problem of Bruns has been discussed in all details 
in the work of Hopfner [1]. 

An attempt has been made by Appell [2] to attack the problem 
of figures of equilibrium by a new method. He made use, namely, 
of the expression of potential in a point of the free surface given 
by Gauss. On writing 



Then, the equation of the free surface takes the form 


i—x 't'j—y a; 

di d')'] dC 

(1.4.19) JJ 9A dJi dX + 

^ ^ +^^52“-const.=0 

d/u dju dfJL 
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On solving this integro-differential equation, one can determine 
the figure of equilibrium of a rotating homogeneous fluid mass. 
One unknown function will be defined by this equation if use 
is made of tangential coordinates. No further development of 
this method seems to be known.^ 

l.S. Conditions that Determine the Figures of Celestial 

Bodies 

The question, of how to explain all those figures of celestial 
bodies which were observed, requires a further specification of 
conditions which might be essential for the formation of a body. 
If we eliminate from the list of celestial bodies all those having 
an irregular form like certain planetoids or nebulae, there will 
be different t37pes of bodies the figures of which can be given a 
satisfactory explanation. In the short survey of methods used 
in the theory of figures of equilibrium given in the preceding 
section, we did not mention aU conditions which must be taken into 
account. These conditions, however, will determine the degree 
of approximation to the actual state of a body, and, for example, 
in the case of the Earth, they may be very important for the 
understanding of the variation of its figure. The theory of figures 
of equilibrium is not able to 5 deld an explanation of many character¬ 
istics in the Earth's structure which must be considered. 

It is evident that the theoretically possible figures are determined 
by a set of conditions which can concern either the physical state 
of the mass of a celestial body, or its kinematical state, or form a 
group connected with the d 5 mamical state in general. Because 
of basically different assumptions which have been made about 
the figures of celestial bodies, several groups of particular problems 
were obtained. 

As to the physical state, the principal condition is that a 
celestial body is or was a fluid body. Most of those observed are 
at the present time in such a state and there are no objections in 


® More references concerning tlie development of methods will be given later. 
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general to the second part of our assumptions concerning the 
past periods. Certain exceptions, like the fact that the Saturn 
ring is composed of small solid particles, as it seems to be, will 
not diminish the importance of results based on the hypothesis 
of fluidity. It is evident also that the change in the physical 
state of a body may be followed by variations in its figure, as it 
may be asserted in case of the Earth. The transition to a plastic 
state or the solidification is a factor which must be taken into 
account. Despite the fact that changes in figures produced by 
these processes are normally small as compared with the size of 
a celestial body, they are very important in case of the 
Earth. 

The second group of conditions refers to the kinematical state. 
In the theory of figures of equilibrium, we had a simple example 
of a postulated state of motion or rest. If a given fluid mass 
rotates about an axis having a direction fixed in space with a 
constant angular velocity, its figure does not vary. The constant 
angular velocity is, however, not a necessary condition for the 
invariance of the figure. Liquid bodies not changing their shape 
can exist if there are internal movements which follow, for 
example, a distribution of velocities given by a law, co = F{s^, z), 
that is, which have a zonal rotation. There may also be other 
important kinematical conditions concerning some regular kind of 
deformation. Postulated pulsations, oscillations of some other 
type, or progressing variations belong to this group of conditions. 

In the third group of conditions, the law of gravitation must 
be mentioned in the first place along with the influence of other 
bodies. We may assume, obviously, any law of gravitation, but 
in all classical investigations the use has been made of Newton's 
law and only this law is here considered. If the influence of other 
celestial bodies is neglected, we have the simple case of an isolated 
fluid mass. Conditions of equilibrium of a system of two or 
more fluid bodies have been, of course, also considered, and in 
more general problems, we have to take into account, besides the 
rotation of each body, its revolution and the variations of both of 
them like those taking place in the solar system. Certain 
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approximations in problems of this kind were possible because of 
very large distances between the members of a system. 

Taking into account the great variety of conditions, a certain 
classification of problems may present an interest and is now 
given. When we consider a figure of a celestial body, its free 
surface is not, in general, the only object of an investigation. 
The stratification is sometimes also of importance and, therefore, 
the term figure will from now on be used in the more general 
sense which includes the stratification. The problems which 
were more or less investigated are as follows: 

A. An isolated mass. 

I. The figures of equilibrium: 

1. of a homogeneous fluid mass, 

2. of a heterogeneous fluid mass, 

3. of a system having a solid (rigid) core and a fluid 
envelope, 

4. of a fluid mass with floating solids, 

5. of a fluid in a sohd crust, 

6. of a plastic mass, 

7. of an elastic solid. 

II. The invariable figures: 

1. of a zonal rotating homogeneous fluid mass, 

2—7, of a zonal rotating mass corresponding to con¬ 
ditions given in L, 2-7, 

8. of a mass under some other conditions. 

III. The vaiying figures: 

1. of a homogeneous fluid mass, 

2. of each system mentioned above. 

B. A system of separated bodies: 

1. a fluid mass and a mass center, 

2. system Earth and Moon, 

3. double-star system, 

4. more complicated systems. 
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Of course, we are able to present only a short selection of 
problems. An attempt is however made to cover the results 
from the whole field. 

The most simple problem on this list is the problem of figures 
of equilibrium which an isolated homogeneous liquid mass can 
have. As mentioned previously, it presented great difficulties, and it 
is not solved yet in all details. The methods of solution applied 
in this problem as well as in the case of an inhomogeneous fluid 
reached such a degree of perfection that certain of them belong 
to the most advanced mathematical tools. However there are 
but few of the problems listed above which are attacked by 
similar methods. It is hoped that more exact mathematical 
methods will be applied in all branches of the theory. 

Two directions in which the problem of equilibrium of a fluid 
mass is developed at the present time should be mentioned here. 
Most of the investigations are dealing with the equilibrium and 
small deformations of gaseous masses. Since stars are such 
masses, their stratification is the main object of the study. Some 
general conclusions reached about the stratification were also 
applied to planets. However, for these celestial bodies and 
especially for the Earth, more attention must be paid to factors 
like increasing viscosity or plasticity and to different periods in 
the solidification. 


1.6. Integral Equation of Liouville 

Since in some of the investigations which will be discussed, 
use has been made of the so-called Liouville formulas (Liouville 
[1], [2]), they will be derived in this section. These formulas 
represent solutions of an equation which belongs to the class of 
homogeneous integral equations. The general form of this 
equation is 


( 1 . 6 . 1 ) 


I 




mC 


where the integral is taken over a certain given surface. The 
element of this surface is dQ' and D is the distance between two 
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points, one of them (M') being a point of the surface. As usual 
C is an unknown function of coordinates of the other point (M) C' 
its value at M'. The factors I' and m will be specified immediately. 
On considering a unit sphere 2 and its element da' we can put 
da' = l'dQ' and integrate in (1.6.1) over the unit sphere. The 
particular case important for the theory is that of an ellipsoid, 
and for this case we will use the notation 


( 1 . 6 . 2 ) 





Tx'da' 

D 


where we have to integrate over the corresponding unit sphere 
or over the surface of this ellipsoid respectively. 

It has been shown by Liouville that the solutions of the equation 
(1.6.2) can be given in terms of Lam6 functions. Since in other 
sections we shall make use of special notations, we introduce 
some of them in the formulae which foUow. The distance D=MM', 
in the expression (1.6.1) or (1.6.2) interpreted as a potential 
will be written in the form D {a, 1). The new parameter a is one 
of those which will be used in order to determine the position of 
a point. The value a' = 1 refers to those points with respect to 
which we integrate. The corresponding complete expressions 
■will be given later. Now, let R, S, M, and N be the Lam6 func¬ 
tions « of the nth. order and R' and S' the values of the first and 
second functions corresponding to a! = 1. Then, the solutions 
of the equation 


(1.6.3) 


%=^ 


J D(a, 1) 


are given by the formulae: 


' The inconvenience that the Lam4 function M has the notation coinciding 
with that of the point M is increased by the fact that the momentum in Section 1 4 
IS denoted by the same letter. It is, however, very difficult to change many 
standard notations to achieve a unification. Since in every method there are 
many concepts involved, it was possible to make only few adjustments. The 
hope can be expressed only that the use of the same letter to denote different 
quantities in different chapters will not produce much confusion. 
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I'M'N'dQ' 

4r7t 

D{a, 1) 

2m+1 

I'M'N'dQ' 

471 

D{a, 1) 

2n+l 

I'M'N'dQ' 

4c7C 

D{a, 1) 

2^-T i 



RSMN iia=l,Mhon the ellipsoid 
RS'MN if ^^<l, Mis an internal point 
R'SMN if 1, Mis an external point 


The first member in every equation (1,6.4) represents the value 
of the potential of an ellipsoid at a point M {a, d', %p), if the density 
of this surface is equal to V M'N', d' and ^ being spherical angles. 
From the integral equation (1.6.3) and its solution (1.6.4) it 
follows that the products MN are the eigenfunctions and the 
expressions 


(1.6.5) 


2w + 1 1 

4:71 RS 


are the eigenvalues. 

For an ellipsoid of revolution, we have 


( 1 . 6 . 6 ) 

or 


m'N'da' 

J D{a, 1) 


4m 

2n + 1 


RSMN 


( 1 . 6 . 6 ') 


I 


y'n.sdo' 


4m, 


D{a, 1) 2w + 1 


RSY„ 


where are spherical harmonics. If we put p, — cos these 
functions are expressed in terms of associated Legendre functions 
j, which themselves depend on the Legendre polynomials P„. 
By definition, it is 


(1.6.7) 
and 

( 1 . 6 . 8 ) 


P 1 - 1)” 

2.4... 2» dp^ 


— PY 


d^PniM) 

dp} 


Then, the spherical harmonics of the nth order are determined 
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as follows: 

(1 6 9 ) Pn,n{ii)cOS'n'ip 

'^n.iKr-x- ■P«.i(/«)sin^,P„_2(^)sin2y,..., „(/<)sin?'Z'VJ 

The eigenvalues in (1.6.6') can also be given in terms of Legendre’s 
functions. It is 

yn.J 9 + /“^ 

y«.s = J [-P«.s(/“) cos Sipfia 

and the variable q will be determined later. A detailed theory 
of Lam^ functions and spherical harmonics can be found, for 
example, in Vol. IV of Appell [1]. 


( 1 . 6 . 10 ) 

where 

( 1 . 6 . 11 ) 



CHAPTER II 


Figures of Equilitrium; Inverse MetLo J 
2 ./. The Fundamental Equation 

A brief account of classical investigations concerning the 
figures of celestial bodies v^ill be given in this chapter. If we 
assume that a celestial body is changing its shape, the problem 
of its motion becomes, in general, too intricate, since in the 
equation of motion the free surface of the body involved is un¬ 
known, and in the case of an inhomogeneous body, even the 
whole stratification is unknown. Thus, the motion will be 
produced by variable forces depending on a varying stratifica¬ 
tion. Even in the simplest case representing the equilibrium of a 
liquid mass, by which the figures of many celestial bodies can be 
explained, the problem has been reduced to the solution of a 
functional equation. New, very complicated, methods were 
required to solve this type of equation, and one can easily see, 
therefore, why the first solutions were found by the inverse 
method. Attempts were made, namely, to investigate whether 
the most simple geometric figures can be free surfaces of a fluid 
mass or not. It had to be proved only that all the conditions of 
the problem will be satisfied by such an assumption. For example, 
a figure of equilibrium of a fluid mass at rest is determined by 
equation (1.1.9). 

. , r 

(2.1.1) —--—= const, 

Jv D 

In this functional equation, the limits of V are unknown and so 
is the density, if the fluid is inhomogeneous. It may be easily 
seen, however, that, on assuming that the figure of a homo- 
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geneous mass is a sphere or that the stratification is represented 
by a set of concentric spheres such that at the surface of each 
of them the density has a constant value, the conditions of 
equilibrium can be satisfied. 

For the figures of equilibrium of a rotating fluid mass we have 
the equation (1.1.8) 


( 2 . 1 . 2 ) 




or for a homogeneous mass 


—— = const. 

2 / 


(2.1.3) 


I. 


dV' 


0 )^ 5 ^ 


const. 


Now, we can see that no sphere satisfies these conditions. But, 
in case of a homogeneous fluid mass, ellipsoids could be among 
the possible figures of equilibrium. Of course, using the inverse 
methods, we can hardly see whether all solutions are found or not. 
Asmentionedpreviously,theexactproofsof thefactthat asphereis a 
unique figure of equilibrium of a nonrotating fluid were given by 
Liapounov and Carleman. For a rotating fluid, however, we know 
many solutions. They represent ellipsoides, ring-shaped figures, 
other figures differing but little from the ellipsoides, and it is 
not yet shown that this list is complete. 


2 . 2 . Maclaurin^s and Jacobrs Ellipsoids 

The basic equations of the theory have been mentioned in 
Section (1.4). We now write the equations (1.4.1)-(1.4.4) of 
that section. If the expression for the potential of a homogeneous 
ellipsoid at a point M {x, y, z) 

(2.2.1) /^ = const 

is inserted into the equation (1.1.8) the conditions for the surface 
of equal pressure are obtained in the form 

(2.2.2) (i, - co^)x^ + (Z„ - a,2)y2 + = C 
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Furthermore, if the surface of an ellipsoid of revolution 

2.2.3 1 

a2 

is identified with the free surface of the fluid 

( 2 . 2 . 2 ') ( L ,- co 2){*2 + 2 /')+ Z ,^2 = C 

the conditions which must be satisfied are: = Ly and 

(2.2.4) ~ co^) = = C 

On deriving further conditions, it is convenient to introduce a new 
variable I determined by the equation 

(2.2.5) a? = c2(l + l^) 

Then, by (2.2.4) and (2.2.5), we obtain 

( 2 . 2 . 6 ) 

For an ellipsoid of revolution having the axes 2a and 2c and 
a constant density x the coefficients = Ly = L and in the 
expression of the potential (1.4.1) are given by the formulae 


(2.2.7) 


L — 27tfxa^cj 


■ 2Tifxa^ c 


0 [a? + A)2(c2 + A)’^“ 
dX 

0 (a^ + X) (c^ + A)"^” 


(see, for example, Appell [1], Vol. IV, p. 48). If a new variable 
(t) is introduced and X = c^t, L and Z, can be expressed in terms 
of elementary functions. It is 


poo 

Z = 27cfx{l + P) 

Jo 


( 2 . 2 . 8 ) 


(1 + i 

1 + ZV 

= — l^ta 


tan “^I 


I+PJ 


1 _L 72 

L, = 47ifx a - tan~i 1) 


The mass of the ellipsoid is also a function of the variable I 
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(2.2.9) 


4 4 

m = — nnci^c — = u— 7^^=-- 

3 3 ^ ^ 3-v/l + 


Now, if the values of m, %, and I are given, we can find from 
this equation a and from (2.2.5) the other axis of the ellipsoid. 
On writing (2.2.4) in the form 

a^L — > 0 


and inserting (2.2.7), we obtain 




XdX 


(a^ + A)2(c2 + X)Vc^ + A 


> 0 


It is obvious that the integral is positive and, therefore, a > c. 
Thus, if an ellipsoid is a figure of equilibrium, it must be an 
oblate one. 

On ehminating now L and from (2.2.4) we obtain the con¬ 
dition 


( 2 . 2 . 10 ) 


co^ (3 -f P) tan~i I — SI 


2}ifx 




^h{l) 


which is to be satisfied by co. If a value of I is given, the cor¬ 
responding CO can be found by (2.2.10). From the last equation 
(2.2.4) the constant C, which will determine the value of potential 
at the free surface, also becomes known. The relationship (2.2.10) 
may be represented graphically. Substituting tan-^ — + 

it is easy to see that the curve h—h{l) which begins at the point 0, 
will reach a maximum and approach as5miptotically to the I 
axis. Some other conclusions concerning <0 were mentioned in 
Chapter I. 

In case of elhpsoids of Jacobi, assuming that c <b < a, the 
coefficients L^, L^, and are given by the formulae (Appell 
[1], p. 64). 


( 2 . 2 . 11 ) 


(*00 


L 271 fKabo j 


(*00 


L^ = 27tfKabcj 


iX 

\a‘^+X)G{X) 

dX 

ic^+X)G{X) 


, Ly=2nfxabc f - — - 

* Jo {b^+X)G{X)’ 

I ^^(A)= (^^-(-A)(&^-f-A) (c^-j-A) 
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Instead of (2.2.3) and (2.2.4) we now have 

/yi2 fl«2 w2 

(2.2.12) ^+l. + L l 

and, by (2.2.2) and (2.2.12) 

(2.2.13) 0)2) ^ j2(£^ __ ^2) _ 

Hence 

a^T 

(2.2.14) 

and 

(2.2.15) 


C 


62 


a^b^L, - L,) + c2(a2 -- h^)L, 
Inserting (2.2.11) in (2.2.14), we obtain 


(2.2.16) 


0)^ 


^TtfK 


= abc\ 
Jo 




{a^ X){h^ X)G{X) 


and by (2.2.15) 


fOO 1 

aH^ 

C2 - 

dA 

Jo 

L(a2 + A)(ft2 +A) 

c2 + A_ 

W) 


0 


This condition is identically satisfied by b = a. Assuming 
that b ^ a the substitution 


^2 ^2 

1=5, 7i = ^ 

62 




yields the expressions 




(2.2.18) 27tfK 
and 

(2.2.19) 0: 


poo 

St 

Jo 


xdr 


h{s, t), 


(1 + st)(1 + tr)A 

= (1 st)(1 tx^ (1 t) 




' rdr r^dr 

ir-'l ipr ') 


The variables 5 and t have positive values according to our 
definition, and the integrals in the last equation are also positive. 
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By (2.2.19) if the given ellipsoid is a figure of equilibrium, it 
must be that s + < 1. 

The mass of the ellipsoid can also be expressed in terms of 
parameters s and t. We have, namely 

4 4 4: b^t 4: s 

(2.2.20) m ~ —Tzxabc = —nx —= — jzx —r: = — nx—r 
3 ^ Vst S Vs ^ Vt 

If the mass, the density, and the values of parameters s and t 
are given, we can compute all axes and the corresponding angular 
velocity, but we have always to take into account that by (2.2.19), 
5 is a function of t. Equations (2.2.18) and (2.2.19) represent 
a curve in the space of variables s, t, h, and the equation g [s, t) = 0 
gives its projection on the plane st. The complete discussion of 
these conditions can be found in Vol. II of Tisserand [1] or 
Appell [1], Vol. IV, pp. 61-69. The function h = h {s, t) reaches 
a maximum at 5 = ^ == == 0.3396 and the corresponding 

ellipsoid belongs to the series of Maclaurin’s ellipsoids. This 
maximum == 0.1871 is the limit value (1.4.9) which determined 
the point of intersection of two curves representing sets of the 
ellipsoids of Jacobi and those of Maclaurin. 


23. Inhomogeneous Figures 

The first investigations of Clairaut concerning the figures of 
equilibrium of an inhomogeneous mass were mentioned in section 
(1.4). His well-known equation representing the connection 
between the average density of a layer and its elhpticity as well 
as some related theorems were for a long time the single base of 
the theory of inhomogeneous figures of equilibrium. To obtain these 
results, we begin, however, with more general assumptions 
concerning the stratification in a mass. We follow the derivation 
given by Wavre [1]. 

If we write the equations (1.1.3) in the form 


(2.3.1) 


1 dj> 
K dx 


X + cd’^x, 


I dp 1 dp 

X oy X ox 


they show that gravity force, i.e., the resultant of attraction and 
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of the centrifugal force has the direction of the normal to the 
surface of equal pressure p [x, y, z) == const. In case >c = cp [p), 
three sets of surfaces: p = const., = const., and the level 
surfaces coincide as it was shown earlier. 

We shall now consider this case and denote the gravity potential 

by 


:2.3.2) 


w 


= 

J X 


+ const, = W [h) 


A stratification is usually taken which corresponds to a density 
distribution such that x can be expressed in terms of a parameter 
which we call d. Its sense will be specified later. 

Then, the gravity potential, as well as the gravity g become 
functions of this parameter W{d) and g{d), respectively. The 
coordinates of the vector g are obviously given by the second 
members in (2.3.1) or by the equations 

dW dW dW 

(2.3.3) g^=ga = —, — — Sz — Sy— 

a, /3, y being the angular coefficients of the interior normal of the 
surface W = const. If dn is the element of this normal, we have 


(2.3.4) 


dn 


^ ^ _i_ ^ #_i_ 

dx dn dy dn 


dg dz 
dz dn 


= n-Vg 


and, by differentiation of (2.3.3), we easily obtain 


(2.3.5) 




[dx dB 
dx dy 


dzl 


The expression in brackets represents the div n if n is a unity 
vector in the direction of the normal. Now, it is 

(2.3.6) div n = —i^ 

where the value of K is the double average curvature of the surface 
W — const. Thus 


( 2 . 3 . 7 ) 


^zw = ^-Kg 
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and by (1.1.6), it is 

co^s^ 

V2 1F = V2/i^+V2 — 

Since the potential fU satisfies the Poisson equation, we obtain 
the equation of Bruns [1] 

(2.3.8) ^^Kg=-4^fx + 20)2 

an 

This is an exact equation which 5 nelds a value of one of the 
quantities involved if the others are known. It may be easily 
seen from the equations (2.3.3)-(2.3.8) that the relationship 

(2.3.8) does not depend on any hypothesis about the connection 
between three sets of surfaces mentioned above. Some other 
conclusions can be drawn from (2.3.8). Suppose that dn and dd 
are both positive at the same time, i.e., we change now the sense 
of dn which was positive inwards and put dnjdd = N. 

Then it is 

(2.3.9) -gN=.W'=~ 

dd 


g being the derivative of W along the normal. Along a surface 
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we can write 
(2.3.10) 

This holds obviously when the corresponding points and 
are displaced along the lines of force and respectively (Figure 1). 
By differentiating (2.3.10) we obtain 


d/Cl J 2 


and from (2.3.8) and (2.3.9) 


(2.3.12) 


dd 


KW + (- 4- 2a)2)iV 


On inserting (2.3.12) into (2.3.11) it is 
(2.3.13) {4^fK - 20)2) (A^l - N\) 

dN\ 


= W' 


ddh 


(kn + 


1 ' 
N ddJz^ 


That this equation can be subject to further useful transfor¬ 
mations has been shown by Wavre, but we restrict ourselves by 
mentioning this formula only. As to the shape of the surfaces of 
equal density, we did not make any assumption as yet. It has 
been proved by Volterra [2] for an infinite number of layers and 
for a finite set of such layers by Hamy in 1887 that these surfaces 
cannot coincide with a set of homothetic ellipsoids. Volterra's 
statement holds when the density is an integrable function. 
Nevertheless, ellipsoids are used as surfaces at least approximately 
representing stratification. Some conclusions drawn from this 
assumption are discussed in the next section. 


2.4. Clairaut’s Equation and Some Theorems 

We shall now assume that (a) W = const, are closed surfaces 
each enclosing all preceding in the set, (b) they have a common 
center, and (c) that they are surfaces of revolution. The equation 
(2.3.13) can be applied to a surface = const. The polar and 
equatorial "radius” of this surface are c,(a), a^{d) respectively. 
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The average curvature at its pole is Kp{d) and at the equator(a). 
Let gp{d) he the gravity at the pole P,. Then we can write (2.3.13) 
in the form 



According to the assumption of Clairaut, the stratification is 
determined by the set of concentric ellipsoids, and the parameter 
d represents their polar semi-axes. Then, the curvatures are 
approximately 


Ke = -, 



where d — a^ — c^is small for ellipsoids differing but little frona 
a sphere. By neglecting the terms of the order of d'^, equation. 
(2.4.1) takes the form 


(2.4.2) 


isifx — 2(0^ 
Sp 


d' 2 6" 

7 "^^^ ~~2 


Furthermore, if cu^ is small, d, d', and d” are also small, the 
product can be neglected and the average density in the 

interior of = const, as well as the gravity can be computed, 
approximately by assuming a spherical distribution of masses. 
Thus we will have 


= 4jr 

Jo 

It follows from this equation that 


4-jt 

3 


r = 


-p = Zclx 
and 

(2.4.3) 1 + = ^ 

c, x^ gp 

since the denominator in the second member ^jifc^x^. represents 
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the gravity at the surface of the sphere. Thus, since = d, we 
obtain by (2.4.2) and (2.4.3) 


(2.4.4) 


2 2 d 1 (3" 

d'^ K d' 2 6' 


The ellipticity e = djd can be introduced as a new variable. Then 
the Clairaut equation takes the form 


(2,4,5) 


a a 


It is also known that for the gravity at a pole Clairaut’s 
theorem can be proved. It is expressed by the conditions 


(2.4.6) 


4 co^d 

— - 

5 g 


^ 2^, 


1 (D^d 

2T“ 


5 

“ — 4 


(o^d 

g 


Two further relationships can be also proved, the second is 
very important in applications. It is, namely, 


(2.4.7) 


^eqna jpole _ ^ CO^U 

g ~ 2 g 


where g is the gravity at some latitude on the surface of a planet. 

All these results are obtained for those stratifications which 
were assumed in this section. We shall also mention here another 
classical relationship known as the equation of D’Alembert 


(2.4.8) 




1 (o'^aX 

I'Y) 


where A and C are moments of inertia of a planet with respect 
to equatorial and polar diameter respectively. 

As mentioned above, when the problem of a figure of an in¬ 
homogeneous planet is reduced to the form given by Clairaut, it 
is assumed that the level surfaces are concentric ellipsoids and 
that the distribution of densities is known. The problem is then 
to determine the variation of eUipticity from the free surface 
toward the center. We mentioned a quite different approach to 
the problem of the figure of the Earth made by Bruns [1]. The 
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basic definitions and assumptions of his theory are as follows. 
The level surfaces represent a set of surfaces perpendicular to the 
force of gravity at each point of the Earth. Because of rotation 
these surfaces are determined hy the function W = fU + Q, 
i.e., by the equation W{x, y, z) = const. Then the gravity 
satisfies the equation (2.3.8). We note first that under certain 
assumptions concerning the mass of the Earth which seem to be 
close to the actual conditions, Pizzetti [1] has proved that the 
level surfaces form a hnear set of closed surfaces each enclosing 
the preceding ones. This fact had been considered at the beginning 
of this section as an assumption. The further investigation of 
level surfaces by the methods of potential theory has been brought 
in connection with geodetic data (see, for example, Hopfner [1]). 

We must restrict ourselves by quoting only few conclusions of 
Bruns. No level surface can be well approximated by a single 
anal 5 d:ical surface having a simple shape. As to the figure of 
the Earth, a level spheroid 

W = const. = TEq = — -b -f 

r 2/ 

should be first considered. It differs but little from an ellipsoid 
of revolution. Next approximation is then the geoid W W-i- 
= Wq. Hence at a point of the geoid the term W should not be 
equal to Wq. Now the expression is the sum of all terms of 
higher orders in the expansion of the attraction potential in a 
series of spherical harmonics. If C is the distance between these 
two figures which has been measured along the geoid’s normal it 
can be proved that it is determined by the equation 



g cos 8 


where g is the gravity at the level spheroid and a the angle 
between the two normals. This is Bruns’ theorem. To find the 
distortion represented by the geoid from the apparent differences 
in gravity g — g is the weU-known problem of Stokes. 
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2,S. Dirichlefs Ellipsoid; Investigations of Riemann et at. 


As the last example of the classical discussion in the theory of 
figures of celestial bodies, we shall take the ellipsoid of Dirichlet. 
Assuming that at an initial instant a homogeneous liquid mass 
has a shape of an ellipsoid 


(2.5.1) 


AC\ AO *^0 

-h ~ H- 

^2 ^ 52 ^ ^2 


Dirichlet considered the motion determined by the condition 
that the coordinates of an element x, y, z are linear functions 
of X, y, z. 

(2.5.2) X ^ lx my -f nz, y = . ,z ^ . 

where Z, m, n are functions of time. As usual the particles are 
subjected to the mutual attraction. The pressure at the surface 
of the liquid can be in general a function of time P{t), Inserting 
expressions (2,5.2) in equation (2.5.1) we see that the outer 
surface of liquid is a variable ellipsoid concentrical to (2.5.1) and 
the equations of motion are satisfied if the pressure 


(2.5.3) 


i> = P{t) + c{t) 




This motion of the liquid can be split into two parts: first, in a 
rigid rotation with respect to an axis passing through the center 
and second, the motion of particles relative to the rotating axes 
of reference The relative velocities have to be propor¬ 

tional to the new coordinates, y, f 

(2.5.4) == 

For the particular case where there is no rotation, Dirichlet 
obtained isochronous oscillations in which the shape of the liquid 
is varying from an elongated to an oblate ellipsoid on taking at 
some intermediate instant the spherical shape. If the initial 
velocity is not equal to zero and does not satisfy some special 
condition, the corresponding ellipsoid becomes indefinitely 
flattened or stretched. 
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If the initial velocity is below a certain limit, the angular 
velocity of the rotating liquid can vary between two limits. 
However, for a larger value of the angular velocity of an elongated 
ellipsoid a sufficiently great external pressure must be applied. 
The existence of ellipsoids of Maclaurin and Jacobi has been 
proved by Dirichlet as a special case of his theory. In an ad¬ 
ditional note, Dedekind derived from the investigation of 
Dirichlet a very interesting result. Each ellipsoid of Jacobi, 
namely, will preserve its shape and position if the internal move¬ 
ments of particles are given by the equations 

(2.5.5) cos o)t+y sin cot, y= —x — sin cot-\-y cos cot, z^z 

0 a 

The constant co is then the angular velocity of the corresponding 
elhpsoid of Jacobi (2.1.18). Each particle, therefore, describes 
an ellipse, the equation of which is given in the parametric form 

(2.5.5) . 

Riemann [1] improved the unfinished investigation of Dirichlet 
and modified the problem. He studied the characteristic changes 
of the axes of the ellipsoid and the relative motion of the liquid 
with respect to them by eliminating the initial moment from the 
equations of motion. Then, it could be proved that there are only 
four particular cases known from previous investigations in which 
the principal axes of ellipsoids do not change. As to the small 
oscillations, there are two types of changes which differ in stability. 
According to Riemann, if the disturbed surface of Maclaurin's 
ellipsoid remains an ellipsoid of revolution, the equilibrium is 
stable. 

In case the circular equator becomes an ellipse and the polar 
axis remains unchanged, there is instability if the eccentricity 
exceeds certain limits. 

The problem of a varying homogeneous hquid ellipsoid was 
further investigated by Brioschi, Lipschitz, Greenhill, Basset, 
Tedone, Love, Stekloff, and Hargreaves. (See, for example, 
H. Lamb Hydrodynamics, p. 722). 



CHAPTER III 


Metliod of Poincare 


3.1. Fundamental Functions 


Poincar^ in his famous memoir [1] suggested a new method to 
solve the problem of figures of equilibrium of a rotating fluid 
mass. This method yields solutions representing also figures of a 
new kind which are usually called the "'figures of Poincar6.'' 
As mentioned previously, in the series of Maclaurin's ellipsoids there 
is one corresponding to the so-called branch point from which 
the series of Jacobi's ellipsoids start. The question arises whether 
there are other branch points of a similar kind. The answer 
given by Poincard was positive. It is necessary to note that 
Liapounov discovered the same figures of equilibrium one year 
earlier than Poincard. He did not decide to insist on their 
existence in his first paper. Nevertheless it seems more ap¬ 
propriate to call these figures "Liapounov-Poincard’s figures of 
equilibrium," on taking also into account that the former in¬ 
vestigator later applied most exact methods to solve several 
problems concerning these and other figures of equilibrium. 

In this chapter we shall discuss the method of Poincard. An 
important part in it is played by the "fundamental functions" 
associated with a certain closed surface S. They are known now 
as eigenfunctions, but the terminology of the theory of integral 
equations was introduced in the investigations of Lichtenstein 
which will be discussed later. 

To see the conditions which fundamental functions must satisfy 
we have to make use of well-known properties of a surface potential 
due to the Newton's gravitation or to Coulomb's electric forces. 
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This potential, namely, is continuous on the surface S, but its 
normal derivative is discontinuous. 

Now, we will determine a set of functions C/j, which have the 
properties as follows: (a) The value in the space bounded 
by 5 is a harmonic function 

(3.1.1) Jf7W = 0 

(b) The value of the function 17*; in space x^ exterior to S 
satisfies also the Laplace equation 

(3.1.2) = 0 

and is regular at infinity, i.e., of the order 0(1/?'). 

(c) At the surface S it is 

(3.1.3) jyU) == jjM 

and, if n is the unit vector of the normal, 

(3.1.4) (grad grad C/W) • n = ^ ^ = 0 

dn dn 

where A*., = 1, 2, 3, .. ., is a set of positive numbers. To deter¬ 

mine these numbers we must have certain conditions. We will 
assume that these conditions are 

(3.1.5) /<») = (grad W*^)^dx is a minimum 

(3.1.6) /<») = (grad W^^^dx = 1 

Then, by the equations (3.1.3), (3.1.6), and Green’s theorem we 
obtain 



IfweputC/^®* = const.,it is by (3.1.6) /{,*’> = 0. 
since /<*> ^ 0. We assume now that hn = 0 
by (3.1.4) 





This is a minimum 
It is, in general. 



hence 
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(3.1.7) 

The functions Uj^ which satisfy all the conditions just given 
are called fundamental functions. The conditions (3.1.7) and 
(3.1.6) can be also written in the form 

Similar conditions represent connections between functions having 
different subscripts, namely, 

r r 

(3.1.9) Ujf^ = 0, dS = 0 k^m 

J s on J s on 

The spherical harmonics yield such functions in case of a sphere 
and the functions of Lam6 in case of an ellipsoid. 

It has been proved by Poincar6 (see also the investigations of 
Stekloff) that in general if a function 0 is determined at a surface 
5 it can be represented under certain conditions by a series of 
fundamental functions 

(3.1.10) 

7c«l 

To find a coefficient Aj, we multiply this equation by 


dn 

and integrate over the surface S. Then, we obtain 

J ' 377 («) 

from (3.1.8) and (3.1.9). This general result has been applied in 
the theory of figures of equihbrium which differ but Uttle from 
ellipsoids, and, therefore, we shall now take into account some 
properties of Lam6 functions. Let q, and v be elliptical 
coordinates and R{,q^), S{q^), and N{v^) the functions of 

Lam6. As is well known, the elliptical coordinates correspond 
to three confocal surfaces passing through a given point P{x, y, z). 
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The equation of a surface of the second order {c < b < a) 


(3.1.12) 


_ 52 + ;i3 _ ^2 


represents ellipsoids, if and > 0 . It will determine 

confocal hyperboloids with one nappe if X^ == and X^ — < 0 , 

;i 2 _ 52 ^ Q a^nd confocal hyperboloids with two nappes, if 
':> X^ = '> 

There are eight points in space corresponding to three given 
values of q, y, v. Let f{Q^) be a polynomial and n ^ 2k its degree. 
There are four classes of Lamd functions: 


(3.1.13) R = 

(3.1.14) R = \/q^- a^fiQ^), - c^f{Q^) 

(3.1.16) R = V{Q^-a^)(Q^-b^)f{Q^). 

V(g 2 _ 52 ) (^2 _ , 2 ) V(p 2 _c 2 )(g 2 _« 2 ) /(g 2 ) 

(3.1.16) R = -v/(^2 - a^)[Q^ - b^) (^2 - c2) /(g2) 

The functions M and iV are represented by similar expressions 
where or is substituted for p^, for example, by (3.1.13) it is 
M = iV = /(v^). 

Expressions such as MN or RMN are called Lam 6 products. 
They are solutions of the Laplace equation. A Lamd function 
of the second kind S is defined as the second solution of the 
differential equation 

d^R 

(3.1.17) = [n{n + l)pu + h]R 


where pu is the Weierstrass function. The constant kin (3.1.17) 
is determined by the condition that solutions are of the type 
considered above. The function 5 must also satisfy the condition 

^dS ^dR ^ ^r^2n + l 

(3.1.18) f?__S_= 2 »+l. 5 = i? —^du 

du du Jo R^ 

For more details about the Lame functions we refer to 
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Vol. IV, 1, Chapter VI of Appell [1]. We mention here only few 
particular solutions and values of Lam^ products. These parti¬ 
cular solutions are as follows: 

(3.1.13') n = 0, i?o == 1 

(3.1.14') « = 1, R^ = Rs = Vq^-c^ 

(3.1.15') « = 2, R,= R,R„ R, = R,R^, R^ = R^ R^, 

R 7 = s^ — [g+iV 3 ^ 2 ], Rg = Q^-[g- ^Vsg^] 
if we write the Weierstrass function in the form 


and put 


pu 


-I 


ds 


°° V4s® — e,s 


u 


gs 


3g2 = Kg + H = h, K = n{n + 1), 
= pu-i- g 


where Zf is a constant. It is obtained when the equation (3.1.17) 
is transformed. Using the variable q this equation takes the form 


(3.1.17') 

where 


R dQX dq) 


= Kq^-^H 


J.... (0^ - a^)(e^ - b^)(g^ - c^) 

Forn = 1, there are three functions of Lam6 and three products 
which are, denoting by C an arbitrary constant, 

(3.1.19) RMN=Cx,Cy,Cz 

For n = 2, there are five R functions and, therefore, five Lam6 
products. Three of them are 

(3.1.20) RMN = Cxy, Cyz, Czx 

For M = 3, we have seven Lamd products. By (3.1.16) one of 
these expressions is 

(3.1.21) 


RMN = Cxyz 
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A new variable I wiU be used later. It is defined as a derivative 
of the vanable u introduced in (3.1.17) along the outer normal 

(3.1.22) ~ = 12^ 1 

’ (> - Q^){v^ - Q^) 

Now, assume that a certain ellipsoid is determined by putting 
Q Qo and that is a harmonic function on this elhpsoid. We 
can put according to Poincare 

A«=0 

we^take^^^ the superscript 0 the value of a function for q = q^, 


(3.1.24) 

Then, the function 




(3.1.25) 


A )=0 


^ the tomonic function in the interior of E. having on 
this ehpsoid the value U,. On the other hand 

(3.1.26) V-ia.ffiMS.Jlf.JV. 

^=0 

is^a^armonic function in the space exterior to and U, = 

Let us assume now that 

(3.1.27) ^ C/» = C/(^) = 

It is obvious that these functions satisfy the condition i/W = 
on ijg. It will also be * * 

if we put 

(3.1.29) Aj, = _ yW j .^(0) 

\3o/o/ *= laoL 


(3.1.29) 


A —A 
dn 


swI 
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Thus, Uj^ are fundamental functions for because of the 
properties of Lam6 functions. 


3.2. Figures of Equilibrium That Differ But Little Prom 

Ellipsoids 


We shall now find the potential of a homogeneous ellipsoidal 
layer in terms of Lam6 functions. Assume that the varying 
thickness of it is t and its constant density is « = 1. The mass 
element can be interpreted as a product of da by the surface 
density C- Then, we have the well-known condition for the 
discontinuity of the force at the simple layer 


(3.2.1) 


dU, 

dUf 


+ 


dU, 

dn^ 


= - 4:71^ 


On taking into account the expressions (3.1.25), (3.1.26) for 
the functions 17^ and we can write their normal derivatives 
in the form 


dn^ du dfii 0 ^ du ^ dn, 

oo 


dU, 


0 (f'U 


on^ du on o dn^ o du 


where 4 is the value of /I at g == go and u the parameter which 
is introduced in the Lam6 differential equation (3.1.17). 

Thus, by (3.2.1) we obtain from the expressions for the normal 
derivatives of U^ and U^ that 


(3.2.2) 


OO 


— 

0 




By (3.1.18) it is 


00 9^ -L 1 oo 

(3.2.3) f-= 

0 4jc 0 

where 


2^ -f- 1 
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00 


0 


is a given function, the coefficients are known. Then, we have 
also and can compute coefficients of Lam6 products in (3.2.3) 
and vice versa. 

If the coefficients C in (3.1.19) are determined to give the 
connection between cartesian and elliptic coordinates, the former 
are expressed in terms of latter as follows 

(3.2.6) x = h^R^M^N^, y = z = h^R^M^N^ 

where 


(3.2.7) h\ = 




K 




[c^ - a^){c^ - 62 ) 

The volume of an ellipsoid can he also given in terms of functions 
It is, namely. 


(3.2.8) F = I ie, i?, i?3 = ^ ^ ^ ^ 

Now the direction of the normal to an ellipsoid E is determined 
by the cosines of the angles as follows 

(3.2.9) cos {nx) = ^ = hJM^N^R^, cos {ny) = h^lM^N^R^, 


cos {nz) = h^lM^N^R^ 

If Ue is the potential of a given elhpsoid E and U' that of a 
lapr on this ellipsoid having the thickness f, the sum Ue + U' 
will represent the potential of the deformed ellipsoid. For an 
incompressible fluid, the thickness of this layer will be subject 
to the condition 


(3-2.10) f = Q 

JS 

In order to express the potential of an ellipsoid in terms of Lame 
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functions, Poincard derived first the following formulae for the 
case where the ellipsoid is displaced in the a;-direction (and in a 
similar manner for other directions). It is by (3.2.9) 

^ = s cos {nx) = sh^l^RlM^N^ 

If the thickness of the layer is reduced to one term in (3.2.3) 
its potential is also represented by a single term (^ = 1) 

(3.2.11) — 


Now, the change in potential of the ellipsoid due to the displace¬ 
ment e is obviously 

BUb 


(3.2.12) 


U, 


E 2 




dx 


since in the second position the same particle has the coordinate 
X — e with respect to point P where the value of the potential 
is considered. By (3.2.11), (3.2.12), and (3.2.6) the first of the 
equations (3.2.13) is readily obtained, the other two follow from 
similar considerations concerning the displacements in the y- or 
^-direction. Then, each of them can be written in the form: 

dU^ 


(3.2.13) 


dU. 


By 


c(0) 

V ^ X 
c(0) 

. V-^-v 


djh 

dz"' 




It is obvious now that at a point in the interior of the ellipsoid 
the potential is 

For the external points it may be shown that 
bUb dUB 


(3.2.14) 


U, 


Ei 


y*2 J-^ /iji 


s'o* \ 


(3.2.15) 




rJ’ 


i?3 


dx By R^'^’ dz 

If we now take into account that, because of notations used in 
this section, the a;-axis corresponds to the small axis of the 
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ellipsoid and that, therefore, in the equation (1.1.8) x- and 2-axis 
must be interchanged, we shall write and insert in 

this equation the expression (3.2.14) for the potential. Then, 
we obtain, by putting / = 1, 

(3.2.16) - 7 ^ 0 : 2 +const. ==o. 

The equation of the ellipsoid = ?(,) can be also written 
in the form 


(3.2.17) 


t{z^, y, — ^(0)2 + —1 = 0 


by (3.1.12) and by the definition of the functions (3.1.13)- 
(3.1.15') given earlier. If this ellipsoid is a figure of equilibrium, 
the equations (3.2.16) and (3.2.17) represent the same surface 
and we obtain the conditions (on omitting the subscript 0) 


S^R^-S^R^ 

V R\ 

S^Ri — _ Sgifg — 

If the volume of the hquid F and the angular velocity co are 
given, these two conditions joined to (3.2.8) will determine the 
length of axes of the ellipsoid since R and S are functions of a, b, c. 
In case of the ellipsoids of Maclaurin it is R^ = R^ and we have 
only two equations, namely (3.2.18) and (3.2.8). The equation 
(3.2.19) which is to be used for the ellipsoids of Jacobi, can also be 
written in the form 


(3.2.18) 
and 

(3.2.19) 


{S^R^ - R^S,)R^R, = - Rl) 

and subject to an important transformation. By (3.1.18) we 
have for « = 1 
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Since it is by (3.1.13)-(3.1.15') == 52 __ ^2 

i? 2^3 = ^4 and by (3.1.18) 


S^ = R 


l*U 

Jo 


Rl 


du 


we obtain the condition 


(3.2.20) 


RiS^ R,s^ 
3 5 


It may be shown that this is equivalent to (3.1.19). 

Now we shall discuss the existence of figures of equilibrium other 
than these ellipsoids. On assuming that the new figures differ 
but little from an ellipsoid of Maclaurin or Jacobi we can represent 
such a figure (£') as mentioned above as an ellipsoid E covered 
by a layer having a varying thickness f. This may obviously have 
positive as well as negative values. The condition 

^0 + ( 2 ^ + •^o) = const. = /7o 

holds for each point P of E, if this ellipsoid is a figure of equilibrium. 
Therefore, at a point P' in the neighbourhood of P we can write 


(3.2.21) 

and, if we put 


(3.2.22) 


n. 


„ dn 


dn 

dn 


calling g the gravity acceleration, we have 
( 3 . 2 . 21 ') n = n,-g!: 

Now, U' is the potential of the ellipsoidal layer and if E' is 
a figure of equilibrium we have the condition 

77 + C7' = iTo - gf + 17' = const, 
on this surface, or 

(3.2.23) U' -g^= const. 

The thickness f which will determine the deformation or the 
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departure of the new figure from the ellipsoid can be represented, 
by the series 

00 

(3-2.3) 

0 


If we insert this series in the equation (3.2.12) we obtain 

0 JS 

Because of the properties of Lain6 functions (see, for example^ 
Appell [1], IV, p. 159) all these integrals are equal to zero 
for k> 0. Thus, the coefficient /Sq must also vanish and we have 

(3-2.24) 

1 

By (3.2.3)-(3.2.5) we have, then. 


(3.2.26) 

1 2m -)- 1 

The value of the gravity acceleration at the pole (y = 0 z = 0) 
is by (3.2.22), (3.2.14), and (3.2.17) 

c(0) 

(3.2.26) 

and 


(3.2.26') 






By substituting (3.2.24), (3.2.25), and (3.2.26') 
this condition takes the form 


in 


(3.2.23) 


(3.2.27) 


OO 

1 


(SkRk 

S,RA 

\2w + 1 ' 

3 / 


= const. 


In this appreximation g is considered as a constant. In the 

r^d a coordinates 

fi and V. Therefore, it must be 
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(3.2.28) 


2w + 1 


5 , 


k — 1, 2, 3, . . ., const. = 0 


These conditions are basic in the theory of Poincar^. They 
were also given without the factors /3j. in the first paper of Lia- 
pounov. The factors are coefficients in the series representing 
the function f (3.2.24). They can differ from zero, because of 

(3.2.28) , only if the corresponding second factor vanishes. Taking 
this fact into account, we have: (a) if w == 1 and ^ = 1 the 
coefficient can differ from zero, since the second factor in 

(3.2.28) vanishes identically. This case does not present any 

interest, the value of giving just the displacement of the 
whole ellipsoid in the a:-direction. (b) In two cases n = 1 and 
A — 2 or A = 3, the coefficients must vanish because of (3.2.18) 
and (3.2.19). (c) If n = 2 and ^ = 4, /S 4 can differ from zero, 

since the condition (3.2.20) holds for ellipsoids of Jacobi, but this 
case is not important either, the value of just giving a pure 
rotation of the ellipsoid. Thus, new figures can be expected for 
k = 5 or for larger values of it. In order to answer the question 
whether there are such new figures of equilibrium Poincard 
investigated the more general expression 




R,S,. 


(3.2.29) F 

' 2# + 1 2j!. 4- 1 

where i ^ k and p can be equal to n, and the function 


F 

' W 


1 


(ir 


(3.2.30) F^- + 1 2p FI R^ 

This function depends on the parameter (F varying in the inter¬ 
val (00, a^) and the problem is reduced to that of the existence 
of zeros of F-^. The detailed discussion of the function (3.2.30) 
shows that for m = 2 and /fe = 5, 6 , or 7, i = ;^ = 1 there are 
no zeros and, therefore, no new figures of equilibrium. In case 
k = 3 we have the branch point for Maclaurin’s and Jacobi's 
ellipsoids. 

The first new figures of equihbrium are obtained for n = 3 
and n = 4:. In general, in order to find such new figures, i.e.. 



54 


FIGURES OF CELESTIAL BODIES 


the nonvanishing coefficients and, therefore C ^ 0, use must be 
made of Lame functions of a new kind. These functions can be 
written in the form 


(3.2.31) M), Ve^-b^fie^), VQ^-c^fie^), V{e^-b^){e^~c^)f(Q2) 
where /{q^) is an entire rational function of degree k. Its zeros 
(a, a,) are between and a^. Let oc be the largest root. Then 
we can put 

(3.2.32) /(^s) = (g2 — cc){q^ — a,) . . . 

The number n can be equal to 2^ or 2^ -f 1. We shall write only 
the equations for two figures of the lowest order. 

If w = 3, we have the function = V— c^{q^ — a), where 
c® < a < b^- Thus, we can write the expression 

C = = IobR^M.N^ 

The factor s is equal to and i?f has a constant value 

on the ellipsoid q = q^. The expressions ilfj and are obtained 
by writting jj, and v instead of q in Ry.. Therefore, it is 

^ = /osV(5^-c2)(,a2_c2)(„2_,2) (p2 _ (^2 _ (^2 _ 

Now, by (3.2.6) this square root is proportional to z and 


(e® — a) — a) {v^ 



both expressions having the same roots because of (3.1.12). 
Therefore, the equation of the new figure of equilibrium takes 
the form 


(3.2.33) 


j=ez 

‘‘0 


-J_ ^ 1 ^ 

\a — 

OC — 0^ 

a — 


where s' = const. This is the so-called ovoidal or pear-shaped 
figure. 

In the case « = 4, we put R^ = /(g^) = (g 2 _ (^2 _ 

where < ag, oci < It is now 


C -1,8(q^- ai) (fi^ - aj {v^ - aj (g^ ~ a^) - x,) (v^ - a^) 
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and this equation can be written in the form 

(8.2.34) f = 

L \oL-i ~ a. — 0^ oci -- 


The constant a' can be either positive or negative, and, therefore, 
two new figures are determined by this equation. 

These examples and the general analysis show that there are in 
general in the expression of ^ as many factors of the form 


as there are zeros of the function Ry,. Thus, it is, in general, 

(3.2.35) !; == Iqs' if n = 2k 

f == Iqs'zE^E^ n === 2k 1 

The conditions here discussed hold for the ellipsoids of Jacobi. 
The conditions of the existence of new figures of equilibrium which 



Figure 2 
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differ but little from the ellipsoids of Maclaurin take a different 
form. We must have first -f-A = 0 (mod. 2), i.e., n and k 
must be both either even or odd numbers. Putting 

(3.2.36) — = tan w, y -r = cos'& = ft 

y 

we can make use of the fact that the Lam4 products in case 

of ellipsoids of revolution, determine the spherical harmonics. 
Taking these functions in the form (1.6.9) we can write 

(3.2.37) f = loWkPn, 7 c ifi) cos ky> + 4, {p.) sin kv) 

If we put /Sj., = jSj, tan w tMs equation takes the form 

(3.2.38) Y = Pn.k iP) cos {kip — co) 

cos CO 

It may be easily seen that the axes of coordinates can be rotated 
to make the angle S == 0. Then, (3.2.38) takes the form 

(3.2.39) C — Pn.kifi) 


N 



Figure 3 
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and for 0 < A < w it will determine a set of new figures of equi¬ 
librium differing but little from the ellipsoids of Maclaurin. 

In two extreme cases /e = 0 and k = n we have 

(3.2.40) C = W.P.(,) 

and 

(3.2.41) f ifl) cos^v. P.,n = (1 - 

Equation (3.2.40) represents a ''zonar' figure of equilibrium. 
The function C does not vanish identically at the equator or at 
the poles. A typical figure of this kind is shown in Figure 2. 
Equation (3.2.41) gives the so-called sectorial figures: C vanishing 
at poles and along a set of meridians (Figure 3). If 0 < A < n 
we obtain a network of meridians and parallels dividing the 
surface of the ellipsoid into cells where C is alternating from 
positive to negative values. 

3.3. Stability of Figures of Equilibrium 

Since the theory of figures of equilibrium also had the purpose 
of explaining the evolution of celestial bodies, the stability of 
different figures had to be determined. This question, however, 
is one of the most difficult in the whole theory. It requires some¬ 
times, as it was the case of the pear-shaped figure, a very intricate 
analysis and the highest precision is necessary in the notions and 
criteria used. 

Forces of three types can be considered according to the general 
discussion of the equilibrium given by Thomson and Tait [1]. 
First, there are forces usually taken into account like universal 
attraction which depend only on the coordinates of particles. They 
represent a system of conservative forces. Now, if a liquid mass 
is in equilibrium, the corresponding equations are the same 
for an ideal or a viscous liquid. However, if it is assumed that 
a viscous mass is deformed, the new state differing (maybe even 
little) from an equilibrium will in general give rise to the second 
type of forces which depend on velocities such as friction'. In 
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this case, the system is a dissipative one. For most of the figures 
that are considered, the equilibrium is defined with respect to 
a set of moving axes, and the third type of forces is represented 
by the centrifugal and Coriolis force. For these forces Thomson 
and Tait used the name gyroscopical forces. The basic statements 
concerning the influence of introduction into conditions of equi¬ 
librium of a new type of forces are as follows. 

An unstable equilibrium produced by conservative forces alone 
can be transformed under certain conditions into a stable one if 
the gyroscopical forces will be also introduced. 

The kind of equilibrium which exists under the action of con¬ 
servative forces alone will not be changed if the forces of both 
other types will be introduced. 

If the equilibrium is stable when forces of all three types are 
acting, the removal of dissipative forces as well as of the gyro¬ 
scopical ones does not affect the type of equilibrium; similarly the 
instability will also persist. 

These conclusions follow from equations of motion of a system 
of particles, written, for example, in the Lagrange form and from 
the fundamental condition 

(3.3.1) ±(^r-U)=:-f<0 

In this condition 2T represents, as usual, the kinetic energy of 
the system, U is the potential of conservative forces, i.e., F = grad 
U and 2f = biiql + , . . where the coefficients are 

not negative is the dissipation function. The condition (3.3.1) 
shows that the total energy of the material system is dissipated 
if the function / differs from zero. 

If in the expression for /, all coordinates are involved the dis¬ 
sipative forces are complete. According to Lord Kelvin the 
stability is a secular one if it takes place under the action of 
complete dissipative forces. It is the temporary stability if 
there are incomplete dissipative forces acting. For further devel- 
lopment of the theory, certain criteria of stability were introduced 
in mathematical form. 
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It is well known that if the equations of motion of a mass system 
are written in the Lagrange form 


A (^\ ^ 

dt \dqj dq^ dq^ 

the positions of equilibrium are obtained from the equations 


(3.3.2) 


(3.3.3) 


0 


i = 1, 2, . . .,n 


The classical investigations have proved that U reaching its 
maximum is the necessary and sufficient condition of a stable 
equilibrium. It is, of course, a new assumption that the results 
proved for a discrete mass system hold for a continuous body as 
well. Since it is made in the theory of figures of equilibrium of 
fluid masses, this assumption enables us to make an important con¬ 
clusion from the equation (1.1.8). 


(3.3.4) 


r it'dV' 

]v D "^2/ 


const. 


If we multiply the expression representing the left-hand member 
by fniV and integrate over the whole volume of the fluid we obtain 


(3.3.5) 


- / 




dV 0)2 

Yj 


s2 xiV 


where W is the energy of the system due to attraction and I the 
moment of inertia with respect to the axis of rotation. It is easy 
to see that in the case of equilibrium, this expression is constant, 
and, therefore, its derivatives with respect to the coordinates are 
equal to zero. In general, in small movements which will follow, 
the distortion of a figure of equilibrium forces of all three types will 
be involved. Lord Kelvin’s criterion of stability is, then, that 

(3.3.6) 

is a maximum. 

In this expression the angular velocity is assumed to have a 
given value. There are, however, cases in which a deformation of a 
figure of equilibrium will not be performed under the action of 
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figure ofPoincar^ indicated a uniform contraction of a 
ODinion J ^ example of such a phenomenon. In his 

case is more^^'^T^ total momentum is constant in such a 

being M = 7' a ® ~ const. This momentum 

stabilitv TV the validity of another criterion of 

(3.3.7) _t a 

2 1 

(3.3.6). 

a varying paramJef^ Thr”^^ quantity can be considered as 
of equilibrium r^PT^ '''^® obtain a continuous set of figures 

eqmhbnum depending on this parameter. This representation is 

Si 


- 


V', 

i ^ 


Figure 4 

ttr ®®ution here only 

he dia^am corresponding to the distribution of stable and in^ 
stable figures of equihbrium. In this diagram (Figure T 
meter is used which is proportional to the angulaf momentum S 
the hqmd mass and s = c^/a^ where c and a are semi-ax^ 
he set of Maclaurin’s ellipsoids is represented by the line 
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It starts at a point which corresponds to a sphere and on the 
branch we have stable Maclaurin'sellipsoids. These become 

unstable for values of parameters giving the part At the 

branch point the stability is transferred to Jacobi's ellipsoids 
the set of which is represented by the line J /g. As mentioned 
in Section (2.3), the new figures of equilibrium differing but little 
from ellipsoids are determined by the zeros of expressions (2.3.29) 
or (2.3.30). Therefore, the next branch points on each of the lines 
MqM-iM or // 1/2 will be determined by the roots of equations 


(3.3.8) 
and 

(3.3.9) 


2^ + 1 3 


2 ? = Ro 


or 


2?4 5 ^ 




== 0 


the former holding for M-^M, the latter for J / 1 / 2 . 

At each branch point, when moving from the first one [M-^) the 
degree of instability of the ellipsoids belonging to the same set is 
increasing. It is obvious that all these considerations present a 
very high interest, especially from the viewpoint of Rational 
Mechanics, but cannot be applied without restrictions to the 
evolution of figures of planets or stars, since the assumption that 
celestial bodies are homogeneous masses is an oversimplification. 


3.4. The Ovoidal Figure of Equilibrium 


Among the figures of equilibrium of a homogeneous liquid mass 
differing but little from an ellipsoid there is one which attracted a 
special attention of investigators. It is given by the equation 
(3.2.33) 


(3.4.1) 


s'z 

\<x- 


, + 


x-b^ 


+ 


and called a pear-shaped or ovoidal figure. If 


(3.4.2) 


a:“ z 

_i_ j_ _ 

e 
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is the critical ellipsoid of Jacobi, by deformation of which the 
ovoidal figure can be obtained and x, y, z are coordinates of a point 
at this surface, we can write 

(3.4.3) C" = {x-xf+{y-yf+ 

If {xld^)jl^,... are the cosines of angles which the normal to (3.4.2) 
makes with the axes of reference, we have 

(3.4.4) £=a;+— 

The intersection of the surface (3.4.1) with the Jacobi ellipsoid is 
indicated in Figure 5. It was pointed out by P. Humbert [1] that 
the meridians of the ovoidal figure do not have points of inflexion. 



Therefore, this figure is actually more hke an egg than a pear. 
The representation given by Poincard [2] has been based on the 
first approximation, and this suggested the name pear-shaped. 
Certain speculations concerning the fission of a celestial mass have 
been made because of this representation. In a theory of forma¬ 
tion of the Moon developed by Darwin [3], [4], p. 317, the possibil¬ 
ity of such a process was connected with the question of stability of 
the pear-shaped figure. If a figure having such shape were stable, 
an evolution of it could be expected such that the narrow belt 
would contract, and the figure would be divided into two separate 
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masses. The first investigations of Liapounov have shown, how¬ 
ever, that the ovoidal figure is unstable and, therefore, the separa¬ 
tion of the Moon from the Earth should not be explained by the 
fact of the existence of such a figure. The question of stability of 
the ovoidal figure has been considered as an open one for a longer 
time (see, Appell [1]). Upon developing a new method of approach, 
Jeans [1], [3] confirmed Liapounov's conclusion. 

In a recent discussion of the question of stability Lyttleton [1] 
corrected several other statements of Jeans concerning the fission 
of a rotating fluid mass. He reached the conclusion that the fission 
process as applied by Jeans to the origin of binary stars is also 
impossible. 

3.S. Some Other Figures of Equilibrium 

To explain the existence of the ring of Saturn and the shape of 
some nebulae, the ring-shaped figures of equilibrium were in¬ 
vestigated. A critical analysis of the results concerning the figures 
of this kind was given by Poincar6 [2]. Although the study of a 
ring-shaped figure of equilibrium of a liquid mass was initiated by 
Laplace, an exact proof of the existence of such a figure has been 
given much later by Sophie Kovalevsky [1]. However, neither the 
hypothesis of a liquid ring, nor that of a solid can explain the ring 
of Saturn as it was shown in investigations of Laplace, Maxwell, 
and Poincard. It was mentioned before that only the hypothesis of 
Cassini (that the ring of Saturn is composed of small solid particles) 
can meet all requirements of observations and theory. 

Some other figures of equilibrium were obtained in the two- 
dimensional problem. These figures investigated by Matthiessen 
in 1859, Jeans [1], and Globa-Mikhailenko [1] are infinite cylinders 
with an elliptical cross section. It was pointed out by Liapounov 
that such problems are nothing but mathematical curiosities. 
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Mettod of Liap ounov 

4.1. Potential of a Body in Which the Stratification Differs 
But Little Prom an Ellipsoidal One 

The solution of a problem concerning the figures of equilibrium 
of a fluid mass depends on the series expansion of the potential of 
this mass. As mentioned previously, Clairaut’s theory could not be 
extended over the first approximation. Although the methods of 
Laplace and Legendre yielded any desired approximation in a formal 
way, the convergence of the series used remained an open question. 
Neither Poisson nor Callandreau could make essential progress 
investigating it. Even the degree of precision which could be 
obtained by the method of Poincar^ (see Chapter III) was not 
quite satisfactory as pointed out by Liapounov. In a series of 
papers he developed, therefore, a new method of attacking the 
problem. After more than 30 years of investigation in this field, 
the final form of the method was given by Liapounov in his last 
work [9] published 1926, seven years after the tragic death of its 
author. This form of the method which can be apphed to a hetero¬ 
geneous body as well as to a homogeneous one will be called here 
the method of Liapounov. 

It must be pointed out that the conclusions are not restricted by 
the a.s.sumption that the departures from the ellipsoidal structure 
which will be taken into account, are infinitely small. They can 
also have finite values, and in this respect as well the method of 
Liapounov is superior to other methods. 

Now, in order to write the expression for the potential of a body 
the stratification in which differs but little from that which is 
given by a set of surfaces of ellipsoids, we sh all m ostly use the 

notations of Liapounov. Let t>e an 
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ellipsoid that belongs either to the set of Maclaurin's ellipsoids 
i^q == 1 ) or to that of ellipsoids of Jacobi [q < I)- Then 
£^(^fVp+l, aVaVq) represents a set of ellipsoids similar 
and concentric to E. Suppose that in a given body B the surfaces 
of equal density differ but little from the ellipsoids E^, 
the density is decreasing from the center of the body towards its 
free surface. We assume that the variation of the density is small 
and that the law of densities is given in the form 

(4.1.1) X = l-\-d(p{a) 

where 9 ? is a given function of a parameter a. This parameter will 
determine the surfaces of equal density {x = const.); (5 is a second 
parameter having small values. The case of a homogeneous body 
is obtained from the formulae which follow on putting (5 = 0 , 

X = const. Obviously the parameter a can be given different 
interpretations. We shall assume that a is determined by the 
condition that the volume of the space bounded by the level 
surface a = is equal to that of the ellipsoid . In the center of 
the body it is (^ = 0, and at the free surface F of B we have <3^ = 1. 

If (p{0) is zero the density at the center is = 1 , i.e., the value 
adopted by Em den [ 1 ] for a standard solution. The form (4.1.1) 
was chosen by Liapounov to make easier the application of the 
series expansion obtained for a homogeneous ellipsoid in the prob¬ 
lem concerning heterogeneous mass as wiU be seen later (equation 
(4.3.4)). In general, however, we could assume that 

(4.1.1') X — (p{a, d) 

and, on expanding this function in a power series with respect to the 
parameter (3 obtain 

(4.1.1") X = 1 + '^8^(pj,{a) 

The results of Liapounov's investigation can be easily general¬ 
ized for this law of densities on adding terms corresponding to the 
sum 

&=2 

Now, according to Liapounov, the level surfaces can be 



METHOD OF LIAPOUNOV 


67 


represented by the equations 

X — a{l + Q + 1 sin d' cos %p 

(4.1.2) y^a{l + q sin '& sin ip 

z == a{l + C)a/ q cos 

where ^ and y are spherical coordinates. The function C must 
have small values according to assumptions just made. In general, 
the formulae (4.1.2) will determine the position of any point in 
space ifO^a^ oo, 0 '^'d' n, 0 ^-^^ 2jr. To each function 
?9', y)) corresponds a set of surfaces On the free surface F 

it is 


(4.1.3) f(l, yj) = C(^, V’) = C 

If we assume that to each pair of values of the angles & and y) 
corresponds a single value of f and that the differential coefficient 
dC/da does exist, the equations (4.1.2) yield the expression for a 
volume element as follows 


(4.1.4) = aVe(e + ?)(e + l)(l+C)2(l+ sin 'd' da d'& dy) 


We put now 

(4.1.5) A - Vi{Q + q) {Q + .1) da = sin dd' dy) 

Since according to the definition of the parameter a the volumes- 
and must be equal we have 


= -—a^J (1+0®«?(T = Ijra® Ve (§+?)({?+ 1) 


By (4.1.5) it follows from this equation that J(1 + = 4:;r 

and, therefore, we obtain the first condition for the function C in 
the form 

(4.1.6) -jc^da-ijc^da 

Let fU be the potential at a point M{x,y,z). At M'{x',y',z')i 
we have 
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u = 


Cx'dV' 

D 


where D = MM'. On substituting (4,1.4) and (4.1.5) we obtain 


(4.1.7) U = A 


y! da' j* 

0 


{i + cr{i + 


da‘ 


9 


da' 


D 


In order to solve the problem it will be necessary to expand this 
potential in a series. Then, the function C must be subjected to 
certain conditions. Liapounov assumes first that this function 
can be expanded in a power series with respect to the parameter <3. 
We shall see that in some other problems, use can be made of cer¬ 
tain other parameters as well. Let l{d) and g{d) be two variable 
numbers that do not depend on a, a!, •&, d'’, %>, f' and vanish for 
d = 0. We put 

(4.1.8) (D)j_f/=o = Dq 

Then two other conditions for C introduced by Liapounov are 


(4.1.9) 


la < 


2 ^ 


Since C is not a given function, the existence of two such num¬ 
bers should be proved and this has been done by Liapounov after 
a formal evaluation of 

For further transformations it is convenient to take a new 
variable f instead of C' using the definition 


(4.1.10) 

Hence 


i + c 




(4.1.11) 


1 + C' 

1 + C' 


1 + f (i + c)-Mi+ 


da' 




1 + 


da' ^ 
da' 


By (4.1.2) we obtain 
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V(a:—. . 

V(tH-Ty[7(T+F) sin & cos y)~a'{l+^') sin cos ipT+- ■ ■ 

(1+f) V (e+lT[" sin ■& cos y)—a'{l+i) sin &' cos v']^+(e+?)[- • •?+?[• • 

If we put a'(l H- I) = 6 and D{a, b) for the second factor in 
(4.1.12) and use new notations we obtain 

(4.1.12') £>==(! + ^)D(a, b) = D{a + af, a' + a’C) = £>[?, C'] 
D{a, a') — D[0, 0] = Dq 


Expressed in terms of the new variable | the potential U takes 
the form 

/ da'i\ 

(4.1.13) Jo J 


da' 


da' 


D{a, b) 


41(1 +f) 


f/ da’b^da' 

Jo +17)i)(«T) 


To expand now the surface integral here involved in a power 
series we will take into account that 


r b^da' 

J 'dJu^T) 


oo ^fn 

= 2 — 
nfo W! 


13* 

mi^da'\ 

19&"J 

D{a, b)i 


since b = a' a' and that 


J D{a, b) 


00 a '2 


^ r 

'6«J 


b^i' 


da'I 
da' 


'o«! \db^J D[a,b) 


da' 


It is 

a'"r da'i _ 1 d{a' ^)"+^ 

wi 9«' (^+ 1 ')- 


and, therefore, the surface integral in (4.1.13) is 
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J D{a,h) JD{a,a 
” d fa'"!” b^da' 


D{a, a') laa"-! db J ^Ml=a' 


-L'Vi U M/U I 

„ti la^ D{a, 

^ [• ^ V J_ ^ 

JD(a,a') m! da'\db^-iJ~D(a7^]^^, 

On substituting the variable r again this expression takes the 
form 

r = f-?^ + _ ^ 

J D{a,a'y l + o D{a.a') 

+ y ^ J_ _Lr«"’(r - !:)HUa'\ 

„=2 (1 + ^)” da'Xdb^-^] 3{a, b) 

Thus, (4.1.13) is transformed into equation 


U = (l+f)M f 

Jo J D(a,a' 


Jo JD{a,a') 

Jo aa'j D{a,a') 

+A y ii+girp... I j- p-K'-O-iV,- ! 

"=2 «! Jo aa'laa^-ij Z)(a, S) 

and we can put ** " 

(4.1.15) C7 = ^ ^ ^ 

The terms of this^ series are now given by the expressions 

U, = A\\'^^'da' 

Jo JD{a,a’) 

= 2C:fo+^ \\’da’~ 

(4.1.16) -D(«. a') 

^2 = Uyu r da' — r , 

•)o da'J D{a,a') 

+4r«'*'i- i-KgifW)!., a 

2 Jo aa'laftj z)(a, a) 
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The values of these terms at the ellipsoid’s center (a 

f da' 

3 

U,{0) 


0 ) are 


(4.1.17) 


fi a ra'H'da' 

i %' da' — I -rrT-—r 
Jo da J D{0, a ) 


For « > 2 it is F„(0) = 0. On making use of the method of 
Cauchy (majorantes functions),Liapounov was able to “at 

(4.1.15) is an absolutely and uniformly convergent series, if th 
conditions (4.1.9) are satisfied and if 

first term m (4.1.15) V, is the potential 

•« wTiiob the ellipsoids £„ are surfaces of equal density. Th 
m which the eiiipsoias , ^ 

be easily seen from (4.1.2) ana 

UndeT Lditions (4.1.9) and (4.1.18) the etpiession 

u' -ty 

can be also expanded in » 

serica, as it has been shown by JeTt «d 

to note that in all “r^Xe^j^t developed will show how 

to he infinitely smal. reduced to a set of integral 

tilt} functional equation (1-1- ) 
and integro-differential equations. 

4.2. Integral Equation of Liapounov 
-me tether development oi the theory is based on the solntion 
t)f the integral equation 

(4.2.1.) P “ 

wliicll has been given by a given contin- 
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on ■& and f. D is, as before, the distance between the points 
¥{■&, ip) and f') of the ellipsoid E, and the integration is 

over the unit sphere Z". Putting in (4.1.12) ^ = 0, a = a' = I 

we obtain 

(4.2.2) i)2 (g-j-i)[sin cos ip— sin cos ip']^ 

+ (i?+?)[sin ^ sin sin sin y']2+^[cos ^9—cos 

Instead of this expression, the distance Z> (a, 1) can be also used 
in the equation (4.2.1) which is obtained from (4.1.12) and (4.1.12') 
if C = r = 0 and a' = 1. 

In order to solve the inhomogeneous hnear-integral equation 
(4.2.1) use has to be made of LiouviUe’s homogeneous equation 
(Section (1.6)). 

If the function S is continuous f will also be continuous. Now, 
equation (4.2.1) will have a solution if the function will satisfy cer¬ 
tain conditions. To obtain these conditions we multiply (4.2.1) by 
and integrate the product over the unit sphere. Then we 

obtain 

Reversing the integration in the second term it is 



Since the expression (4.2.2) is symmetrical with respect to the 
variables, the formulae (1.6.6)-(1.6.6) yield 

Hence 

From (1.6.5) and (1.6.10) it follows that the expression 
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depends on q. Then, we have 
(4.2.4) =0 if v = 

i.e., for certain values of g only. Ellipsoids which correspond to 
these values are the ellipsoids of bifurcation. For the ellipsoids of 
Maclaurin there are conditions 

(4.2.6) Ti,o = 0 Vi,o = V 

and 

(4.2.6) ^v,2?c = ® '^v,2k-l = 0 

Another set of conditions is required for the ellipsoids of Jacobi. 
Humbert [1] and Orlov [2] calculated the elements of ellipsoids of 
bifurcation or critical ellipsoids up to those of the seventh order. 

Now, taking into account (4.2.3) and the equation preceding it 
we can easily see that the function S must satisfy the conditions 

(4.2.7) = ^ 

for all values of n and s such that the equations (4.2.4) hold. It is 
always 

(4.2.8) 

We assume now that S can be expanded in a Laplace series 

(4.2.9) ~ 

which converges absolutely and uniformly. Then, it is 

(4.2.10) = 
where 

(4.2.11) 

since the functions y„,« are orthogonal. A particular solution of 
(4.2.1) can be readily found. It is 


(4.2.12) 
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where the summation is extended over all values of n and s except 
those for which it is = 0. Substituting now (4.2.9) and 
(4.2.12) into equation (4.2.1) we obtain 


or 



y 


n,8 


This is an identity. 

Thus, the general solution can be obtained if the terms corre¬ 
sponding to the conditions ^ = 0 are added to the particular 
solution (4.2.12). These terms are composed of the spherical 
functions Y^ ,.. It is, then. 


(4.2.13) 


- Y 

w, r m,r 


n, s 

where « are arbitrary constants. Their number does not exceed 

three The expression (4.2.13) is a general solution. This is 

proved by the fact that a continuous function (f) is completely 

determined if the values of all inteerals f AV ri • ^ 

vu , integrals ] are given. 

can h. "^^^^^ned by the coefficients , or 

can be found if we assume that the constants , have some 

Z n l ‘‘>rT *“ “Editions 

fuetaM'aUqrtiOT*^n “tdilions the 

theory of homogeneous figures of eatuKhr°™°'''i”°f“‘^’ “ 
little from an ellipsoid. 

at the free surface by the SS ® 
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X = + C + 1 sin # cos y) 

y =z V ^ + C + 1 sin # sin ^ 

^ == f cos ^9’ 

instead of (4.1.2) and assuming that C can be expanded in a series 
2 Ci Liapounov obtained integral equations of the form (4.2.1) for 
the functions Ci- The integral equation of Liapounov (4.2.1) has 
been also investigated by Crudeli [1]. He could show that it can 
be reduced to a Fredholm equation of the second kind having a 
finite kernel.’^ 


4,3. Transformation of the Fundamental Equation 


By writting the equations of a level surface in a parametrical 
form we shall have small values of the function C(^, '9', w) if we 
assume that the stratification in a figure of equilibrium differs but 
little from the ellipsoidal one. In each case, however, the function 
can be expressed in terms of coordinates x,y,z, i.e., 

= y, z). 

The equation of level surfaces is, then, by (4.1.2) 

(4.3.1) —^ + _ = 

Q + I Q + q Q 

and the free surface will be given by 


(4.3.2) 




+ -^ + —= + 
Q + i Q 


where C* is the value of C sit this surface and a = L 
Since the density is a function of the parameter a and x = 
the pressure f can also be expressed in terms of this parameter, 
p = p{a). Therefore, in the equation (1.1.5) the right-hand mem¬ 
ber is a function of a and the equation of level surfaces which fol¬ 
lows from equations of motion of the liquid is 


^ To our knowledge Crudeli is the only investigator who called the solutions 
of this equation Liapounov’s figures of equilibrium. 



76 


FIGXJKES OF CELESTIAL BODIES 


(4.3.3) U + ^ix^ + y^) = funct. [d) 

The potential U can be expanded in series (4.1.15). Thus, the 
left-hand member of (4.3.3) becomes a function of C This equa¬ 
tion is now the fundamental equation from which the function C 
can be found. If the law of densities is given in the form (4.1.1) a 
transformation 

(4.3.4) U=^^+d0 

•will introduce two new expressions W and 0 which can be derived 
from U. To obtain these functions the factors 1 or (p(a') must be 
taken in (4.1.7) instead of Then, if we expand the function 
IF or <5 in a series 

OO 00 

(4.3.5) 

<=o j=0 

the functions and 0{ are determined by the formulae (4.3.16) 
where, obviously must be put equal to 1 or to f{a'). The first 
term represents the potential of a homogeneous ellipsoid E. 
Now, in the second equation (4.1.16), if = 1 the integration with 
respect to a' yields 

(4.3.6) = 

For further transformation we restrict onrselyes by an ellipsoid of 
revolution. Let M{x, y,z) be a point of an ellipsoid It is 

(4 3 7) ^ ~ ^ + 1 sin ^ cos y = aV ^ -f 1 sin sin y) 

z = a\/~Q cos A = (q + 1) y^Q 

It is easy to see that the variable I in (2.1.5) can be expressed in 
terms of ^ if one takes the axes a ^/^ + 1 and ( 1 \/q used in (4.3.7). 
It is / = 1/ a/q. On writing now the constant in (2.2.1) in the form 

2KfA tan“^-^ 

Ve 

and using the equations (4.3.1) and (2.2.8) the function ??ois given 
by the equation 
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Wo^27tA 


tan~ 


(4.3.8) 

We put now 

(*.3.9) 

The integral 


Vq 

— cos^ 


(^ +1) sin^11an“^ 




Ve 


dt 1 

--;;-- = tan-’- —p 

p it \)-\/t ■\/Q 


I 


_ r 


1 ) 


Vq \ 

Q + V 


represents the potential of a layer on the ellipsoid E at an interior 
point M. It is equal to its value at the center of E where D (a, 1) 
= D{0, 1). By formula (4.1.12) we have in case of an ellipsoid 

D{0, 1) = (g H- sin^i?')^ 

On putting cos # = /i we obtain 

f _ r si n -&'d{^'dxp' 

J D{0, 1) Jo Jo y/g -{■ sin^i9' 

, 1 

== 271 _: = 471 sm-’- ,1,=-= 

•i-iVe + 1 -Ve + 1 

Since 


sin-^ - = tan"’- 

Ve +1 Vq 

and by (4.3.9) / = 47rC it is 

(4.3.10) = (2Fo - 47rCd)f + d J 

Because of (4.3.4) and (4.3.5) we can write the equation (4.3.3) in 
the form 

CO® °° 

(4.3.11) IPo + !Pi + — (a:® + y®) = funct. {a) — (5cP — 2 

2/ 1=2 

Since £ is a figure of equilibrium, the condition 



78 


FIGURES OF CELESTIAL BODIES 




^,+~[xl + yl)=^K 

must be satisfied on each surface the sum in parentheses being 
the square of the distance of a point on from the axis of rota¬ 
tion, and ha is a number not varying on this surface. It depends 
however, on a. The expression for h^ can be determined from 
(4.3.8) by putting ?? = 0. It will be 


ha = ^TtA 

If we put now 

(4.3.12) 

we obtain 

(4.3.13) 

By (4.3.7) it is 


tan~^ —p 

Vq 
2? = Ve 


(Ve-etan->i)_ 


Q tan~^ —p 

Vq 


^ta = 27tA{C - a^R) 


*a + 2/a = + 1) sin^# 

and the condition for the surface Ea takes the form 

.2 


(4.3.14) 


0)^ 


+~7^"(^ + 1) sin2^ = 27tA{C — a^R) 


The expression for Wq determined by this equation and Wi from 
(4.3.10) can be inserted into (4.3.11). Then, we have 




27zA{C-~a^R) - sin^'d'- 

2f 


co^ 


4:71 A a^R-i ~—(^ +1) sin^ '& 


where 


0)2 r t'da' ~ 

+-57 J 0(^1) = ^ ,5 (») 

s2 = ^2(1 _|_ ^)2(g 1) 


since the condition (4.3.3) holds for level surfaces Fa which differ 
from Ea- Rearranging the terms we obtain 
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(4.3.15) 


r Cda' 
J D{a, 1) 


27iA{a^R-C) 


a)“ 


— funct. (a) 

2/ «=2 


This equation is of the type 

1 r I'da' 

- 1 - 


(4.3.16) 


4;ra^J D[a, 1) 


W + P(a) 


We can put, for a ^ 0, 


W = 


(0“ 


(4.3.17) 


8/jrZl 


(g+1) sin^i9'C^ 


+ 


ijta^A li=2 


2 W,-W^{O) + [0-0{O)]6 


(4.3.18) = — {a)-\-27cA[a^R—C) — 'Fi{^)—0[^)^} 

The terms ^^^(O) = and <5(0) = [0]„=o do not depend on 

■& and ^p. As to the function W it depends on and terms of a 
higher order. 

On taking into account all expressions for W,Wi, and 0 we can 
easily see that equation (4.3.16) belongs to the class of integro- 
differential equations. The fundamental equation (4.3.3) is now 
written in the form (4.3.16) for figures of equihbrium which satisfy 

the conditions mentioned above. 

The solution of this equation C is a function of two variables « 
and ■& since it is written for an ellipsoid of revolution, and of a 
parameter 5. The general case where f = v) was also 

solved by Liapounov. There is a second unknown function in 
the equation (4.3.16). It is denoted by the symbol funct. («). 
Nevertheless it is possible to determine the function C since we 
have the second condition. This is the equation (4.1.6). 

Now, to solve the problem we have, according to Liapounov, 
assume that the function C is expanded in a senes 

i=l 


(4.3.19) 
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where the term Co is omitted to eliminate the case of a homogeneous 
hquid. If we put 

OO 

(4.3.20) 

i=i 


equation (4.3.16) yields the set of conditions 


(4.3.21) 




1 r Ty_d(T' 
4jra^J Z)(«, 1) 




for the coefficients C«- 

Besides these equations, substituting (4.3.19) into (4.1.6) we 
obtain the conditions 


Cxd<r = 0 U,dc=^- Ulda 

(4.3.22) . ^ , 

j^sda = — 2 jCiCa^o- — .. .,J = iV,- 

where depends on Ci, C21 • • In order to eliminate the 

functions in (4,3.21) we multiply this equation by do and inte¬ 
grate over the unit sphere 2*. Then, 


(4.3.23) RNi 


47ia^J 




In the surface integral , the denominator D{a, 1) is the distance 
of a point on to some point on E. If we put f f' = 0 in 
(4.1.4), (4.1.5), (4.1.7), the last-named equation will represent the 
potential of a homogeneous layer between the surfaces and 
Ea~hda ^ point of £. The density of this layer is equal to unity, 
and we obtain the expression 

r do" 

a^daA [q) J J, 

On changing the notations in the well-known expression for the 
potential of a homogeneons ellipsoid (see, for example, Appell [1] 
Vol. Ill) at a point M {x, y, z) outside it we obtain this potential in 
the form 



t+1 t + q~T)j(r) 


(4.3.24) U, = nA {q) 
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where e is a positive root of the equation 


(4.3.26) 


t + 1 t + q T 


which corresponds to the ellipsoid confocal to and passing 
through M. Now, the potential of the layer just considered will 
be 

/ 7 A f J c z^\ de 

DU = 27tA (Q)ada —tiA (^) U^--77T 

" Je A(t) £+1 s + q s/A(8) 


{Q)ada^ 


— £+1 B + q 8 JA{b) 

since and £ are functions of the parameter a. The expression 
in brackets vanishes because of (4.3.25), and we obtain comparing 
two expressions for the potential of the layer 

r do 2jr r°° dt 
(4.3.26) ZW 

For an ellipsoid of revolution (?= 1), (4.3.26) can be expressed 
in terms of elementary functions. It is 


(4.3.27) 


431 1 

— tan~^ —ji 
a, ■y s 


where e is a positive root of the equation 

(4.3.28) t Z sin^ +-^cos^&' = 

if “b 1 i 

This equation is readily obtained from (4.3.25) and (4.1.2) if 
those equations are written for the case f = 0. 

Eliminating P, by (4.3.26), (4.3.23), and (4.3.21) we obtain 


(4.3.29) 


jL fei-j_1 

TicP J ^ xDia, 1) 2ciJe d (i)/ 


= W, 


Wider + — RNi 
43r 


The right-hand member of this equation depends on functions 
Cl, • • •, Ci-v If functions are determined (4.3.29) can be used 
to find the next function Ci- 

Let us put now in (4.3.29) a = 1 and denote by a bar the cor- 
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responding value of a function of this parameter. For Ci we have, 
then, the equation 


(4.3.30) 





Wi + const. 


This is an integral equation of the type (4.2.1), and we have just 
followed the method used by Liapounov to solve it. For the 
ellipsoids of Jacohi the expression will depend on both angles, 
# and hut we will consider an ellipsoid of revolution only. The 
formula (4.3.17) shows that is a function of The last two 
terms in (4.3.29) yield constants. If now a function li is de~ 
termined from (4.3.30), the corresponding function Ci iBay be 
directly computed by (4.3.29). Thus, the solution of the integro- 
differential equation (4.3.16) is reduced to the solution of an infinite 
set of integral equations (4.3.30). The terms of the series (4.3.19) 
can he computed successively from this set. The effective com¬ 
putation of the first three terms in case of a heterogeneous liquid 
required, however, several hundreds of pages in the quoted me¬ 
moir of Liapounov. It should be noted that in order to exclude 
any doubt about the existence of these figures of equilibrium, 
Liapounov gave a proof for every intermediate step. The im¬ 
portant progress brought by this method consists of the possibility 
of using any desired approximation. On the contrary, as mention¬ 
ed before, Poincare’s analysis was restricted to the first approxima¬ 
tion only.® 

We can give here only the final expressions of Liapounov for 
the first approximation, i.e., the function Cv It is like all 
represented by a regular Laplace series 


(4.3.31) Cl =2 

«=1 ( 8 ) 

(Liapounov [9], p. 170). The symbol S denotes the summation, 
■with respect to the index 1. 


® A simplified computation in case of an ellipsoid of revolution was given by 
Jardetzky [3] for the problem of a zonal rotation. 
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The first step to obtain this result is to expand the function 
in a series 

oo 

(4.3.32) Wj_ = 'X S + funct. [a) 

n=l {1) 

Avhere 

cp{av){l-v-y-^^dv 
^n,oW 


,io.™ < -.J" 


and by (4.2.11) 


y2n,4Z ■ 




The coefficients denoted by an array are by definition 


(2n, 2l\ 

(ai, 2/) = - 


{4n -f l)d 


y2e,4j®2z,4i^2w,4l 


where 


®’2ri,4z(^0 •^2?z,4z(^'V^^^) 

+ 1 

^ ® 2 w, 4 z(^^)J^ 


(4.3.35) ^ ^ 2n 

S2w,4z(^) 

and 


y Y 

-*■ 4j 2n,4Z j 

-H— 




+ 1)(^ + ?) 


(4.3.36) H = {Q+p^)(e+v^) = {e-{-q)esin^^cos^'P 

siii2 5?sin^^+(e + l)(?i-?) cos^^ 

® and % are old notations for Lain6 functions (see, for example, 
Heine, E., Theorie der Kugelfunctionen, 1878, Vol. I, p. 358, 385). 
It is E{/u) = M, E{v) = N, etc. in notations used in Chapter HI 
of this took. 

Now, the function i3 used in (4.3.33) is a solution of the differ¬ 
ential equation for a hypergeometric series. It is 

(4.3.37) = ^[n 4- 1 — «, + w + |, ■+ f, .z) 


Finally it is 

, f1 \ 4^7 

(4.3.38) 

where T„,z depend on eigenvalues of the integral equation as 
given by (4.2.3). 
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These few formulae show the extreme complexity of Liapounov's 
theory. It was apphed by Liapounov in his memoir to the case of 
a homogeneous liquid, and the existence of figures of equilibriura 
which differ but little from an ellipsoid was confirmed once more. 

The last conclusion of the theory is as follows (Liapounov [9] 
p. 436): in all cases (determined by the conditions mentioned 
above) the surface of a figure of equilibrium is determined by the 
equation of the form 


oTi 


+ 


y2 

e + ? 


2 % oo 

+ —=1+ 2^.(X, Z)a^* 

q 4=1 


(A) 


where is an entire function of X, Z of the degree {m — 2)i q- 2. 
It is an even function of Z and an even or odd function of X whe'ri 
the product mi is even or odd respectively. As to a in this equation it 
is a parameter in terms of which the function C could be expressed. 
In the investigations concerning the homogeneous liquids, 
Liapounov made use of the expression a = \/J^da and also of 
the second parameter ri = (co® - (o\)l2f. 

The fact proved by Liapounov that are all even in Z yields 
immediately the symmetry of every figure of equilibrium witli 
respect to the equatorial plane. The general form of this state¬ 
ment was suggested by Lichtenstein (Chapter V). 

The importance of the last result of Liapounov's theory is 
evident. Equation (A) represents the most general known solu¬ 
tion of the problem on figures of equilibrium. 


4.4, Clsiruut^s EquBtion aad M.ore Qeaeral EquHtioas 

The formulae quoted in the preceding section can be easily 
transformed for the case of spheroids. On omitting the factors 

+ Ij Vg q, Vg in (4.1.2) we obtciin level surfaces which, 
differ but little from spheres having radii equal to a. We -^^ban 
denote the radius of the free surface of an undisturbed hetero¬ 
geneous mass by A. Let a: = >c{a) be the density and co the 
angular velocity of rotation having a small value. The formulae 
used by Liapounov (1903, 1904) to derive the equation of Clairaut 
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are particular forms of (4.1.7), (4.1.15), (4.1.16), and (4.3.3). 
There is no more small parameter d in the law of densities, and the 
function f representing the departure of a level surface from a 
sphere was expanded by Liapounov in the series 

(4.4.1) C = Cl or + + • • • 


where 

(4.4.2) 


coM" 

fm 


It is evident that this is the ratio 


A 


fm 

'~a' 


where w == 4:rr is the total mass. A set of equations for 

coefficients Ci resulting from (4.3.30) can be resolved if we assume 
that 

where P, are Legendre’s polynomials and p, = cos The func¬ 
tions must satisfy, then, Legendre-Laplace equations 


IJa 


2j + l* 


Clairaut’s equation (Section 2.4) is obtained for i = 1, f - 2, 
if we put 

S7t x„ w = — 

(4.4.4) 15 (7 ' ^ 15A® 


where a„ is ellipticity of the fl-layer. From (4.4.3) the empticity 
can be determined in terms of a. On differentiating (4.4.3) twice 
the Clairaut differential equation is readily obtained. To solv 
these equations a law of densities must obviously be given. 

The solution of Clauraut’s problem is now represented by 
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(4.4.5) C — ^ ti2‘^2(/^) 

In order to extend this theory to those cases where the density is 
not a continuous function of the parameter a, Liapounov con¬ 
sidered also equations of the form 

(4.4.6) aW{a) 

for all values of a in the interval (0, ^). The symbol 5 in this 
equation denotes a Stiltjes integral, and, therefore, (4.4.6) is an 
extension of the known integro-differential equations in which the 
Riemann integrals are involved. 



CHAPTER V 


Liclitenstein’s Investigations 

5.1. Homogeneous Mass 

The complexicty and extreme length of Liapounov’s theory are 
due to the difficulties of the problem of equilibrium of a liquid 
mass and, of course, to the desire not to leave out the proof of any 
point of the theory. A modification of the method of Liapounov 
has been given by Lichtenstein [2] in a series of papers. His work: 
Gleichgewichtsfigure^i rotierender Flussigkeiten [1] represents the 
synthesis of these investigations. Having discovered an exact 
method for the solution of a functional equation (4.3.3) which 
reduces it to an integro-differential one as well as to a set of 
integral equations in the special case considered in Chap. IV, 
Liapounov did not present his method in a more general form. 
He did not even make use of these names for the types of equations 
just mentioned. A generalization has been found by Lichtenstein. 
Adding some new statements, Lichtenstein was able to discuss the 
general problem of figures of equilibrium in the vicinity of a known 
one. Poincard postulated the existence of new figures of equilib¬ 
rium which differ but little from any given figure of equilibrium of 
a homogeneous liquid mass. Making no assumption about the 
shape of this primary figure, Lichtenstein derived the fundamental 
integro-differential equation for such a general problem, and this 
equation is of the type considered by Liapounov [1] (Chap. IV). 
The method of successive approximations and Fredholm’s theory 
provided the base for the solution of this equation. All these 
generalizations, however, concern in the first place the mathema¬ 
tical treatment of the problem, and no new types of figures of 

equilibrium are actually computed. 

According to Lichtenstein, it is convenient to make use of Gauss 
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parameters i and rj and to write the equation of the free surface 
of a homogeneous liquid mass in the form 

(5.1.1) x=^F{^,rj), z = H{^,7j) 

This surface is not restricted to be a single boundary of a contin¬ 
uous mass. It can have several separated parts, but it is assumed 
that 


(5.1.2) 


rd{F, G)Y r9(G,i?)T 


- d[H, F)Y 
- - 


0 


i.e., there are no singular points on it. 

Let S now be the surface of a figure of equilibrium of a homo¬ 
geneous liquid mass corresponding to the angular velocity o. 
We assume that for a new value of angular velocity coi, differing 
but little from co we obtain in the neighborhood of S another sur¬ 
face Si- The equation of the new free surface of liquid can be 
written in the form 


(5.1.3) C = 

where t is a small distance of a point on Si from the surface S 
measured along the normal to the latter. We assume that 
If I < Ho, where Sq is a certain small number. The equations of the 
surface Si are 


(5.1.1') z = HS>ri) 

and, if we write conditions of equilibrium (1.1.5) for S and 5^ and 
subtract, we obtain the equation of the form 


(5.1.4) 


U^{F^,G^.H^)-U{F. G,H) 

2 2 


If the pressure at the boimdary surface is zero and the volume of 
the liquid does not vary, the term C' vanishes. Lichtenstein 
assumes that now: 


|fl. 


dC 


9: 

di 

i 

d7] 


^ e* < 8(, 
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These conditions are equivalent to Liapounov’s inequalities 
(4.1.9).^ The fact that the latter hold has been proved by Lia- 
pounov, as mentioned above. If rj, are the values of para¬ 
meters corresponding to any given point P{x, y, z), the potential 
U and Ui respectively can be expressed in terms of these para¬ 
meters and we have 

(5.1.5) U{x, y, z) - W{1 rj, C*), U^[x, y, z) - rj, f*) 

Then, a unit mass on S will be subjected to the force of attraction 

(5.1.6) f^W{§,ri,0) 


and the gravity acceleration will be 



~VJ 


U y = dsjdn is the cosine of the angle which the perpendicular 
to the axis of rotation makes with the normal to S, we obtain 
by (5.1.5) 

(5.1.7) g= ca^sy = fi 


Let D denote the distance between two points: (|, tj) and (!', r]'). 
According to Lichtenstein we can assume that the left-hand mem¬ 
ber in (5.1.4) can be expanded in a series 


(5.1.8) f{U^-U) = .., 


{7(«)=/f —K^da' 

‘is D” 


where is a polynomial of the nth degree in C'. dCI3t, 
dC'Idrj'. Then, 




® It is not even necessary to postulate the existence of according 

the last investigation of Lichtenstein [4], When ^ is determined it can be shown 
that these derivatives exist. 
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are proved to be convergent series as well as (5.1.8) even for com¬ 
plex values of It can also be shown that if is the density 

(5.1.10) ?7(i) =/f 

on Js D 

If we now put^® 

(5.1.11) == 2/A, = 

and insert the expressions (5.1.7), (5.1.8), (5.1.10), (5.1.11) in 
(5.1.4) we obtain the equation of the form 

(5.1.12) = C'~sH+ W{X, f) 


Since aF hG — sy the function 


co‘ 


(5.1.13) 


{a^+b^)i^-2syX^~ {a^+b^)XC^ 
- j (C/<2> + C7(3) + . ..) 


It represents the sum of terms of the second and higher order with 
respect to the variables C and A, while the coefficients | and K'jD 
of C and f' are finite. Because of (5.1.8) the equation (6.1.12) is an 
integro-differential one and similar to (4.3.16) which has been 
derived by Liapounov for a particular kind of figures of equilibrium 
namely, for ellipsoids. If we take a certain definite value of 
there is a set of values of the parameters ip, ef which under con¬ 
dition 0 ^ e ^ 1 corresponds to a set of bodies considered by 
Liapounov. To rj, dC, when 0 ^ ^ 1, it corresponds to a set 
of bodies between a given figure of equilibrium and a new one in 
Lichtenstein’s theory. 

In the equation (5.1.12) g has a negative value. As in Liapou- 
nov’s theory, the solution of (5.1.12) is based on the discussion of 
the homogeneous linear integral equation 


10 This parameter A is identical to tj introduced by Liapounov and mentioned 
at the end of Section 4.3. 
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(5.1.14) 



ri(i' = o 


This can be easily transformed into an integral equation with a 
symmetric kernel by the substitution u — ^V—gh. There are 
two trivial solutions of this equation. They correspond to an 
evident property of each figure of equilibrium. When it is trans¬ 
ferred, namely, along the axis of rotation or it is rotated by some 
angle with respect to the axis, it remains a figure of equihbrium 
corresponding to the same values of potential and angular velocity. 
These two solutions are u-^ = const, and = const. {Fb — Ga) 
and correspond to the conditions Ti q = 0, ^ = 0 of Lia- 

pounov’s theory. (See (4.2.4).) 

As is well known from the theory of linear-integral equations, 
the eigenfunctions [u^) of the transformed equation (5.1.14) can 
be normalized and orthogonalized. Then, it is 

(6.1.15) = 0, for j # I, f^xgufda = —1 

If 1 ? and ip are spherical polar coordinates and we have a figure 
of rotation, the eigenfunctions are either of the type F{^) or 
Fid) coskip or F{&) sin kip. 

If we put D{P, Pi) for distance between points P and P^ and 

1 m 

(6.1.16) -jj = N{P, Pi) + 


the equation 

(5.1.17) + = 0 

has no solutions representing eigenfunctions as it has been proved 
by the theory of integral equations. 

On the other hand (6.1.12) can be written m the form 

m 

(5.1.18) K +J = C' — s^X 'Fix, 0 "1“ 


where 

(5.1.19) 
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Lichtenstein proved that this equation has a simple solution if A, 
C, Afj,.. have sufficiently small values. This solution as 
well as its first derivatives with respect to the parameters f and r] 
are continuous functions and satisfy the conditions (6.1.4). Itcan 
be evaluated by successive approximations. Let Q[P, Py) be the 
resolvent. Then (6.1.18) takes the form 


(5.1.20) 


1 mi 

C = - (C'-s2A) + I M,Uy+ - F(A, c) 

5 i 


r n m 1 - 


da' 


Since (see, for example, Courant and Hilbert [1]) 

( 6 . 1 . 21 ) j^Q{P,P^)u'yda' = 0 


on writing 


(6.1.22) C ~ s^A + W{X, C) = W(C', A, f) 

(6.1.20) takes the form 


1 m /* 1 

(6.1.23) C=-iy(C',A,f)+2^»«.- Q(P,Pi)^W(C',A,C')da' 

6 1=1 J S ^ 


Now, as in Liapounov’s theory for a special kind of figures of 
equilibrium (Section (4.3)), the successive approximations are 
obtained from a set of integral equations. The first two of these 
equations are obtained from (5.1.18) by writing 


/K +1 — C — s^A + ^iMyUi 

(5.1.24) 

= C' - S^A + + 1F(A, Cl) 


the second and other equations follow by an obvious transforma¬ 
tion. 

The equations (6.1.24) are nonhomogeneous Hnear-intearal 
equations of the form 
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m 

(5.1.25) + x'NC*'da' — C ~ -f- + W{X, C*) 

where is a given function. 

According to the general theory of equations of this type, the 
solution of (5.1.25) is of the form 

(6.1.26) C* = P + 0(A,^) 
and in case of (5.1.24) we can put 

1 r 1 

(6.1.27) P=--(C'-s21)+ Q{P.P^)\;{C'-s’n)da' 

S !=1 Js g 

On taking into account (5.1.20), (5.1.21) we can easily see that 
by (5.1.26)—(5.1.27) the solutions of (5.1.24) are given as follows: 

C^ = p 

(6.1.28) C, = P + e{X,Q 
fa = P + 0{X, fa) 


where 0 is the sum of second order and higher terms. It is by 
(5.1.18) 

(5.1.29) |0 = -j^x'N'0'da' + W{1 f) 

The equation (5.1.18) is equivalent to the fundamental equation 
of the problem (5.1.12) if the conditions (5.1.19) are satisfied. 
On writing these conditions in the form 

(5.1.30) Mi +j^Pi'i'u'iC'da'= 0, 1=1,2, ...,m 

it can be easily seen that Afi = Mg = 0, because of the vanishing 
of Ui and other values of subscripts these conditions are 

obviously analogous to the conditions (4.2.4) of Liapounov’s 
theory, i.e., correspond in general to some figures of bifurcation. 
Thus, a new set of figures of equilibrium can be expected in the 
vicinity of a given figure if this figure of a homogeneous liquid 
mass belongs to a special set determined by the conditions of 
bifurcation. 




94 


FIGURES OF CELESTIAL BODIES 


5,2. Heterogeneous Mass 

In Chapter IV of his work, Lichtenstein considered some new 
heterogeneous figures of equilibrium in the vicinity of a given 
figure which is not necessarily homogeneous. These investigations 
have a general character and have never been carried out to some 
actual evaluations of new figures like those determined by Lia- 
pounov’s theory. The density law used by Lichtenstein is of the 
form = h{x) + (3;^:(a) i.e., similar to (4.1.1). It is also assumed, 
that the parameter 6 takes small values, i.e., this is the case of 
the so-called “weakly inhomogeneous” figures. When the first 
term of the density is a certain function of a parameter a the 
primary figure of equilibrium, i.e., that in the vicinity of which, 
new figures are sought, will be a heterogeneous one. Moreover, the 
density must not be continuous in the whole liquid. If there are 
some surfaces of discontinuity, Stiltjes’ integrals instead of the 
Riemann integral must be used. As mentioned above this has 
been done by Liapounov in the generalization of Clairaut's problem. 

Now, we have to show the principal steps in the derivation of the 
fundamental equation of the theory. 

Let us assume that all surfaces of equal density intersect the 
.?-axis. Then, the distance q of the point of intersection from the 
origing of coordinates can be taken as a parameter which will 
determine these surfaces 5^ if they are not ring-shaped. The 
equation of can now be written in the form 

(5.2.1) x = x(i,t],q), y = y[^^ri,q), z=^z{^,ri,q) 

If V IS the volume of a certain initial figure of equilibrium and its 
density is represented by « = h{a), we have the condition 

(5.2.2) + = = const. 

to be satisfied on S^. If the law of densities is changed to 

x^h{a.) Jr dxix) for the new set of surfaces of equal densities S', 
we shall have ® 

(5.2.3) a?! = a: af, y^ = y ^ 


— z ci^ 
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and, since the angular velocity has also been changed, 

dV^ + —= const. = C^{q) 


(5.2.4) 

a, 

Hence 


4- 

f 

t 

'v, Di 

(5.2.5) 



COiSf 0)^5^ 


Q - C = tC'{q) 


It was shown by Lichtenstein that the difference represented by 
the first two terms of the left-hand member can be expanded in a 
uniformly convergent series. 

(5.2,6) W^-W= + If<2) _f„ _ . 

The properties of this series are quite similar to those of the series 
(5.1.8)~(5.L9). It has been obviously assumed that new surfaces 
of equal density differ but little from the initial ones and that, of 
course, the difference -- co is small enough. Because of the 
variability of h, equation (5.1.10) must be replaced by a more 
general equation as follows: 


(5.2.7) Tf(i> 


dw r 
C + / 


A' 


dn 


Sq D 


C'da' - / 


V on D 


This expression holds for a continuous h. However, if there is a 
discontinuity of density along a surface we have to put 
[A'] == A(g* -|- 0) — A( 5 '* — 0) and 

(5.2.8) Pf W = Pf - / f [A'] 4- da' 

Js* D 

By (5.2.5)~(5.2.8) we obtain 


(5.2.9) 


dW r h' r dh' 1 

-j,dV' = fC ^ If(2)--pf(a) 
on JsD Jv on D 


_/H L dV'^-fsn~ ^ . 

Jki 2 


■ (co^+2/A) {ax-^ by) 
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On putting again 


+ Oi)^{ax 4- by) = /| 


the fundamental equation of the problem takes the form 

C h' r C' C dh' 1 

(5.2.10) =C'-5f 

Jvi £>1 2 / 

_i (ip(2) + p^(3)_p . . .) 


In the general discussion of this equation given by Lichtenstein, 
there are no points which in principle would differ from those given 
for the case of a homogeneous liquid. In order to obtain the equa¬ 
tion (6.1.12) which holds for this simple case, we have to put in 

(5.2.10) 

^=0. [4'] = 0 


The equation of Clairaut is obtained in this discussion by putting 
CO = 0(A ^ 0).^^ Lichtenstein's suggestion is to take four para¬ 
meters: d, X, C' [q], and Mg, where q is the value of q corresponding 
to the free surface of liquid, and to assume that these parameters 
take small values, i.e., that their moduli are smaller than a certain 
given small number. Then, for C, its expansion in a power series 

(5.2.11) C = 2 

can be used in order to compute new figures of equilibrium. 


5.3. Some Other Results of Lichtenstein*s Theory 

The general theory outlined in preceding sections could be 
obviously applied to a set of problems that present a certain in¬ 
terest from the viewpoint of Celestial Mechanics. Lichtenstein 

i.e., starting with a figure of absolute equilibrium. 
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discussed, therefore, the figure of a liquid mass rotating in the field 
created by some external masses, Laplace's problem of the shape of 
a liquid surrounding a rigid spherical core, as well as the liquid 
cylinder and the exact shape of the ocean. 

The basic point in which the treatment of all these problems 
differ is the shape of the given surface S in equation (5.1.14). In 
Liapounov's theory, the integral involved in this equation has been 
taken over the unit sphere and, therefore, the eigenfunctions could 
be expressed in terms of spherical, or Lam6 functions. In case of 
some other surface S, the main difficulty will probably consist in 
the kind of functions by means of which the eigenfunctions can be 
found in the form suitable for further discussion of the problem. 

A new extension of Liapounov's theory has been given by 
Lichtenstein [1] (Chap. V) in the application of the former's 
method to problems of this type: Suppose a given figure of a liquid 
mass satisfies conditions of equilibrium in first approximation; an 
exact figure of equilibrium in the vicinity of the first figure is to be 
found. 

It is obvious that the difference of potentials of two such bodies 
can be expanded in a series similar to (5.2.6), since for this part of 
analysis no condition of equilibrium is to be used. Then, a func¬ 
tion C can be defined to show departure in the shape of the second 
body from the first. Now the second figure satisfies the exact con¬ 
ditions of equilibrium while for the first only the first approxima¬ 
tion holds. Let co be the angular velocity of the figure of equilib¬ 
rium. Assume o) ((5) is a certain function of a parameter d, and con¬ 
sider a set of figures of equilibrium corresponding to it. Then, in 
the way shown in Section (5,1) an equation analogous to (5.1.12) 
can be obtained with three parameters C', X, and d. 

Obviously attempts can be made to solve problems of this kind 
for values of C', and d which are small enough. An essential dif¬ 
ference between this application of the method and the theory 
given in preceding sections must be kept in mind. The equation 
like (5.1.12) is obtained by subtracting two exact equations. As 
for the new case, one of these conditions holds only in the first ap¬ 
proximation as mentioned above. 
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The problems in which this approach could be of interest are 
according to Lichtenstein as follows: (a) a ring-shaped figure of 
equilibrium, (b) a ring with a central bodp, (c) a revolution of art 
almost spherical hodp around a mass point (central body); as a. 
primary figure an ellipsoid of Roche could be taken, (d) in particu¬ 
lar cases of the problem of n bodies when they are rotating as a. 
rigid system about a common mass center and all have almost 
spherical shapes, (e) double- or multiple-star systems with a central 
body, (f) a system of concentric rings. Investigations of several 
students of Lichtenstein concerning these problems and related, 
questions are listed in the Bibliography. 

5.4, Symmetry of Figures of Equilibrium 

All figures of equilibrium, the existence of which has been de¬ 
monstrated, namely, the eUipsoids of Maclaurin and Jacobi, the 
Poincard-Liapounov figures possess at least one plane of symmetry. 
This plane is perpendicular to the axis of rotation. The existence 
of it is a general property of figures of equilibrium and this fact has 
been proved in its general form for the first time by Lichtenstein 
[ 2 ]. 
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6J, Fundamental Equations 


The difficulties concerning the series expansions of the factor 
IjD in the expression for potential represented by the formulae 
(1.4.13) as well as the desideratum of Tisserand [1] (Vol. II, p. 317) 
were mentioned before. A method to overcome these difficulties 
has been suggested by Wavre and has been called the ''Uniform 
Method.'' In order to obtain the fundamental equations the 
solution of which, by this method, yields the figures of equilibrium 
we write first the equation (1.1.5) in the form 
( 6 . 1 . 1 ) 

The right-hand member of (1.1.5) obviously depends only on the 
density x it k = (p[p). For Q we shall have a more general ex¬ 
pression in the next chapter, but in the case of figures of equilibrium 
(2 is a linear function of and, therefore, = 2co^. We assume 
that 0 and its differential coefficients of the first two orders are 
continuous in the space V filled in with the fluid and that 0 is an 
analytical function in the part of space exterior to F. On taking 
into account the Poisson equation for the potential U, it follows 
from (6.1.1) 

(6.1.2) V20 = - 4:nfK + 2(0^ 

We assume now that 0 at the boundary S has a constant value 
05 . To find the expression for the value 0^ at some point M, we 
have to make use of Green's formula in two different cases. Let M 
be an exterior point and D = MM' its distance from a point ikf' 
of the body. Then, we have 


Hi 


V2 0 _ 0^2 




d<l> 

s \D in 


0 


i 


in D 




dS = 0 
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taking the inner normal. Since V2(l/D) = 0 and in the last term 
it is 0 = ^^5 == const., Gauss' integral is equal to zero and the 
equation takes the form 


(6.1.3) 


I 


D 


V2<5iF + 


r 1 m 

Js dn 



It holds for the outer part of space and along the boundary sur¬ 
face. If M is an interior point and a an arbitrary small sphere with 
the center at P, the surface integral in Green’s formula is taken 
over 5 and cr. It will yield two terms, Gauss’ integral for each 
of them being equal to 471 :. Thus, we obtain 

(0.4) 


and this equation holds in the interior and on the boundary. Sub¬ 
stituting (6.1.2) in the first term of (6.1.3) or (6.1.4) and taking 
into account, that 



xdV 

D 


we obtain 



- infU + 2co2 




By ( 6 . 1 . 1 ) we have and, therefore, (6.1.3) and 

(6.1.4) can be written in a S3mmietrical form 




C02 f 1 r 1 ^ 

On taking any level surface S*, the potential U can be rep¬ 
resented by the sum 


( 6 . 1 . 6 ) 




= f 

Jr. D 


xdV 


where is the volume bounded by S„. and L is the space between 
S and 5^. It is evident that if 5,^, is an external surface the density 
H in the layer L must be taken equal to zero. On applying the 
formula (6.1.5) to the first term in (6.1.6.) we obtain the equation 
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r X co^ r dv 1 r 1 d0 


<^S. 


The necessary and sufficient condition of equilibrium has been 
presented by Wavre in this form. As to the density distribution, 
it has been assumed that the density in a figure of equilibrium does 
not decrease toward the center. We have discussed this assump¬ 
tion in Section (1.3). The second equation (6.1.7) is a condition 
which must be satisfied on each level surface S* in the interior of 
a rotating liquid mass, and its solution will determine the stratifica¬ 
tion in the mass interior as well as the free surface that is the figure 
of equilibrium. In order to solve this equation, Wavre suggested 
its transformation as follows. Since d0jdn = g, we can write the 
second equation (6.1.7) in the form 


, ^ 1 r ^ , r ^ .r. co^ c dv ^ 


0^=0 


The first term can now be interpreted as the potential of a simple 
layer on the boundary of a "'cavity'' having the density equal 


tUz 



to g. In the interior of all the terms in (6.1.8) are analytic 
functions. Therefore, the left-hand member of (6.1.8) will be 
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determined by its expansion in a Taylor series, for example, for its 
value in the vicinity of the center. To obtain this expression, it is 
assumed that the volume of the cavity is divided in two parts. 
The third term in (6.1.8) can then be represented by a sum. At 
first we take the integral over the sphere of a radius equal to the 
distance of the point of intersection of the surface and the axis of 
rotation. It remains the integral over the space K enclosed be¬ 
tween this sphere and the surface S*. If i? is the distance of the 
point P in the ''cavity,'' and 6 the angle (R, co) we have 


(6.1.9) 


I 


v,1d 


2n 




Now, if the point M is close enough to C for a given (in 
general an arbitrary one), one can use the second Laplace expan¬ 
sion of IjD in a series (1.4.13). In fact, in 




Pnifi) 


for R < R' the latter condition holds for all three integrals in 
(6.1.8) the expression (6.1.9) being inserted in (6.1.8). In these 
equations R' is the distance of any element of integrals from C and 
fi = cos MXM, If we now take the conditions that coefficients 
of all R'^ in (6.1.8) vanish, we obtain 


1 

4jr 


Png ;,o , A Pn'X' 


dV 


( 6 . 1 . 10 ) 


2:7r 




for n=Q 

2 

w=2 


since 


Qp — 




0 w=l, 3, 4, 5,... 


p2(C0S 'd) = -f COS^ d' — 


This set of equations can be completed by an equation which is 
readily obtained from (6.1.2). If this latter equation is integrated 
over the whole volume of the liquid mass, the right-hand member 
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becomes 4,nfm + where m is the total mass. Since the 

left-hand member can be transformed into a surface integral by 
Green's formula, we obtain the equation of Poincar6.^^ 

= 4jr/m — 2a)^F 

If a similar transformation is applied to (6.1.2) and the volume 
V:^ taken, the equation 

(6.1.11) ^ f = 

47r Js* JL 27 c jk 

is readily obtained, as a supplementary condition. 

6.2. Outlines of Solutions 

Written in the form (6.1.10) and (6.1.11) the equations are not 
yet expressed in terms of those vaiiables which must be determined. 
In the choice of such variables Wavre's method differs from those 
discussed in preceding sections. The level surfaces are in this 
method represented in parametrical form by the equations 

(6.2.1) a; == i? sin ^ cos ^p, y = i? sin sin f, z = R cos 

R = R{q:^, yj) 

Thus, the assumption is made that the set of level surfaces is 
determined by the parameter q:^. It corresponds to a set of con¬ 
centric spheres with a radius varying from zero to a certain value 
given for the free surface of the liquid mass. 

Now, the conditions (6.1.10) and (6.1.11) can be transformed to 
introduce the variables R, f, but no essential progress can be 
made with a general form of the function i?(y*, yj). Therefore, 
Wavre assumed that 

(6.2.2) R = ^,,(1 4- f) 

From this equation, Poincar^ obtained an upper limit for the angular velocity. 
If the particles on the free surface of a figure of equilibrium are not separated, 
g has no negative values on S. Calling = mjV the mean density we obtain 
from the equation A lower value for the limit {co^ < Tcfn^) under 

certain conditions concerning the regularity of the surface 5 has been found by 
Crudeli [2] and Nikliborc [1]. Crudeli also proved [3], [4] another inequality 
restricting the highest admissible value of co. See also Prasad [1]. 
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where the function f = can be expanded in a series 

(6.2.3) C = + . . . 

These assumptions enable us to compare the results that can 
be expected by applying Wavre's method and those suggested by 
other investigators. On substituting (6.2.2) in (6.2.1) we obtain 
a parametrical representation of level surfaces different from 
(4.1.2) used by Liapounov or (5.1.1) of Lichtenstein. The para¬ 
meter corresponds to the case of spheroids. 

If in (4.1.2) instead of taking unequal axes, we put the radius of 
the boundary surface equal to it follows from (6.2.1) and (6.2.2) 

that The equations (4.1.2) represent the level surfaces 

in the final form adopted by Liapounov after he had tried some 
others in his first investigations. He realized that this form is very 
convenient for the exact proof of the existence of figures of equilib¬ 
rium differing but little from ellipsoids. For the law of density 
(4.1.1), the function t representing the departure of such level 
surfaces from ellipsoids could be expanded in a series (4.3.19). 
Thus, Liapounov's results will hold for any angular velocity which 
will determine the shape of an ellipsoid and for certain small but 
finite departure from homogeneity since <3 is not assumed to have 
infinitely small values. A similar approach is in Lichtenstein's 
method in which another parameter (A) is used. It is defined by 
the first equation (5.1.11), and this parameter can be used in case 
of a homogeneous liquid, when d of Liapounov and of Lichtenstein 
in the density law is identically equal to zero. Wavre, however, 
on making use of the series (6.2.3) restricts his considerations to 
figures which differ but little from a sphere. As mentioned in 
Section (2.3) only the terms of the order co^ were retained in 
Clairaut's theory and this assumption has been made in several 
further investigations. On considering this as a first approxima¬ 
tion, Wavre made use of the function f as represented by the equa¬ 
tion (6.2.3). On using the transformed equations (6.1.10) and 
(6.1.11) he has shown that by successive approximations the set 
of functions (similar to of Liapounov’s theory) can be 

expressed in terms of the parameter and of spherical angles, if, 
of course, the density k = H{q^) is a given function. 



METHOD OF WAVRE 


105 


Wavre’s function^^ C represents the deformation of a level surface 
with respect to a sphere of the radius and the effective calcula¬ 
tion of two first approximation yields figures of revolution. By 
computing the first approximation, Wavre obtained the classical 
results. We can easily see from (5.1.11) that, if the primary figure 
of equilibrium is a sphere, i.e., if we put co = 0, the parameter A 
becomes proportional to and co^ is the actual angular velocity 
of the figure of equilibrium to be determined. Thus, in this re¬ 
spect, the series (6.2.3) becomes a special case of more general ex¬ 
pansions used by Lichtenstein (see, for example (5.2.11)). In 
Wavre's method, the effective computations are to a certain degree 
simpler than those in Liapounov's or Lichtenstein's method. 
However, the latter methods are more general and exact. 

6,3. Other Investigations 

Certain modifications of the method of Wavre were suggested by 
Mineo [3] in order to simplify the successive approximations. The 
potential of a stratified heterogeneous body has been expanded 
again in a series of even powers of the angular velocity. Other 
expansions for potential, gravity force, and radius of the free sur¬ 
face of the Earth in terms of Legendre functions were also used by 
Mineo [1], [4] to determine the figure of a planet when the values 
of gravity at the free surface are given. The proof of impossibility 
of a homothetic stratification in a rotating fluid mass in equilib¬ 
rium given by Wavre was also simplified by Mineo [5]. The case 
of a strictly ellipsoidal stratilication in which the ellipticity varies 
from the center to the surface was discussed by Wavre in a more 
recent publication [2] and extended by Dive [3]. An exact ellip¬ 
soidal stratification is impossible in a fluid body, it x ^ cp [f) and 
the rotation is permanent. Other results concerning the figures of 
celestial bodies were obtained by Wavre not in direct connection 
with his method. We shall mention here only the extension of the 
theorem of Stokes to the case where the angular velocity is a func¬ 
tion of the distance s from the axis of rotation. In this form 


^ is e in Wavre’s notation. 
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Stokes' theorem is as follows: Newton’s potential in the space 
exterior to a celestial body depends only on the free surface (S) on 
the angular velocity co (s) and the total mass [M). Using this 
theorem, Crudeli [3] derived an integral equation from which the 
gravity at the surface can be determined. 



PART II 


Other Invariable or Varying Figures 




CHAPTER VII 


Zonal Rotation 

7.L General I^esults Concerning Zonal Rotation 

The theory of figures of equilibrium of a fluid mass yielded an 
explanation of many characteristics of celestial bodies. By the 
highest degree of precision in mathematical methods used to show 
the existence of a set of figures of equilibrium, a most solid base has 
been given to this theory. However, there are phenomena known 
in celestial bodies which require an extension of the theory for 
states other than an equilibrium. As mentioned before, the so- 
called zonal rotation is observed at the present time on the Sun, 
Jupiter, and Saturn in our solar system. It is well known that in 
certain nebulae the rotation follows a law different from that hold¬ 
ing for the Sun. Faye’s law (1880) represents very adequately the 

(7.1.1) CO = a + & sin^^?;^ {a == 14^.44, 6 = — 2^.31) 

zonal character of rotation of solar spots up to the heliocentric latitude 

= ± 35° discovered by Carrington. This empirical law applies 
also to other movements in the atmosphere of the Sun, showing 
that in general the angular velocity co is increasing at the surface 
with the distance of a particle from the axis of rotation. There is 
no such law given for a similar phenomenon observed on the Jupi¬ 
ter or Saturn. 

Many theoretical investigations dealt with the problem of zonal 
rotation.Initiated to explain Faye’s law of the zonal rotation 
of the Sun, they were extended to various other aspects of the 
problem as well as to the conditions of a general circulation of 
matter in a celestial body. The latter problem has been discussed 


Most complete lists of these investigations were given by Jardetzky [4] 
and Wasiutynski [1]. 
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ty Jeans [3] and in the work of Wasiutynski [1] in all detail. 
Such a detailed discussion is beyond the scope of this investigation. 
However, some further data will be given in Section (7.4). It is 
important to mention that there is as yet no general agreement 
about the physical phenomenon which could produce and maintain 
zonal rotation in a celestial body (see, for example, Jardetzky 
[13]). Nevertheless, it is quite natural to consider it as a possible 
second approximation in the case of a rotating fluid mass, when¬ 
ever the first approximation based on the assumption of a simple 
equilibrium cannot provide a complete explanation. 

In the first attempts to give a theory of zonal rotations, the 
equations of motion of a viscous fluid were used. Wilsing (in 1891, 
1896) suggested the h37pothesis that the core of the Sun rotates 
with a constant angular velocity and that in the outer layer, the 
distribution of co is given by the law co — funct. (s, z ), where s is 
the distance from the axis. The possibihty of the laws co = a (s^) 
or CO — co(z, s^} as well as the circular motion in a viscous fluid 
were investigated also by Sampson (1894), Harzer (1891-97) and 
Wilczynski (1896-97) (see Jardetzky [4]). 

The permanent and zonal rotations in a mass of an ideal fluid 
were treated by V6ronnet [l]-[3], Wavre [1], Dive [1], [2], 
Jardetzky [l]-[4] and others. 

The most important problem discussed in many papers was to 
find the stratification. In some of these investigations, the con¬ 
ditions were derived under which one could expect the existence of 
a given stratification. Thus, Harzer reached the conclusion that 
ellipsoids of revolution can be surfaces of equal density in a viscous 
fluid having a zonal rotation. (See also. Dive [2], Danoz [1].) 
The direct way to determine the figure of a fluid mass, when the 
law of rotation is given, i.e., the angular velocity is expressed in 
terms of coordinates, has been used by Wavre and Jardetzky. 

7.2. Fundamental Equation 

The function Q = {co^s^]j2 introduced in equation (6.1.1) is a 
particular form of the potential of acceleration as may be easily 
seen from equations (1.1.2) and (1.1.6). It is, therefore, evident 
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that the equation 

(7.2.1) /y + <3 = J—+ const. 


will hold for the general case where there is a potential of accelera¬ 
tions (a = — grad Q). As long as there are no other conditions 
imposed, one can postulate any distribution of accelerations in a 
fluid mass. For example, we can assume that a rotation occurs 
according to a law 


(7.2.2) w — co{s^) or (o{s^,z) 

In the first case on multiplying equations (1.1.3) by dx,dy,dz 
respectively and adding them, we obtain the term — co^ds^ 
(s2 = *2 -f 2/2). 

Thus 


(7.2.3) 


Q = |co(s2)is2 


We call the equation of the form (7.2.2) the law of rotation. For 
any form of this law which differs from that corresponding to the 
relative equilibrium, we can apply one of the methods discussed, in 
preceding chapters. The problem of figures of masses, the motion 
of which is subjected to a given law of rotation, can be reduced to a 
fundamental equation of the type considered before with a single 
change affecting the term representing the centrifugal force. Thus, 
for example, on restricting ourselves by the first law (7.2.2), the 
functional equation, the solutions of which will determine figures of 
a given fluid mass, is in its general form 


(7.2.4) fU + j(o{s^)ds^ = 0{>() 

There are no methods to resolve this equation without some 
new assumptions. We have to make use of some known represen¬ 
tations of the potential U and take a definite law of rotation. 
If we assume that the distribution of angular velocities in a mass 
does not differ essentially from a constant value of co, we can expect 
that a new figure of the fluid is a slightly deformed figure of equihb- 
rium. Whether such an assumption is correct or not, it shoidd be 
proved, of course, by a discussion of all conditions of the problem. 
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One of such conditions is immediately derived from the properties 
of the fluid. If the density k is constant or a function of the pres¬ 
sure alone [% = cp{^)) the equation (1.1,2) shows that the accelera¬ 
tion (7.2.1) has a potential. Now, if the motion assumed is re¬ 
duced to a pure rotation about a fixed axis ( 2 ') it follows from the 
equations (7.2.2) that 


dx 




dy 


= — 


dz 


= 0 


i.e., CO cannot depend on z. The law of rotation is of the form of the 
first equation (7.2.2).^® 

Wavre and Danoz applied the uniform method discussed in 
Chapter VI to zonal rotations determined by the law 

(7.2.5) co2 = a + hs^ 

where both constants, a and h, have small values. The case con¬ 
sidered by Jardetzky is more general because of the form of the law 
of rotation. He has assumed that 

(7.2.6) CO = coo + tui(s2, A) 

where: (a) coq can have any finite value; (b) only the parameter ^ 
takes small values; and (c) the function co^, is such that the func¬ 
tion Q can be represented by an absolutely and convergent series 

Ci)^ ^2 00 

(7.2.7) (?(s^A)=-^ + /2<?.(5V‘ 

The domain of convergence is given by the conditions 
0 ^ A ^ ;ii 0 ^ s ^ 

where is a certain limit value and a maximum distance from 
the axis of a particle at the free surface. 

Besides this, the function 


15 For the second law (7.2.2), no equation of the form (7.2.4) will hold 
(Bjerknes [1]). Attempts were made to determine a law of rotation co(x, y, z) 
such that it will correspond to a given characteristic of density distribution, for 
example, to an ellipsoidal stratification (method of Dive [2], for permanent 
rotations). 
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(7.2.9) 

&=2 

as well as its first derivatives with respect to coordinates should 
not be larger than where is a certain constant. 

Obviously, Faye's law (7.1.1) corresponds to the simplest case 
of (7.2.6), namely, to 

CO = (Uq + 

since odq is not necessarily small. This is also a more general form 
than the law (7.2.5). Under such assumptions, the figures of a 
rotating fluid mass which differ but little from an ellipsoid because 
of the zonal character of rotation could be determined by the 
method of Liapounov for a homogeneous as well as for slightly 
heterogeneous masses. Inserting the new value (7.2.7) of the 
acceleration potential in the equation (4.3.3) of Liapounov's theory, 
we obtain 

oo ^2 oo 

(7.2.10) 2 + - 57 = funct. (a) 

On assuming again that the level surfaces have the form deter¬ 
mined by a function f and represented by equations (4.1.2) or 
(4.3.1)-(4.3.2) we can make use of the expansion of Newton’s 
potential in a series (4.1,15). Then, the transformation of this 
series which yields the equation (4.3.16) applies to (7.2.10). The 
final result is, then, the equation 


(7.2.11) 


where 

, 1 

(7.2.12) 

4:7ta^A ^=1 


and W is given by (4.3.17). The function P {a) is determined by 
the formula (4.3.18). 

In these equations it cannot be seen whether an ellipsoid of 
revolution or that with three unequal axis can be used as a figure 
of reference. The assumption that the rotation is permanent 
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restricts, however, the choice of the figure of reference (Pizzetti 
[1]). For a fluid, we have the equation of continuity 

+ % div V = 0 
at 


Since the motion is permanent dvldt = 0 and dxjdt = 0. Thus, 
we have 


V • grad jc -(- % div v = 0. 

If the law of rotation is co = co(s^) and v is given by (1.1.1) 


(7.2.13)18 (iivv = 


doj dco lD(s^,co) 


dx 


dy 2 D{x,y) 


and we obtain the condition 


(7.2.14) v-grad« = 0 

The surfaces of equal density must be surfaces of revolution. 

The function C which yields the stratification in the fluid mass, 
must satisfy the equation (7.2.11) and the condition (4.1.6) 

(7.2.15) jCda= 

This condition is necessary in order to determine the arbitrary 
function of the parameter a introduced in the expression of P.. 


7.3. Figures That Differ But Little From Ellipsoids 

To solve the equation (7.2.11) we assume that C can be expanded 
in the absolutely and uniformly convergent series 

oo 

(7.3.1) f 

i+J>0 

where are functions of a and Since the right-hand member 
of (7.2.11) depends also on both parameters 6 and X, it has to be 
expanded in a series similar to (7.3.1), and we put, therefore 

(7.3.2) 


Thus, the equation (7.2.13) is necessarily satisfied by an incompressible- 
fluid, but the inverse statement is not true. 
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where are functions of a and ^(l^oo = 0)- On inserting (7.3.1) 
and ( 7 . 3 . 2 ) in ( 7 . 2 . 11 ), we obtain a set of equations 


(7.3.3) 




1 

J D (a, 1 ) 


= + /«(«) 


which must be satisfied by the coefficients Cii- The functions 
f^j-(a) are unknown but another set of equations derived from 
(7.2.15) will help us to find these functions. We have, namely, 

= - j (2C,^Sn^)^da - (I 


Hence 

(7.3.4) I = 0, jiarda = 0, JC 2 o^<^ = “ J^ 10 '^°'’ ■ • • 


Each integral 


(7.3.5) 



can be expressed in terms of those functions Csz for which k 1 
is less than i + j. Thus, N,, can be subsequently computed. 

If (7.3.3) is multiplied by da and integrated over the unit sphere, 

the equation 


(7.3.6) 


RN,, - 


1 

4:na^ 



da 

D{a^) 


Wda + (^) 


wiU be of the form (4.3.23). On taking into account the vaMe 
(4.3.27) of the second integral in the second term and eliminating 
from (7.3.3) and (7.3.6) we obtain the equation 


(7.3.7) 


RCi^ 


1 



1 

BiaTY) 


—^ tan“^ - 7 ”) da' 
a V e/ 

4:7C j 


1 

4:71 


R.N 0 


Now, the expression W ( 7 . 2 . 12 ) begins with the term of the 
order The coefficients wiU therefore depend only on those 
functions which have subscripts subjected to the same condition 
i J we had for N ,-j. 
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In the equation (7.3.7) we have the function Cij which is to be 
found for the whole range of values of the parameter a and the 
angle d'. In the second term the value is involved. It is taken 
for a = 1, i.e., for the free surface of the fluid. Thus, the natural 
way to determine this function is to put in (7.3.7) ^ = 1 again, to 
solve the integral equation 


(7.3.8) 


4:71 J Z)(l, 


-f const. 


and to substitute the function into (7.3.7). The bar denotes as 
before the value of an expression for a = 1. It is easy to see which 
terms in (7.3.7) after the substitution a == 1 yield the constant 
term in (7.3.8). are represented by the integrals (7.3.5) and 
are functions of a only. The same holds for terms — 1/471 J 
Now, to transform the third term of the left-hand member of 

(7.3.7) for a = 1, one has to put in (4.3.28) a = 1 and q = t = e; 
this term is then 


— Cij tan“^ - da' = const. Cij — const. = const. 

47% J Y Q J 

The function t (7.2.16) can be, therefore, calculated by successive 
approximations on resolving a double set of equations (7.3.7) and 

(7.3.8) which yield the functions 

(7.3.9) 

According to the first conditions (7.3.4) we have 

(7.3.10) iV,„ = iVo, = 0 

By (7.2.12) and (4.3.17) we can easily see, that 


(7.3.11) 


w. 


10 ■ 


4:7ia^A 


[0, - 0„(O)] 


^01 = 


4m^Zl 


CoiK 


The latter expression is obtained if we make use of the fact that 
in (7.2.12) are functions of s^. Because of == -j- = 



ZONAL ROTATION 


117 


^ 2(1 ^ _|_ 1 ) sin^ '& on a level surface Qj, depend on i.e., on 

the parameter 6 too. Thus, we have to put 

OO OO 00 

(7.3.12) 2 

^ ^ )c=l p=0 v=l 

As in the problem of the figures of equilibrium, the expressions 
and 0o(O) obtained by substituting = (p{a') in the cor¬ 
responding terms of Newton’s potential U (see, (7.2.15)) and 
(4.1.17). The effective calculation gives then (Jardetzky [4]) 

(7.3.13) W^o = W (?5(a')a'(tan-Ve-tan-V«)^^«' 
where s' is the positive root of the equation in t 

(7.3.14) sin^ ^ + 7 cos^ ^ = ( 7 ) 

In general, if the density law. i.e., the function (p{a) is given, it is 
possible to expand in a series of sphencal functions with 
coefficients which will depend on a. It was shown by Liapounov 
[9] (Part I, pp. 27-34 and 66-98) that under these conditions, the 
function Cio can be represented by the series of the form 

(7.3.16) ?io==2^2» 

where the spherical function Z^„ is given by the equation 


(7.3.16) 


Z^n 


y^n ^ 


4:7V 

4w 4-1' 


^ = COS & 


and PaJ/S) are Legendre polynomes. Jardetzky [3] has proved 
that the function Coi can be represented by the senes 


(7.3.17) Coi = 2 i^2«-Pan (/*) 

and the last two results show the difficulty and extremely 
long calculations needed to find the expression Ca in terms of 
coordinates. 
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If we restrict ourselves by the first approximation 

(7.3.18) ?i = fio<5 + CoiA 

we can make use of two known general expressions (7.3.15) and 
(7.3.17) on substituting certain particular laws of density and 
rotation. Liapounov calculated case where 99 (< 2 ) is a 

polynom in of degree k. In order to find a figure of a fluid mass 
in case of a zonal rotation Jardetzky made use of the density law 

(7.3.19) = 1 "t" dcpioCj = 1 -j- 7 ^ 2 ) 

where and are constants, and assumed that 
(7.2.6') CO = coo + A 52 

i.e., that the law of rotation is of the type of Faye's law. On taking 
the single term (7.3.18), oneassumes obviously that the parameters 
d and X are of the same order. Otherwise the terms in d and X 
should be combined in a different manner. After extensive trans¬ 
formations (Jardetzky [4] pp. 48-59) the equation of the free sur¬ 
face of the liquid mass which in general is given by the formula 
(4.3.2) could be written in the form 

(7.3.20) ^ + (7 

where , K^ are constants; Vq 1, Vq 1, and Vq are the semi¬ 
axes of the ellipsoid from which the figure (7.3.20) differs but little 
because of the heterogeneity of the liquid and of the zonal rotation 
represented by the laws (7.3.19) and (7.2.6'). The constants L^, 
can be expressed in terms of g, and coq , but we omit here 

these complicated formulae. 

From (7.3.20) we see that the free surface is represented in a 
first approximation by an algebraic surface of a simple shape. 

The density is a linear function of decreasing from its value 
Zo = 1 + at the center to = 1 + 5(i^, -f H^) at the free 
surface {Hi < 0, it d > 0). Maybe this law is too simple to be 
applied to celestial bodies, but it does not mean that is a linear 
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function of the distance from center. It is possible that it would 
be better to find the connection between a and 

The characteristics of the surface (7.3.20) depend on the 
numerical values of constants involved. The free surface will be 
either compressed or stretched along the equator of the ellipsoid of 
reference according to the conditions (^ = 0) 

(7.3.21) 2Lo<3 $ 0 

To see the deformation at the poles, one has to compute at 

= 0 and & n. For a negative f*, the figure (7.3.20) will be 
compressed at the poles since according to (4.1.2) we have 

(7.3.22) ± Ve(i + fn 

The shape of the meridian is easily seen from (7.3.20). 

7,4. Connection Between Zonal Rotation and Convective 

Currents 

The hypothesis of a zonal rotation of the kind considered in 
preceding sections represents also a certain simplification of the 
state of motion which is actually observed in many celestial bodies. 
It is highly probable that such rotation is only one part of a general 
circulation of matter in a celestial body due to rotation as well 
as to convective currents. In this section the broad field of 
internal movements in stars will be discussed only in brief, since 
as mentioned before, our aim is just to elucidate the part played by 
the hypothesis of zonal rotation in a very advanced period of the 
life of a celestial body, namely, during its solidification. In order 
to be able to develop a mathematical theory for this period, we 
can, of course, consider the zonal rotation as a given form of motion 
determined by certain initial conditions. Actually, a similar as¬ 
sumption has been made when the problem of a uniform rotation 
of a fluid mass has been formulated and solved. Nevertheless, for 
a better understanding of the evolution of a celestial body, the 
causes of a zonal rotation as well as the physical phenomena which 
are able to maintain this kind of motion should be determined. 

We have mentioned some earlier investigations which deal with 
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zonal rotation as a phenomenon characterized by Faye's law of the 
rotation of the sun. The general problem of the internal constitu¬ 
tion and internal movements in stars has not yet been solved, but 
there are several models which approximate the actual state with 
more or less precision. The simplest among these models is, of 
course, that of a gaseous star in an absolute equilibrium. The 
stratification is then represented by the density expressed in terms 
of the radius. The factors involved the pressure and the density, 
are connected by a certain differential equation. Then, some solu¬ 
tions of it have been found for different cases of an adiabatic equi¬ 
librium (Emden [1]). The assumption usually made is that in this 
case, i.e., in the case of a "'convective equilibrium” the heat transfer 
by convective currents from the interior to the surface of the star is 
included in the general state of the body. But other factors must 
be taken into account, and Jeans [2] pointed out the role played 
by the" "viscosity of radiation. ” According to one of his conclusions 
(Jeans [3] p. 91), there will be no convection currents in a model 
corresponding to a configuration of equilibrium of the type dis¬ 
cussed except maybe near the surface. 

We do not discuss here the sources of energy in celestial bodies. 
The discovery of nuclear processes overshadowed all previous con¬ 
siderations about such sources. Models concerning rotating stars 
are obviously more important than a static model. Jeans, in the 
paper just mentioned, [2], has shown that the inner parts of a star 
must rotate more rapidly than the outer layer. One of his con¬ 
clusions (Jeans, [3] p. 265) deserves our special attention. In his 
opinion, namely, throughout the whole life of any star whatever, 
""the equalization of angular velocity produced by viscosity is 
negligible in the central regions of the star.” In a simplified model 
he shows that, if the angular velocity of rotation at an initial 
moment depends on the radius only, i.e., if the rotation occurs by 
spherical shells, it will not be destroyed by the regular viscosity or 
by the radiation. In the second approximation, however, Jeans 
considered Faye's law of rotation. The fact that the ellipticity of 
spheroidal layers increases towards the center should yield, namely, 
the explanation of the equatorial acceleration. The factor pro- 
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ducing this phenomenon is the drag exerted by inner layers. Thus, 
according to Jeans, a zonal rotation is not necessarily connected 
with convective currents. However, the model of spherical shells 
each rotating with its own angular velocity seems to be too simple 
an approximation. Among the general assumptions concerning 
the Sun, we will find the conclusions based on numerous observa- 
tions as follows (V. Bjerknes [1]): (a) the Sun is near a state of 
stable internal equilibrium, (b) conditions in general seem to be 
favorable for the formation of stratified circulation, (c) the origin 
of all that happens on the Sun is in the solar radiation, (d) heat 
lost by the photosphere is restored from the deeper layers by inter¬ 
nal radiation, convection, and conduction, (e) internal radiation is 
believed to be the dominant factor especially in the deep layers of 
very high temperature, (f) the importance of convection increases 
near the surface, and (g) the rotation of the Sun tends to produce 
general zonal symmetry. 

Obviously, all these facts are beyond the limits of the theories ex¬ 
posed in preceding sections. We could introduce in the equations 
of motion the law of rotation of the form co = o{s^) and apply one 
of the methods to solve the resulting integral equations. But tak¬ 
ing into account new physical factors just mentioned, the equations 
of motion have to be modified as well as supplemented by a set of 
new conditions. Under new assumptions it is, for example, shown 
(Krogdahl [1]) that if a celestial body is composed of a viscous 
fluid its rotation would give rise to meridional currents, which 
would cause the zonal character of rotation. In Bjerknes' theory 
the existence of such currents was just postulated. 

We know several special conditions of equilibrium derived for 
different kinds of physical phenomena. For a stable equilibrium 
of an incompressible ideal fluid when none but mechanical 
phenomena are considered the condition 

dr 

must hold as proved by Rayleigh in 1917,^’ It concerns the cir- 


See also Randers [1]. 
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culation along any parallel. An extension of this condition is 
found in the Solberg-Hoiland criterium of stability for compressible 
fluids (Wasiutynski [1]). The formation of large-scale currents 
could result in case of gravitational instability. However, if the 
conditions of gravitational stability are satisfied, the convective 
current may arise because of the fact that the gravitationally 
stable layers are not necessarily also in a thermal equilibrium 
(Vogt, Eddington [1], Jeffreys [1]). In such a case the resulting 
pressure gradient on a level surface would tend to make the matter 
flow from equator to poles, or vice versa. The equilibrium could 
be maintained if the temperature distribution does not give rise to 
this component of the pressure gradient, but when the radiation is 
also taken into account, this will lead to a new form of conditions 
of equilibrium. According to v, Zeipehs theorem [1] the genera¬ 
tion of energy 4:7zE (per second and gram) within a rotating gaseous 
mass must follow the law 



where 5 is a constant. It was pointed out by Eddington [1] that 
the generation of energy from nuclear sources ''can scarcely follow 
a law of this kind.'' Among possible internal movements which 
would result, since the equilibrium cannot be realized, the primary 
currents in planes through the meridians seem to be the most likely 
phenomenon according to Eddington. These currents must be, 
then, deflected (East or West) by the rotation of the body, and 
different periods of rotation will result in different latitudes and at 
different dephts. In this general characteristic both cases are 
included—an angular velocity increasing in a layer with the in¬ 
creasing distance from the axis of rotation as well as a decreasing 
angular velocity. 

The regular viscosity and the radiation (or the radiative viscos¬ 
ity) are obviously factors acting in opposite directions. The 
former tends to even out the differences in angular velocities while 
the latter can establish currents and, therefore, such differences. 
From the equations of motion of a viscous fluid Eddington [2] 
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derived that the currents which will be maintained by the radia¬ 
tion in the interior of the mass will correspond to the circulation 
near the surface from the poles to the equator. 

The actual complexity of the problem of convection in celestial 
bodies is, of course, greater than it was assumed in the cases in 
which an approximate solution was attempted. Even some con¬ 
clusions drawn from earlier observations seem to be taken in too 
simple a form. For example, Tuomiiien [1] found that the sunspots 
at latitudes lower than 16° have a drift towards the equator while 
the spots at higher latitudes move towards the poles. Thus, there 
can be two opposite directions in the surface layer for convection 
currents and there are no reasons to exclude the existence of a still 
more complicated distribution of velocities. 

At the present stage of the theory of internal motion in celestial 
bodies we must point out that as long as pure dynamics is in¬ 
volved, both cases, the angular velocity increasing or decreasing 
with the increasing distance from the axis of rotation, are equally 
possible. In many investigations in which the problem of con¬ 
vective currents is considered, under special conditions a decreasing 
angular velocity resulted. Nevertheless, sometimes an angular 
velocity increasing in a layer near the free surface from poles to 
equator was also derived. As to the actual state in celestial bodies, 
the outer layers of planetary nebulae are rotating usually less 
rapidly than the core. Zonal rotation is characterized in the case 
of the Sun, Jupiter, and Saturn by an opposite law. The well- 
known experiments of Belopolsky and Riaboushinsky showed that, 
if in a glass sphere filled with a liquid differences in angular velocity 
are produced, the currents along the axis of rotation towards the 
center and moving along an equatorial radius from it are connected 
with an angular velocity decreasing as the distance from the axis of 
rotation is increasing. These currents and those considered by 
Eddington have opposite directions. In the experiments just 
mentioned, there is a new factor involved, namely, the friction at 
the surface of the fluid which obviously does not play any part in 
stars. It seems, however, that the effect of it does not change 
essentially the general character of internal movements, and the 
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experiments of Belopolsky [1] and Riaboushinsky [1] illustrate 
well the zonal rotation. 

On the other hand these experiments correspond better to con¬ 
ditions for convective currents which one can expect to find in the 
fluid core of a planet or of a star already covered by a solid crust. 
Then, the friction between a very viscous liquid interior and the 
crust can become an important factor. 

The studies of such an advanced stage as a solidifying celestial 
body will certainly meet new theoretical difficulties, since, for 
example, the form of the equation of state for plastic and solid 
bodies is much more difficult to determine than that for gases or 
liquids. 

A very important general conclusion has been reached by 
Joukowski [1] in the problem of motion of a rigid body having a 
cavity filled with a viscous liquid. In this problem the friction at 
the wails of the cavity was also taken into account. According to 
Joukowski if there are convection currents in the viscous liquid 
in the cavity, they will be directed from poles towards the equator 
along the surface when the mean square of angular velocity is a 
function increasing with the distance from the axis of rotation. 
The opposite distribution of angular velocities will take place if 
the direction of convective currents is reversed.^® 

7.5. Zonal I^otation in a Core 

In order to study the conditions of the problem of internal 
movements in a viscous core of a planet, we shall write in a vector 
notation the equations used by Joukowski [1]. To simplify the 
problem an incompressible liquid is considered, and, therefore, we 
have 

(7.5.1) divv=0 

where v is the velocity. Let S2 be the angular velocity of a par¬ 
ticle. Since curl v = 2£2 and curl curl v — grad div v — V® v, 
we obtain 


This statement is sometimes quoted in the literature in a not quite exact 
form. The paper of Joukowski [1] has been published in Russian. 
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(7.5.2) V == — 2 curl £2 

Now, the equation of motion of a viscous liquid is 


(7.5.3) 


1 

a = p-grad p -f- vW v 

K 


We assume: (a) that P== grad U, i,e., is derived from Newton's 
potential, (b) the kinematical coefficient of viscosity is constant, 
and (c) the liquid is homogeneous {k = const.). The differential 
equations of motion can be written in terms of components of the 
angular velocity, and the Helmholtz equations represent such a 
form for perfect fluid. It is known that to obtain these equations 
one has to apply the operation curl to (7,6.3). Then, if 

1 

curl P = 0, curl — grad ^ 0 

H 

we obtain 

(7.6.4) ^ == (J2 • V) V + — curl v 

dt 2 


To transform this equation we make use of the known formulae as 
follows 


dSl 

dt 


an 

+ (vV)J3i 


(7.5.5) [^xS2] = Vx[vxn] 

= v(V-a) + (n-V)v-n(V- v)-(vV)n=(a-V)v-(vV)£2 


because of 


V • V = div V = 0, V • SI = I- div curl v — 0 
Now, by (7.5.4), (7.5.2), (7.5.5) and putting 

(7.5.6) curl SI = 2Slj^ curlSi^ == 2Slf^ 
we obtain 

, dSl 

(7.6.7) = curl [vXSJ] — 4i;S22 
ot 

Thus, the forces are eliminated from the equations, and the 
initial and boundary conditions for the liquid part must now be 
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given. Jonkowski made use of a simple form of boundary con¬ 
ditions. They express the fact that friction acting at a point of the 
internal boundary surface of the crust in the direction of the rela¬ 
tive velocity (?y^) of the liquid is balanced by the internal friction 
of the core. If r is the coefficient of external friction at the 
boundary of the liquid and ^ the coefficient of viscosity, we have 
the first boundary condition in the form 

(7.5.8) xv^ = 

on 

The second one must express the continuity of the normal 
velocity component at the boundary of the cavity 

(7-5-9) = [^'Jorust 

The general scalar equations which can be obtained from the 
vector equation (7.5.7)^® were used by Bondi and L 5 d;tleton [1] in 
an investigation concerning the internal movements in a liquid 
core. One of their conclusions is of great importance at least in 
case of the Earth, and we shall present an abbreviated form of the 
mathematical analysis of the problem. 

Different generalized orthogonal coordinates q^, can be used 
and the corresponding scalar equations derived from (7.6.7). On 
representing the line element by the formula 

ds^ = + S^dq\ -f- S^d(^ 

and denoting by y,, v^, and cdi, co^, coj the components of velocity 
V and angular velocity respectively in directions corresponding 
to ^ 1 , q^, q^, we write the first equation of motion, which follows 
from (7.5.7) in the foiun 


dcoi 

w 

(7.6.10) 


•S'a 


I ^Ss{vj_co2—v^o)^)~~Sz{v3COj^—V]^a)s) 
Loq^ dq^ 

a r 9 a 

-(■5'2"2)-^ (•S'l<«l) 


SgSg 13^2 ^ 1^2 


d S, 


'^3 




IB Jardetzky [10]. 
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The second and third equation are readily obtained by changing 
the subscripts. For further investigations the cylindrical coor¬ 
dinates r, y are used, and we have = 

dz^ + dr^ + Hence == Sg = 1, Sg == r. Bondi and 

Lyttleton introduced an operator to simplify the form of equations 
and for the cylindrical coordinates this operator is defined as^® 


9^ 9^ 19 

(,.5.11) 


02 1 0 
9^^ ^ 9 r 7 dr 


It can be used in case of axial symmetry. Assuming that all the 
factors involved in (7.5.10) do not depend on the angle y (longi¬ 
tude) these equations can be reduced to a much simpler form. 
From the expression of div v in terms of generalized coordinates 


div V = 




[vxS.S,) + {v,S,S,) + ~ {v,S^S,) 


it may be easily seen that the equation (7.5.1) for cylindrical co¬ 
ordinates and axial symmetry will be satisfied if we put 


(7.5.12) 


dip 


dip 

dz 


The well-known formulae for the components of curl v = 2£2 
yield the expressions 


(7.5.13) 


CjO-i 


CDo 


(Ug 


2S2S3 

Idq^ 

1 

■ 0 

2S,S^' 

[0^3 

1 1 

■ 0 

“ 251 S2 

Ls?! 


^ {S,v,)~-^{S,v,) 

( 53 ^) 3 ) 


On putting 

(7.5.14) S 3 2^3 rv^ ^ X 


the expressions (7,5.13) take the form (for cylindrical coordinates) 


20 The symbol is used to avoid the confusion with the distance D in preceding 
chapters. 
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(7.5.15) = 


1 


h 

dr 


COg = 


J_ ^ ' 3^2 

2r dz ~ T 


By (7.5.12), (7.5.15), and (7.5.11) we have 


dv-^" 
dr _ 


(7.5.16) 0)3 = ~ 


2r 


d^f 9 1 df' 

. dz^ dr r dr. 



W 


Thus, six unknown variables Vf, co^ are expressed in terms of two 
functions y)(z,r,t) and x(z,r,t). To determine these functions, 
we can make use of the first and third equations ( 7 . 5 . 7 ) on sub¬ 
stituting (7.5.12), (7.5.14), (7.5.16), and (7.6.16). These equa¬ 
tions were linearized by Bondi and Lyttleton. They obtained first 
the equations 


(7.5.17) 


/ r Ldz dr dr dz 


= 0 


2% dy 
r^ dz 


and assumed that the motion in the liquid core differs but slightly 
from the rotation of a rigid body around an axis fixed in space. 
Under this condition let 

C^-b.lS) ^ = (n + kt)r^ -f- z) 

where n and k are two constants and (p a function having small 
values of the order oi y. Bondi and Lyttleton considered the 
action of tidal friction in the case of the Earth and assumed also 
that n has a small value while A < 0 and \k\ < n^. We make 
similar assumptions which wiU correspond to very small changes 
m G)i, 0 ) 2 , 0 ) 3 . Then, on neglecting in (7.5.7) the products of other 
small factors and even the term having a factor kt, the equations 
(7.5.7) become 

-\-2n~= kr^ 

(7.5.19) 

To simplify more these hnearized equations a very slow depar- 
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ture from the rotation of a rigid body is considered, the function cp 
is substituted in the first term of (7.5.19), and the time variation 
of 'ip is neglected. Then, we have 

- d'w 

vD^ <p + 2n — = kr^ 

uZ 

7.5.20 . 

A. d(p 

vD^ f — 2n — = 0 
oz 

Obviously the assumptions which made possible the derivation 
of (7.5.20) can hardly be adopted for gaseous bodies like stars. 
They seem to be plausible nevertheless when we consider very slow 
movements, for example, in the interior of the Earth during dif¬ 
ferent geological periods. 

The solutions of (7.6.20) obtained by Bondi and Lyttleton re¬ 
present those internal movements for which the speed of circula¬ 
tion about the axis of rotation is of a higher order of magnitude than 
the speed of circulation in meridian planes, i.e., the former is 
much larger than the latter. This is a conclusion, very important, 
for the theory of the Earth's figure, since according to the hypoth¬ 
esis of zonal rotation, the principal features of the Earth's crust 
were produced by internal movements of similar kind. 

7,6. Internal Movements and Solidification of a Celestial 

Body 

There are but few planets the surface features of which can be 
observed, and the structure of the surface layer of only the Earth 
and the Moon is known in more detail. The variety of causes which 
can affect the formation of a crust on a cooling celestial body can 
result, of course, in many different structures. The varying field 
of forces of attraction and resulting tidal action, the falling meteors 
or the fission of a body, the nonuniform distribution of matter 
may play a more or less important part in deformation of a solidi¬ 
fying outer layer of a body. Special attention, however, must be 
paid to internal movements as to one of the most important factors 
acting during the formation of a solid crust on a celestial body. 
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To take into account such internal movements one has to assume 
that the interior of a planet displayed the properties of a fluid. 
Whether the usual assumption that the planets were formed from 
a very hot fluid which was cooling during the whole evolution or 
the suggestion that the Earth, for example, has formed a liquid 
core during geologic time is correct is not essential for the hypoth- 
esis of internal movements. The latter possibility, namely, that 
the Earth’s core is formed from a state of a nearly uniform mixture 
of the iron and silicate phases and that in the earlier period of its 
existence our planet increased its temperature by contraction of 
the gas and dust mass is suggested by Urey [1]. We do not specify 
here the chemical and physical processes in the Earth’s interior, 
but at the present time the theory as well as observations reduced 
the probability of existence of a celestial body without currents of 
matter in its interior practically to zero. The condition of a rigid- 
body rotation of any celestial body has been very useful in the 
theory as long as an explanation of figures on a large scale was 
gained. The statement that a factor like viscosity could even out 
all the differences in velocities before a planet solidified was always 
an assumption suggested only by the lack of knowledge about the 
sources which could maintain these differences. As mentioned 
above, the discovery of sources of energy such as the carbon cycle 
of nuclear reaction made this assumption still less probable than it 
was earlier. The single criterion in this case must be taken from 
the observation. If there are no traces of deformations in the 
solidified portion of a celestial body the internal motion in liquid 
parts obviously did not play an important part. If the structure 
of the crust shows that there were certain displacements we cannot 
yet assert with certainty that they are due exclusively to the action 
of internal movements. Several theories involving the existence 
of currents in the Earth’s interior were suggested. 0. Fisher 
seems to have been the first to point out in 1889 that convection 
currents could produce certain tectonic phenomena. This idea 
was applied in different forms by Ampferer, Andr^e, Schwinner, 
Holmes, Griggs, a.o. Solutions of particular problems were obtained 
dealing with thermal disturbances like the difference in tempera- 
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ture distribution under the continents and suboceanic parts of the 
crust or with cellular convection of B6nard type. According to 
Urey [1], the assumption that the cooling at the poles or beneath 
the oceans is the single cause that produces the currents in the 
mantle leading to mountain chains folding does not appear possible. 
He admits, however, (Urey [2]) the possibility of a planetary 
convection which could occur at the time of a high viscosity in the 
Earth. Jardetzky [13] pointed out that none of these theories of 
subcrustal currents were able to explain the characteristic facts 
known about the formation of the Earth's crust. In general, how¬ 
ever, according to Vening-Meinesz's opinion [1] the best hypoth¬ 
esis explaining the geological processes has to make use of cur¬ 
rents in plastic subcrustal layers. These currents could be produced 
by the temperature-gradient maintained by other sources of 
energy or due to internal rearrangement of bodies of different 
densities. Urey [1] also pointed out the role of such rearrange¬ 
ments. On discussing the problem of convection in a sphere^ as. 
well as in the mantle and representing the distribution of current 
system by spherical harmonics of first five orders, Vening-Meinesz 
[2] makes these currents responsible for the formation of the 
protocontinent and for its later distortion, i.e., reaches the same 
conclusion which has been made by the author for the first time 
in 1929. (See also, Gutenberg [1]). The problem of thermal in¬ 
stability in a planet heated within is considered again by Chan¬ 
drasekhar [4], In order to understand the generation of convec¬ 
tion currents, we should know in a more exact way in which fac¬ 
tors produce such differences in temperature in any celestial body 
as well as in the Earth. The simplest assumptions we used before 
are as follows: in an advanced stage of a celestial body at which the 
formation of the solid crust began, one can expect that the differ¬ 
entiation of substances is more or less completed, the stratification 
corresponding to the distribution of forces is settled, and the tem¬ 
perature field is determined by the condition of a thermal equilib¬ 
rium. By rejecting the latter assumption, the way was open to 
locate convection currents at any part of a celestial body where 
they were desired. This explains different aspects of the hypoth- 
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esis of convection currents and also possibly that fact that none 
of them are generally accepted. However, a great deal of truth 
seems to be in this hypothesis. It should be given, therefore, such 
form that the conclusions drawn from it can be verified in a more 
precise way than it has usually been done. 

The fact that internal movements could occur in every celestial 
body until it is solidified is self-evident. These movements could 
possibly have an arbitrary distribution at some earlier period. 
Taking into account the effect of other factors, one can expect, 
however, that a certain regular distribution has been imposed on 
these currents later. It will be very difficult to give a mathemati¬ 
cal theory which will prove that a certain type of internal move¬ 
ments must occur, In many cases we should be satisfied by a proof 
that this type of movement is possible. The investigation of 
Bondi and Lyttleton provided such a proof for the zonal rotation 
in a planet. As it was shown in the preceding section, in a system 
of convection currents having an axial symmetry, the circulation 
about the axis can be the dominant part while the displacements 
in meridians can be neglected. It must be pointed out that this 
character of internal movements is not necessarily imposed on the 
displacements in every celestial body, Whether the relative move¬ 
ments in the Earth's interior had such a zonal character or not, 
can be inferred only by comparison of the logical consequences of 
this assumption to the actual distribution of the principal features 
in the crust. 

7 . 7 . Formation of Principal Features of the Barthes Crust 

It has been shown by the author (Jardetzky [1], [4], [5], a.o.) 
that the h3q)othesis of zonal rotation yields an explanation of 
many of the principal characteristics of the Earth's crust and does 
not contradict others for which complementary causes must be 
suggested. The complexity of deformations during geological 
history as well as the present shape of the crust require an agree¬ 
ment as to which facts are to be considered as the principal ones. 
By neglecting such a definition, much confusion has been produced 
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in the discussion of some geophysical theories, and, therefore, and 
attempt has been made to present a list of most important facts 
which should be taken into account in the first place (Jardetzky 
[14]). (a) The outer layer of the Earth differs but little from a 
spherical shell; the density is increasing with depth, (b) The com¬ 
plexity of the structure is characteristic to the uppermost layer 
(33 km. thick). A layer now termed the crust is petrologicaly 
separated from the rest of the solid part of the Earth which is 
called the mantle. Here we use the term crust in the sense of a 
solid shell, (c) The core should be characterized as liquid even at 
the present time according to seismic information. Its radius is 
3400 km. 

The hypothesis of zonal rotation is compatible with this group 
of data. From (a) we can draw the conclusion that at the time 
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when the mantle became thick enough to resist any deformation 
on a large scale, the rotation of the Earth was almost a rigid-body 
rotation with a small angular velocity. Slow internal movements 
did not essentially change the stratification which differs but little 
from that in a heterogeneous figure of equilibrium (corresponding 
to the law of densities which is represented in Figure 7, Bullen [2]). 

Obviously, the most essential effect of internal movements could 
be provided during the formation of the crust and eventually of the 
first thin sohd layer in the mantle. The data concerning the core 
at present support in general the idea that parts of the crust and 
mantle were fluid during some geological periods. Now, the most 
important general facts concerning the crust are (d) the continents 
are mostly granitic, while the bottom of oceans is basaltic; (e) the 
actual distribution of land and sea for which probably one of the 
most striking characteristics is the existence of geographical 
homologies; (f) the changes in this distribution in different geo¬ 
logical periods; (g) the formation and distribution of mountain 
chains; and (h) the phases in the formation of mountain chains. 

The hypothesis of a zonal rotation yields a very simple explana¬ 
tion of the horizontal differentiation mentioned in (d). It shows 
also in which manner the Old World could be separated from the 
Western continents and dragged by the underflowing magma to 
the East. Australia, the East Indies, and many islands in Pacific 
could be broken away from Asia and Africa and also displaced to 
the East by the equatorial belt of fluid substratum accelerated in 
this direction. Thus, the hypothesis of a zonal rotation of the 
Earth's interior seems to be in very good agreement with principal 
characteristics of the actual distribution of land and sea. As to 
the changes in continental blocks and in ocean basins, a special 
investigation is necessary to estimate the effect of the magma 
currents during different geological periods. There is no doubt that 
some changes were due to other causes than the drag along the 
crust-mantle interface. The deformation of the whole Earth’s 
body because of its displacement with respect to the axis of rotation 
could also be an important factor. Tides in a still fluid interior of 
the Earth when the crust was relatively thin as weU as departures 
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from a simple law of zonal rotation (discussed by the author) 
probably affected the consolidation of the crust at many places. 
In the formation and distribution of mountain chains, however, 
the stresses produced by the zonal rotation of the magma again 
played the most important part. The Tertiary mountain belt, for 
example, has been easily explained by this hypothesis as well as 
some other mountain ranges (Jardetzky [5], [7], [8]). In order to 
understand the formation of remaining ranges of mountains (like 
the Appalachians) more data are needed to find the field of forces 
resulting in those folds. The last point mentioned above, namely, 
the cyclic character of the mountain building, i.e., the phases in 
the deformations can be explained by the interference of two fac¬ 
tors: the increasing viscosity of the underlying magma and the 
advancing thickening of the crust. The second factor shows that 
larger stresses are necessary in later geological periods to produce 
the folds, while the first yields increasing horizontal stresses which 
can reach a value required for such a folding. However, it must be 
pointed out that without the hypothesis of the zonal rotation it 
cannot be seen so easily why the field of these forces existed 
exactly at the places which are determined by the actual distri¬ 
bution of mountain ranges. 




CHAPTER VIII 


Varying Figures 

8,1, Small Oscillations 

A set of figures of an isolated fluid mass the existence of which 
could be proved by different methods includes stable figures of 
equilibrium or figures of a fluid in a permanent zonal rotation as 
well as unstable configurations. The oscillations about a stable 
figure and progressive deformations represent two types of prob¬ 
lems which help to elucidate the behavior of a celestial body during 
a more or less long time interval. More attention has been paid in 
celestial mechanics to the former problems. The periodic changes 
investigated were either the free oscillations of a fluid mass or 
small changes due to the action of external periodic forces. 

Poincar6 [1] studied the free oscillations of a liquid ellipsoid, but 
the question whether it is possible that all particles of a liquid mass 
have the same period of oscillations needed further discussion. 
It was shown by Appell in 1920 and Cartan [1] that the conclusions 
of Poincar6 are correct despite the doubts voiced by Globa- 
Michailenko [1]. 

Hill found in 1884 that the integral of kinetic energy in case 
where the particles of a fluid mass have the velocity 

(8.1.1) v = 'V(p + mVy) 
can be written in the form 

( 8 . 1 . 2 ) 

This equation can be used now instead of (1.1.5). 

On assuming that the level surfaces are determined as before by 
a parameter a, (8.1.2) takes the form 
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(8.1.3) /[7=-^+J(V 95 ) 2 -—(Vy)^ + funct. {a) 

Thus, for example, the series expansion given by Liapounov 
(Chap. V) can be applied as long as these surfaces differ but little 
from ellipsoids and the equation can be solved for the function f 
which will determine the deformation. This would be possible, of 
course, if three functions 9 ?, m, yj were known. They must, however 
satisfy certain general conditions given by Clebsch [ 1 ], and the 
problem becomes much more difficult than problems concerning 
figures of equilibrium. Jardetzky [4] has calculated the first 
approximation for periodic changes of a meridian of an oscillating 
homogeneous ellipsoid using the equation (8.1.3), a given law of 
oscillation, and the method of Liapounov. 

The periodic variations of a liquid ellipsoid such that at every 
instant the shape of the liquid is exactly ellipsoidal were deter¬ 
mined by Dedekind [ 1 ]. (See also Lamb [ 1 ].) 

It is well known that methods different from those discussed in 
preceding chapters were applied to the problem of small oscilla¬ 
tions of a liquid mass under the action of periodic forces. The 
theory of tides represents a broad field, the results of which are not 
used in this investigation. 

8,2. Progressive Deformations 

The problem of progressive variations of a celestial body has 
first been met in a particular form. It has been shown, namely, 
that there are continuous sets of figures of equilibrium, and a 
natural conclusion was drawn from this fact. It has been assumed 
that during its evolution, a celestial body will be subjected to such 
changes that its free surface will coincide subsequently with figures 
belonging to one of these sets. It is evident that this is too idealiz¬ 
ed a case. In that aspect, in which this problem was posed in 
order to draw some conclusions about the evolution of stars, no 
essential progress has been made to our knowledge about the pro¬ 
gressive deformation (see, for example, Jeans, [ 3 ]). There is no 
doubt that the case of an isolated homogeneous mass is of second- 
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ary importance. In most cases there are factors like action of 
other members in the solar system, which will change such an 
undisturbed evolution. 

Even if we cannot expect to have more or less important applica¬ 
tions of the theory of a varying figure in case of every celestial 
body, this theory should be more developed because of the de¬ 
formations of the Earth. 

We can see two classes of deformations in case of the Earth. 
To the first belong all relative displacements of the parts of the 
outer layer. There seems to be no disagreement now about the 
assumption that the folding of greater mountain ranges involves 
horizontal displacements to 100 km. (Jeffreys [1]). There arises 
only the question about the horizontal displacements of the order 
of several thousands kilometers, as it is required in the theory of 
continental displacements. All such variations in the mass distri¬ 
bution would affect the shape of a geoid, and this figure in Car¬ 
boniferous was different from that determined at the present time. 

The horizontal deformations in the crust as well as some epiro- 
genic displacements could result from internal movements, for 
example, from the zonal rotation in the Earth as it was pointed 
out in Chapter VII. 

A deformation of another type could be due to displacement of 
the axis of rotation in the body. The migration of the poles of the 
Earth has been discussed in geophysics and related sciences and 
according to one interpretation of geological data, the poles moved 
progressively over the surface of the Earth during all geological 
periods. One of the most striking phenomena in the past, the trans¬ 
gressions and regressions, changing the distribution of land and sea 
in many different ways could be at least partly explained by the 
displacement of the axis of rotation in the Earth. But, in general, 
the rotation of a deformable body is a problem which is still far 
from a complete solution. The basic equations of motion of a body 
which can change its shape during the motion were given by Liou- 
ville (see, for example, Tisserand [1]). Darwin [2] attempted to 
determine the influence of geological changes on the Earth's axis 
of rotation by assuming that the body is slightly deformable and 



140 


FIGURES OF CELESTIAL BODIES 


the changes are slow. Different conditions concerning this prob¬ 
lem were investigated also by Picart [1]. Darwin and Schiaparelli 
introduced the hypothesis of continuous adaption of the body to a 
new figure of equilibrium which results from the displacement of 
the axis of rotation with respect to the axes of inertia, since from 
the equations of Liouville one cannot determine completely the 
motion of a deformable body. In general, as is well known, in such 
a body the concept of the axis of rotation must be introduced in a 
more precise manner and this can be done in different ways (see, 
for example, Bilimovitch [1]). The most useful seems to be the 
axis of rotation of an equivalent rigid body which is defined by 
the condition that it has the linear momentum and angular 
momentum equal to those of the given body (Jardetzky [4], [7], 
[8], [9]). We restrict ourselves by mentioning here a theorem of 
Milankovitch [1] which seems to be thus far the best approach to 
solving the problem of migration of the Earth's poles. We can 
consider the deformation of the entire body of the Earth by assuming 
it to be elastic, plastic, or in any state differing from the rigid. 
On the other hand, we can assume that a liquid (or plastic or 
fluidal) core is covered by a solid crust which may change its shape. 

In both cases the asymmetric distribution of masses in the upper¬ 
most layer is the factor which plays the most important part in the 
displacement of poles. Let us consider the moments of inertia of 
this asymmetric shell. Every point at the free surface can be as¬ 
sociated with the value of the moment (/) taken with respect to 
the central axis passing through this point. Then the free surface 
becomes a field of /. Now, it has been proved by Milankovitch 
that the progressive displacement v of the Earth's pole occurs at 
each instant in the direction of the gradient of /, i.e., 

y=kVJ 

The progressive displacement of the pole is, of course, combined 
with the well-known periodic changes in the position of the pole 
(precession, nutation). 

The general conclusions of this theory concerning the rotation of 
a deformable Earth can obviously be applied to every celestial 
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body in the period of its solidification. The most important of 
them are: (a) the progressive displacement of the crust over the 
core during the period of solidification is due to the simultaneous 
rotation about an axis in the equatorial plane and adaptation of 
the mass to a continuously changing field of gravity forces (i.e., 
there exists a component of the angular velocity in the equatorial 
plane of the body); (b) the primary cause of such displacement is 
the action of the moment of centrigufal forces arising from the 
asymmetrical distribution of masses in the crust; (c) the adapta¬ 
tion is possible if the crust is flexible and the core displays the 
properties of a liquid; and (d) the currents in the interior can play 
in general an important part in the sliding of the crust over the 
fluid interior, but in such cases as that of the Earth, their influence 
on sliding seems to be of secondary significance. 

To get more insight into the character of these changes, small 
and slow as they are, many problems must yet be solved. 




CHAPTER IX 


Systems ComposeJ of Fluid and Rigid Parts 

9.1. Rigid Core and Fluid Envelope 

After probably more than four billion years of its existence, the 
Earth is still composed of solid and fluid parts. If we do not take 
into consideration the hydrosphere and atmosphere, the most ap¬ 
propriate model of it is a solid body with a cavity filled with a 
liquid. Of course, to simplify it, it is assumed that the mantle and 
crust are rigid. However, in some theories the existence of a rigid 
core in a liquid Earth is admitted for a certain period of its history 
The problem of a figure of such a system is reduced to that of strati¬ 
fication in the liquid part under given conditions. This problem 
has not yet been discussed in more detail, A short discussion of it 
is given in Poincar6's lectures [2] and in another investigation of 
Lichtenstein [5]. In the latter only, the general equations are set 
up to determine the exact form of the sea covering the Earth. No 
solution in terms of any definite functions is obtained. 

Even, if we restrict ourselves to the figures of equilibrium of a 
system composed of fluid and rigid parts, it is easy to see that there 
will be a great variety of problems due to the fact that different 
assumptions can be made concerning the composition of each part 
as well as the shape of the rigid core or crust. The simplest 
problem of this kind is obviously the absolute equilibrium of a 
fluid layer overlying a spherical rigid core. One solution is 
evident, namely, a stratification in concentric spheres, if the 
pressure at the free surface is zero. It seems that there are no 
essential changes required in the proof of the uniqueness and stabil¬ 
ity of the figure of equilibrium given for the case where there is no 

In a recent theory, the solid inner core is supposed to be separated from 
a solid mantle by the liquid part of the core. 
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rigid core. New complications arise, however, in the case of a 
relative equilibrium. For the fluid part of the system, the equa¬ 
tions (1.1.1)-(1.1.8) hold under the condition that U represents the 
sum of potentials Ui and of the fluid and rigid part respectively. 
The equation (1.1.5) must be written in the form 


(9.1.1) 


fU, + fU,-^ 




dp 

-h const. 

J 


and this condition holds for the volume Vi of the fluid. At the free 
boundary of it the condition p z= o or 

2 2 

(9.1.2) /C7^ +const. 

must be satisfied, but there are no such simple conditions for the 
inner boundary of F^. Whether we can assume that the gravity 
center of the core coincides with the center of gravity of the whole 
system depends on the shape of the volume of the rigid part F^ 
and on the density distribution in it. Only under the simplest as¬ 
sumptions about the density distribution and the shape of the 
core, can one expect to obtain a solution. Suppose that the core 
is a homogeneous sphere and that both the core and the fluid 
layer are rotating with a constant angular velocity. Suppose also 
that the whole system is isolated in space. The permanent exist¬ 
ence of a plane of symmetry, through the mass center of the whole 
system seems to be very plausible. Besides the symmetry with 
respect to this equatorial plane, the axial symmetry can also be 
expected. Maybe, it can be proved that figures of the type of 
Jacobi's ellipsoid are impossible in the case of a spherical core. 
However, a rigid ellipsoidal core with three unequal axes seems 
not to contradict the existence of such figures. For a uniform 
rotation of a mixed system (liquid and rigid or solid parts) a 
necessary condition is, as mentioned in section (1.2) that the axis 
of rotation coincides with one of the axes of inertia. This condi¬ 
tion is obviously satisfied in case of symmetrical figures defined 
above. Moreover, the spherical rigid core once set in a uniform 
rotation of this kind will continue to rotate about an axis fixed in 
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space, since the resultant moment of pressures at its surface as well 
as that of forces of attractions will be zero. At its surface r = 
and the potential 


(9.1.3) 


U. 




r 


0 


where is the mass of the core is constant. By (9.1.1) it is at 
r = r. 


(9.1.4) 


fU,+ 


co^ r dp 

2 j ^ 


+ const. 


but there is no reason to assert that the pressure (and density) have 
constant values along this surface. Thus, (9.1.4) becomes a new 
boundary condition which must be satisfied by a function which 
will determine the stratification. 

If we want to apply the method of Liapounov to a problem of 
this kind, the attempt should be made first to further simplify the 
conditions. For example, let us consider a very slow rotation of a 
slightly heterogeneous liquid in a layer around a homogeneous 
spherical core. Then, the angular velocity being very small, the 
level surfaces will probably differ but little from a set of con¬ 
centric spheres with the radius Let a = 1 be the value 

of parameter a corresponding to the free surface, i.e., to the sphere 
£. In the law of densities >c I + {a) (4.1.1) the parameter d 

has very small values again and the equations 
X ^ a{l sin '& cos ip 

(9.1,5) y = a{l + sin d' sin rp 

z ^ a{l 'd' 

represent the level surfaces They differ from concentric 
spheres and, therefore, one part of the set of will intersect 
the boundary of the core. Thus, at different points P{x, y, z) of 
the sphere f f ^ unequal values of a should be taken in express- 
sions (9.1.5). If the definition of the parameter a given in Section 
(4.1) is used, we have to assume now that the volumes of F ^ and E^ 
are equal for a ^ where must be determined by a special 
condition. The volume of F^ is composed of the volume of the 
core and that of the liquid layer. Changing the notations for 
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axes in (4.1.4) and other formulae in section (4.1), we obtain the 
condition 

3^0 + — «l) j (^ + == ^a^Q-\/Q 

Hence, if the condition (4.1.6) has to hold in this case as well, we 
must have a\Q^/Q = r\. The potential Ui in (9.1.1) could now be 
taken in the form (4.1.7) with the limits of the integral equal to 
and 1 . This would make it difficult to apply Liapounov’s 
method. We know, however, that the problem of figures of 
equilibrium in case of a liquid mass could be reduced to the equa¬ 
tion (4.3.16) when the potential is represented by the sum (4.3.4). 
The first term of this sum {W) corresponds to the constant density 
7 i = \. Since the potential due to attraction of a homogeneous 
body does not depend on its aggregate state, we could assume that 
W corresponds to the homogeneous figure composed of the rigid 
core (its density being 1) and the liquid of equal density 
filling in the rest of the volume. The remaining part of the sum 
Ui -f would be due to the excess of density dq>(a) in the volume 
Ftotai ~ ^oore • This would be equivalent to the modified law of 
densities in which we take 95 (a) = 0 for a < and x = 1 + dq3{a) 
for a ^ This condition will affect the term 0 in (4.3.4). 
Accordingly, we then have to assume that the function C = '0') 

which represents the deviation of level surfaces from the 
spheres E^, wiU be determined by the equation (4.3.16) 

in the interval (%, 1 ), i.e., in the liquid part of the system. Now, 
in (4.3.17) and (4.3,18) the 0(0) = 0 because of the vanishing of 
(p{a) at the mass center. Therefore, the difference between this 
problem and that concerning a liquid mass (Chap. IV) will be 
expressed by new values of the function 0 . We can assume again 
that the function C is expanded in the series 

oo 

i=l 


(9.1.7) 
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and apply transformations similar to those in Section (4.3). Lia- 
pounov's results show how difficult it is to calculate effectively the 
functions Ci which yield a set of approximations to an exact solu¬ 
tion. We will not make an attempt here to find even the first 
functions, since we do not expect that the solution of this problem 
will have any important applications. 

In the theory of tides, the free surface of the ocean in the state of 
equilibrium on the solid body of the Earth is taken as the surface of 
reference to determine the oscillation with respect to it. M. Bril- 
louin [1] pointed out, however, that according to his investiga¬ 
tions, this surface is unstable. His conclusion is that there is no 
stable permanent configuration of a liquid which has no relative 
motion with respect to a solid having a uniform angular velocity. 

We have mentioned certain very restrictive but necessary con¬ 
ditions of equilibrium of mixed systems, but the question of stabil¬ 
ity of such systems should be investigated in more detail. 

9.2. Liquid Mass and Floating Rigid Bodies 

The ratio of masses as well as that of densities are factors which 
determine whether a problem concerning the rotation of a system 
composed of liquid and rigid parts will be important in the theory 
of figures of celestial bodies or not. If, for example, the liquid 
mass is much larger than the mass of one or several rigid parts of a 
system and these parts are less dense, we can consider different 
configurations which will have this liquid mass with floating bod¬ 
ies. A period during which a planet or a former star can be treated 
as a large liquid mass differing but little from a figure of equilib¬ 
rium and partly covered with floating bodies seems to be one of 
the natural links in the evolution of celestial bodies. An objection 
is sometimes made that those substances that are present in the 
outer shell of the Earth became more dense through solidification 
and, therefore, had to sink in the remaining liquid mass when the 
Earth was cooled at its free surface. Since the increase in density 
with depth is a well-known factor, the question arises as to how far 
these solidified parts of the outer layer can sink. The stratification 
in a celestial body even if it is in the gaseous state seems to be an** 
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other well-established factor. Can we postulate that the law of 
density distribution is invariable during the whole history of a 
celestial body and that such solidified parts of the outer layer 
would always sink to a large depth? This assumption is made 
implicitly in some investigations dealing with the solidification of 
the Earth. Let us assume that after the time when small solids 
were moving toward the center in the vertical direction, melting 
and mixing with lower layers, the differentiation could reach such 
a stage that no more deep sinking would be possible. Then, there 
could be a larger part of an outer layer solidified such that its 
mean density was less than that of underlying layers. In case of 
the Earth such an assumption is supported by the isostatic mass 
distribution in the crust. Even if imperfect, the isostasy speaks 
in favor of the concept of floating continental blocks on a still hot 
and liquid planet. 

The solidification in the whole outer layer of a rotating liquid 
planet can occur simultaneously, and if no other factors will oppose 
this process, one can expect that the planet will be covered by a 
uniform crust. This is the viewpoint of many investigators who 
expect that the solidification to a great depth is a process progres¬ 
sing very rapidly. It could even be accomplished in a few- 
thousand years. On the other hand, the formation of a single 
continental block that floated on a liquid substratum and by 
bursting gave rise to present continents was the basic idea of 
Wegener's theory [1] of formation of continents. The existence of 
such a pangea or a protocontinent was not explained in this theory, 
but the author (Jardetzky [1], [4], [7], [8]) has shown that the 
hypothesis of a zonal rotation yields a simple reason for the forma¬ 
tion of a single continent on a rotating liquid planet. Since many 
features of the crust can be created by the displacement and de¬ 
formations of continents, as has been shown in Section (7.7), the 
motion of a body floating on a large liquid mass becomes an im¬ 
portant problem. The existence of continents and their deforma¬ 
tion affects, as it was pointed out earlier, the figure of a planet as a 
whole and is responsible for local departures from a more or less 
regular shape. 
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Unfortunately, because of the great complexity of the problem, 
there are but few attempts to treat it in mathematical form. Some 
remarks of general character can be found in investigations con¬ 
cerning this very remote part of the Earth’s history, (W. Thomson 
[1] Vol. IV, p. 189; Jeffreys [1]). They deal, for example, with 
the question of stability if there are many small floating bodies 
collected on one side of the planet. An attempt to discuss and solve 
the problem in an exact form was made by the author (Jardetzky 
[4]), but despite many simplifications, it remained intractable. 
It is obvious that in solving problems of this type, an explanation 
is sought of those characteristics of figures of celestial bodies which 
are due to the relative displacements of their parts. Therefore, 
the figures of equilibrium of a liquid mass with floating bodies are 
of secondary interest. It can be shown that very special conditions 
must be satisfied to secure the existence of a figure of equilibrium 
of this kind. One of these conditions concerning the position of 
the axis of rotation with respect to the axes of inertia was men¬ 
tioned in Section (1.2). 

It is easy to give some examples where this condition will be 
satisfied. Spherical polar caps homogeneous and regular in shape 
placed on a denser (nonrotating) liquid mass could remain in these 
positions but probably not indefinitely. The problem of a non¬ 
rotating liquid planet with a floating continent has been investi¬ 
gated by Jardetzky [3], [6], but as pointed out, the dynamical 
problem is actually of more importance. Less precise methods 
than those discussed in the preceding chapters were used to get 
more insight into the phenomena connected with the behavior of 
a material system composed of liquid and solid parts. By such 
methods, the existence of the so-called "Polflucht”-force in 
Wegener’s theory was analyzed. Several proofs were given in 
geophysics in order to show that, if there is a continent floating 
on a rotating liquid planet, a component of gravity force acting 
towards the equator must tend to displace this continent in this 
direction. Thus, if we assume the isostasy which is defined by the 
condition of floatation, such a component must appear because of 
the higher position of the floating body as compared with the 
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stratification in a totally liquid mass in equilibrium. That not 
all the conditions under which such a component will be directed 
toward the equator are always satisfied is apparent. The evalua¬ 
tion in case of the Earth yields a very small numerical value of this 
force. 

Future investigations will clear up the conditions for the exist¬ 
ence of horizontal as well as vertical forces acting on a floating 
continent. Perhaps even a rough solution will be sufficient in this 
case, but one has to keep in mind that this solution must satisfy 
the equations of motion of a liquid part and six equations for a 
rigid body. Of course, the chances of finding such a solution will 
be increased by the assumption that the motion of each part of the 
material systems differs but little from a pure rotation and that the 
configuration is almost a figure of equilibrium. 

It is evident that a better approximation to the actual figure of 
solid celestial bodies must be obtained by taking into account vis¬ 
cosity and elasticity of bodies, the factors like the drug produced 
by internal movements, etc. Therefore, we should be able to 
estimate the role played by each factor in a more precise way than 
it has been done before. This is, however, a very difficult task. 
As an example, we can mention an investigation of Prey [2] con¬ 
cerning the theory of land connections between the existing con¬ 
tinents. According to this theory, some parts of the crust have 
sunk in a viscous substratum and were covered by the ocean. 
Before that time they could represent connections between other 
parts of land. The reality of such process in some cases can hardly 
be doubted, but the mathematical analysis in the form given by 
Prey leads to some inconsistencies in the order of the magnitude of 
viscosity and in the time required for sinking. The problem dis¬ 
cussed concerns the radial displacement of a circular continent 
floating on a viscous liquid planet; the inertial terms in the equa¬ 
tions being neglected. The difficulties met in this problem may be 
easily understood, if one takes into account that to represent an 
island by means of a spherical function, terms up to the 16th order 
are required. 
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9,3. Rigid Bodies Having Cavities Pilled with Liquids 

A new type of problem concerning the motion of a system com¬ 
posed of liquid and rigid parts different from those discussed in 
preceding sections includes all cases where liquid fills the cavities 
in a rigid body. Mass or density ratios can be arbitrary in these 
problems. 

A problem of this kind approximates the present state of many 
planets probably quite fairly, but it must be assumed that the 
figure of the system has been determined by the conditions cor¬ 
responding to the time interval when deformations in the solid 
crust (or mantle) were still possible on a large scale. The adaption 
to a figure of equilibrium has been one of the most important 
among these deformations. The changes in the stratification in the 
core followed obviously the variations of the form of the crust, 
but all the deformations seem to have been small unless some 
catastrophe occurred. The fact that even the results of Clairaut's 
theory, which is nothing but a first approximation, agree fairly 
with the observations confirms this statement. On a large scale, 
the mass distribution in the Earth corresponds to that in a heter¬ 
ogeneous figure of equilibrium. In the crust, the isostatic layering 
is not quite perfect, but if, for example, the reduction for isostasy 
is taken into account when the ellipticity of the Earth is calculated 
from observation, the value agrees with Clairaut's value. The 
actual level surfaces of the Earth differ but little from the ellipsoids 
and even from a set of spheres. If certain corrections and a more 
precise theory are desired, the effect of different factors must be 
accounted for in further investigations. One of these factors is 
the probable variation of the angular velocity during the solidifica¬ 
tion of the mantle. 

The flexibility of the crust as a whole and its plasticity are well 
known to geophysicists. The tendency of some parts of the con¬ 
tinents to restore the isostatic layering disturbed by an excessive 
load during the last ice age is observed even at the present time and 
is a very important phenomenon. It proves the fact that the 
figures of planets must be better approximated under the assump¬ 
tion that their crusts are deformable. 
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Because of extreme complexity, one can only hope that in future 
investigations, solutions will be found in which use will be made of 
methods similar to those of Liapounov (see Chap. IV). The series 
expansions given for the potential of a body are not restricted to a 
heterogeneous liquid. It is evident that they depend on a given 
law of densities in a body and the fact that the stratification differs 
but little from an ellipsoidal one. The representation of the poten¬ 
tial can be used in problems concerning rigid, solid, and plastic 
masses. 

As to the motion which differs from a rotation about an axis 
fixed in space and in the body, the investigations of Joukowski 
[1], Volterra [1], and Stekloff [2] yield basic results. In these 
investigations, the shape of a rigid body and the form of a cavity 
filled with a homogeneous liquid are given, and the ''deformations'' 
are admitted only in the sense of internal movements. 



CHAPTER X 


Fluid Mass and Centers of Attraction 

10 J General Problem 

The small changes in the figure of a rotating fluid mass due to 
the existence of some other bodies can also be treated by the method 
of Liapounov if the basic conditions required by this method are 
satisfied. In the problem of the figure of a liquid planet which is 
affected by the attraction of other members of a given ''solar'' 
system, the most precise method has been applied by Lichtenstein 
[2]. In earlier investigations, certain particular cases of the prob¬ 
lem have been treated by other methods. The double-star problem 
in its simplest form when the dimensions of the second star are 
neglected, that is, it is assumed that this star acts as a simple mass 
center, or the tidal problem are the well-known examples. 

The equations of motion of a rotating fluid mass attracted by 
several centers have the form (1.1.3) again, but the potential U is 
now the sum of the potential of the fluid mass itself and C/" due 
to the masses located at z^). The integral of these 

equations which will yield the figures of the fluid mass is, therefore, 
(1.1.5) or 

, co^s^ f dj> 

(10.1.1) /([/' + C/") H-= Ml + const. 

2 J 0^ 

We assume that the motion of the fluid mass is a rigid-body type 
of rotation, but on changing the second term of the left-hand 
member, as has been done in Chap. VII, we can also account for a 
zonal rotation. Only the simpler case is to be considered here, but 
it must be pointed out in advance that even a single, far-located 
mass center will be able to change the state of a uniform rotation 
of either a perfect or viscous fluid. In the first approximation, 
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equation (10.1.1) holds for each instant. However, the term 
can depend on time, and, in this case, the function f which will 
determine the departure of a level surface from an ellipsoid 
becomes a function of time. Moreover, in case of an isolated fluid, 
we could assume that the axis of rotation passes through the mass 
center of the whole mass. In the new problem, the mass center 
of the fluid is moving in general with respect to the mass center 
of the whole system. It is obvious that different assumptions can 
be made about the order of magnitude of the masses involved as 
well as about the ratio of the largest dimensions in a fluid mass and 
the distances between the mass centers. Accordingly, different ap¬ 
proximations must be used, and we shall mention few examples. 
In the double-star problem, the case of a uniform rotation of both 
components with respect to the axis passing through the mass 
center of the system had been considered. If the masses are m and 
and the x axis is directed toward the mass center the dis¬ 
tance between and C (the mass center of the first component) 
will be where I = CC^. Since in (10.1.1) s^is the 

distance from the axis of rotation, the second term in this equation 
takes the form 


2 


lx -— i\ 4-j 

A m + m I 


However, for celestial bodies, the periods of revolution and rota¬ 
tion are seldom equal, the angular velocity due to revolution being 
of a high order than that of rotation about the axis of the figure 
(exceptions: Mercury (?), Venus (?), Moon in the system Earth- 
Moon). 

Let us assume now that the mass center of the whole system is 
in the close neighborhood of the mass center of the liquid mass 
(the figure of the Sun affected by planets or the figure of a liquid 
planet affected by the satellites). Then in (10.1.1) the distance s, 
can be neglected. Another simplification concerns the expansion 
of the potential due to the mass center in a series. If are 
distances of from C (which coincides with C,) and l[ = 
where P is a point of the fluid mass, we have 
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(UM 

( 10 . 1 . 2 ) 

Since 

l? = — 2lir cos 

for the mass we obtain 

(10.1.3) u, = -^\l~2jCOsy,^--^ 

In spherical polar coordinates 

(10.1.4) cos = cos '& cos 'd'^ + sin d' sin cos [xp — xp^) 

and the last factor in (10,1.3) can be expanded in the series (see 
equation (1.4.13)) 

oo 

(10.1.5) 2 - P„(cosy,) 

n-0 

where are the Legendre polynoms and 



By (10.1.2)~(10.1.5) we write, if there are q mass centers 

<7 ^ yyi, 

(10.1.6) U"=% '^-~P^{cos^•cos'&^+sin'&sm'&iCOs{^t>—^|>i)) 

i=l n'=»0 f'i 

where (4.1.2) 

(10.1.7) 

{(^+1) sin®?? cos2v'+(^+S') sin®?? sin®v;+e cos®??} 

The form of the series expansion of the potential U' of the liquid 
mass is not affected by adding the term U" in (10.1.1) and all 
transformations of Liapounov given in Chapter TV hold again. 
Therefore, the functional equation (10.1.1) can be reduced again to 
an integro-differential equation and subsequently to the set of 
integral equations of the type investigated in Chapters IV and VII. 
This integro-differential equation is by (4.3.16) 
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( 10 . 1 . 8 ) 


1 r ^'da' 

4jra® J Z)(a, 1) 


= W + P{a) + 


1 

4^a^A 


U" 


where W and P(a) are expressions (4.3.17) and (4.3.18) or the 
more general expressions given in Liapounov’s theory for ellipsoids 
with three unequal axes. In order to obtain from (10.1.8) the 
integral equations yielding functions which can be used for calcu¬ 
lation of different approximations, we have to determine some 
new parameters. In the problem of equilibrium of a single fluid 
mass, the parameter S in the density law or k in case of a zonal 
rotation have been considered as small quantities of the first 
order. In (10.1.6) however, the order of magnitude of three new 
factors: m, l/^,-, and rjl^ can be taken in an arbitrary manner. 
The first term of this series expansion is injl^ and is, of course, in 
general variable. If we assume that the paths of masses around 
the point C, are circular or that the centers are fixed in space, 
the ratio mjlf is constant. It is small if either the masses nif are 
small or the distances large or both. This factor does not affect, 
however, the relative magnitude of terms of the sum (10.1.6) which 
depends on the ratio r/^. If we want to apply the method of 
Liapounov to the equation (10.1.8) we have to assume that the 
function f can be expanded in a multiple series of the form 


(10.1.9) 




where the small parameters are chosen in such a way that they 
can represent the effect of attraction of the mass m^. To avoid 
unnecessary complications, and d must be of the same order of 
magnitude and ?oo...o = ^' Because of (10.1.7) from which it 
follows that rjl^ is of the order Vg + Ijl^, Vg + 1 being the largest 
half axis of the free surface, it seems to be natural to assume that 
for the first approximation, the factor 


( 10 . 1 . 10 ) 




Wi Vg -f 1 

k li 


is of the order of d. The assumptions about the order of magnitude 
of three independent factors m^, 1/4-, and r/4 are combined in one, 
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but it is evident that this is not the only possible condition. Now, 
by substituting (10.1.9) in (10.1.8) the integral equations of the 
type (7.3.14) or (7.3.7) for the coefficients can be derived by 
transformations similar to those in case of a single parameter d. 
To find effectively these coefficients very extensive calculations 
are required even for the first approximation 

(10.1.11) C = + Coi...o<52 H- . . . + 

where the number of terms is g + 1, if there are q mass centers. 
The coefficients foo...i = is given by Liapounov (see Chap. VII) 
and the others can be calculated in the way used by the author in 
the problem of the zonal rotation of a homogeneous mass (Jardetz- 
ky [3]). 

10.2. Recent Investigations 

The classical problem of Roche concerning the figure of a homo¬ 
geneous liquid mass subjected to the attraction of a far-removed 
mass center was recently generalized by Agostinelli [1]. He as¬ 
sumed that there is a set of such mass centers and that the changes 
in their positions can be neglected when a time interval is considered 
corresponding to several rotations of the liquid mass. All mass 
centers are supposed to be located in the equatorial plane of the 
liquid mass. It could be proved that under these conditions, a rigid 
body type of rotation is impossible. However, the free surface of the 
hquid can have the shape of an ellipsoid with three unequal axes, 
which do not change the magnitude or direction with respect to the 
centers of attraction. Then, the particles of the liquid move along 
elhptical paths. 

The fundamental equation of the problem has the form 
Uf y^) -f- Cl (a:® — -f- Cgjy® — z^) = const, 

where the terms in C^ are due to the centers of attraction. Agosti¬ 
nelli [3] had also considered the case of a rotating heterogeneous 
liquid taking Roche's law for density distribution. To treat this 
case, he assumed that there is a complementary motion due to a 
dilation of the liquid. If the stratification differs but little from an 
ellipsoidal one, the potential can be expressed in terms of the polar 
and equatorial superficial elhpticities. This theory applied to the 
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Earth which is subjected to the lunisolar attraction yields for the 
equatorial elhpticity the limit 1.67x10“®. 

On making use of Poincare's analysis, Agostinelli [3] found also 
that there are some critical ellipsoids which can give rise to new 
figures of equilibrium if slightly deformed. These configurations 
are obviously not figures of equilibrium in the strict sense, since 
the results hold under the assumption that the centers of attraction 
do not change their positions. 

A generalization concerning the number of centers of attraction 
leads to another problem which presents a great interest. If we 
assume, namely, a very large number of such centers distributed in 
a ring concentric to a liquid mass, we obtain the problem of the ring 
of Saturn or of a ring-shaped nebula with a central star. The 
classical investigations concerning the first problem were discussed 
in the work of Appell ([1] Vol. IV) in more detail. We shall only 
mention a recent report by Nadile [1]. He could show, namely, 
that the free surface of the central liquid mass can take the form of 
an ellipsoid if the radius of the ring is very large and its cross sec¬ 
tion very small. This conclusion holds when the terms of an order 
higher than are neglected, K being the ratio of the maximum 
dimension of the central body to the radius of the ring. 




CHAPTER XI 


Figures of Compressible Masses 

11.1. Stratification in Stars 

The theories of figures of celestial bodies discussed in preceding 
chapters reached a very high level of mathematical perfection. 
There is no doubt they became a solid base for explanation of 
many features or properties of any kind of celestial bodies (sohd, 
liquid, and gaseous). We have seen, however, in case of the zonal 
rotation, that the approximation based on equations of Mechanics 
only requires an improvement. We know, for example, that the 
equilibrium in a star is maintained by the gravity, the pressure 
gradient, and by the radiation and propagation of energy in it. 
From the mechanical viewpoint, the theories that should yield a 
better explanation of stratification in stars, must take into account 
the much higher compressibility of gases. On considering such a 
new factor, the investigators of figures of compressible masses 
made use of different methods. 

Usually, it is assumed that stars and some planets or their at¬ 
mospheres are composed of an ideal gas. The equation of state has 
the form = RT. There are two cases in which some problems 
can be also treated by using the theories mentioned in preceding 
chapters. For adiabatic processes, we have ^iy'>'=const. and for poly¬ 
tropic variations it is ^ [n — index of polytropy). Even 

if we have a so-called generalized polytropic gas ^ = **('')) 

the relationship between the pressure and density is of the form 
■p = (p{ii), and the motion of the mass is barotropic. The funda¬ 
mental equation of the theory of equilibrium and of a zonal 
rotation (1.1.5) holds for adiabatic as well as for polytropic states, 
and the stratification is determined by it. However, there is a 
difference in one point. We no longer can make use of the condi- 
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tion of equal volumes when deforming a figure of equilibrium, 
since in a compressible gas a change in the pressure field can pro¬ 
duce an essential variation of density and volume. Thus, as a 
supplementary condition, equahty of masses must replace the 
equality of volumes. 

Instead of these equations in the problem of stratification in a 
gaseous celestial body, other conditions can also be used, e.g., for 
the absolute equihbrium, there are three equations for three 
variables: f, U (assuming a spherical layering). 

( 11 . 1 . 1 ) 1 ^ = /.^, V^U=-4:nK, f = 

' % dr dr 

The first equation is a special case of (1.1.2), the second is the 
Poisson equation. The substitution 

^ —I. ,. J y 

( 11 . 1 . 2 ) « = = — 

‘ktf a 

transforms then Poisson’s equation into the equation of Emden 



which yields a function 6 = d(i), i.e., the law of densities x ~ ;^(r) 
for a given n(p{) and boundary conditions. 

In the next approximation, as is well known, the pressure of 
radiation p* is added in the first condition of equilibrium (11.1.1). 
Then, if is the mass bounded by a sphere of radius r, 


(11.1.4) 


1 dP 
K dr 


1 d 
H dr 


[P + P^) - 






This is equivalent to the assumption that an additional force, 
grad^*, is introduced in the equation (1.1.2). It is not our pur¬ 
pose here to discuss all the hypotheses about the generation of 
energy and its propagation in the gas and absorption, which would 
make it possible to express in terms of r. Different models were 
considered in Astrophysics. Each definite form of p^{r) is con¬ 
nected with certain assumptions about physical processes involved 
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and approximates more or less the actual state in stars. The 
nuclear processes are at the present time the most important point 
in this theory. The problem of stratification in stars in this form 
is not a question of how to determine the shape of level surfaces or 
surfaces of equal density or pressure. This shape is given and the 
density (or pressure) is to be found. 

Further investigation of the stratification meets a new difficulty 
in the following fact. If in the equation of state a new variable, 
as, for example, the temperature T is introduced, the equations of 
Hydromechanics do not represent a more complete system of com 
ditions that will determine all the variables {T included). It has 
been pointed out by Duhem [1] that the required sixth equation 
will be given by Thermodynamics. Nevertheless, as far as is known 
to the author, no general use has been made of this suggestion. 
Instead of this, several general conclusions were drawn directly 
from the equations of Hydrodynamics. If x is not a function of a 
single variable 'p, the right-hand member in (11.1.4) is not an 
exact differential and the equation (11.1.5) cannot be obtained 
even if o) is constant. The three sets of surfaces h = const., 
p = const., and U -\-Q ^ const, now no longer coincide, and the 
term stratification loses its exact meaning. The motion of such 
a mass is called barocline, but it must always be specified what is 
meant by stratification. Some useful relationships can readily be 
obtained. Writing the equations (1.1.3), for example, for cylin¬ 
drical coordinates and assuming the rotational symmetry (which 
is proved to exist), we have the first equation in the form 


= / 


dU 

ds 


1 dp 
K ds 


Eliminating p from this equation and the third (1.1.3), the con¬ 
dition for CO is 


(11.1.6) ?tJco2 


/ r idK dU 

5 J \ 0s dz 


dx dU 
dz ds 


j + funct. (s) 


(Pizzetti [1]). If we eliminate U from (1.1.3) it is, for example, 
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dco^ _ ^ (-] 

^ dz dx dz \k/ dz dxxH/ 

This and two other equations show that 

(11,1.6) grad^co^ x r = xVp 

where grad, is the partial gradient in the ^ direction. Thus, the 
isobaric and isosteric surfaces coincide only if dco/dz ~ 0, The 
condition (11.1.5) shows that in general co == (jo{s^, if the sym¬ 
metry with respect to the equatorial plane is taken into account. 
The conclusion drawn from the condition (11.1.6) seems to be 
rather surprising. The fact that isobaric and isosteric surfaces do 
not coincide in a fluid mass is due to a special form of a physical 
law concerning the property of this mass, namely, to the equation 
of state, for example, p — p{x, T), As interpreted in many in¬ 
vestigations, a rigid body rotation of such a mass is a condition 
sufficient to make these two sets of surfaces coincide again, without 
any assumption about the temperature. Of course, there will be 
no contradiction, if isothermic surfaces correspond to the tempera¬ 
ture distribution, say, T = const, on the surface k = const. The 
surfaces p = const, would belong to the same family, and each of 
the variables could be expressed in terms of a parameter like the 
parameter a of Liapounov: % = x{a), p — p{ci), T — T{a). Then, 
the rotation would be barotropic. It seems to be necessary, how¬ 
ever, in order to avoid a contradiction in case of barocline rotation, 
to take care of the fact that such zonal movements represent only 
one part of a general circulation, the second component of which 
consists in meridianal currents. In theories concerning the strati¬ 
fication in gaseous stars, it is usually assumed that the actual state 
is a convective equihbrium. The rising warmer masses expand 
adiabatically, the sinking are adiabatically compressed, but ac¬ 
cording to the assumption just mentioned, the shape of layers 
K = const., p = const, is not violated by these currents. If this 
is a precise description, the components of acceleration due to 
these radial movements should be included in the fundamental 
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equations, even if we neglect certain differences in motion along a 
parallel. Thus, (1.1.3) will take the form 




(11.1.7) 


dx 
a, ~ 


1 dp 
•K dx ’ 


-co^y+a^ = 


dU I dp 
dz X, dz 


or in case of rotational symmetry 

dU 


— co^s + «, = /■ 


(11.1.7') 


ds 




I df ^ du 

H OS OA 

dU 1 dp 


% dz 


dy K dy 


1 dp 


A being the longitude, if we eliminate now U from (11.1.7) we ob¬ 
tain 


-1 

f-1 

1 

-f-1 


9a)2 

day 

da. 

ay' 

\k] 

dz 

dz \h) 

dy 

^ dz^ 

dz 

dy 


f-1 

1 ^P 

-f-1 

dp 

d(o^ 

"17 + 

da. 

da,, 

\k] 

dx 

dx\x) 


dx 

dz 

T| 

f-1 

.dp 




da^ 

day 

dx 


’ dy 

dy\>cj 

' dx 


dy 

dx 


These conditions will coincide with (11.1.6) only if 

(11.1.9) curia = V X a = 0 

In general, curia ^ 0 and instead of (11.1.6) we have 


(11.1.10) V — X Vp ^ grad^o)^ X r + curia 

K 

The vanishing of grad^ co^ now would not lead to the coincidence 
of isobaric and isosteric surfaces. This conclusion also holds for 
viscous liquids. 

We mentioned the method of Bondi and Lyttleton for the in¬ 
vestigation of internal movements in a viscous liquid enclosed in 
a rigid core (Chapter IX). To discuss the problem of such internal 
movements in a gaseous star, one must also take into account the 
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fact that its stratification and the boundary surface are unknown. 
For compressible masses certain approximations have been made 
only for the case of small deformations of polytropic spheres due 
to a very slow rigid-body rotation. The assumptions are made, 
namely, that 

(11.1.11) j> = = H{r, n = const. 

If we put cos ?? = and make use of the substitution represented 
by second and third equation (11.1.2), the density h == 

Since the stratification, i.e., the surfaces k = const., is changing 
when the angular velocity co is varying, it is assumed that for small 
values of o) the density can be expanded in a power series 

(11.1.12) + 

Suppose that in the first approximation on neglecting and 
higher powers of co we can take 


(11.1.13) 




0 ^^ 4 - w - 

^ 27tf 


where 6 is the solution of (11.1.3), = nd^~-hp and y == 'ip{$y p). 

If U and j) are eliminated from (11.1.1) an equation for the 
function ^ will result. This equation can be written in the form 
(Milne [1], Chandrasekhar [1] and other papers) 


(11.1.14) 



+ 


1 A 

dp L 


(1-/^2) 


dip" 

dp_ 


— I 


if we make use of the expression for in spherical coordinates. 
A solution ip = p) of this equation will determine the strati¬ 
fication (11,1.12). The free surface that corresponds to the bound¬ 
ary condition = 0^ = 0 is a spheroid. To solve the equation 

(11.1.14) Chandrasekhar assumed that ip can be represented by the 
series 

oo 

i^l 

where are Legendre polynomials. On inserting (11.1.16) in 
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(11.1.14), a set of equations is obtained to calculate the functions 
(^). In order to find the coefficients A ^ Chandrasekhar introduc¬ 
ed two new conditions. He puts equal expressions for the interior 
and exterior potential. The first derivatives of potential are also 
assumed to be equal at the boundary of the first Emden sphere 
I == which should be the figure of a nonrotating mass. Under 
these conditions, Chandrasekhar’s solution is 

0)2 r ... 6 


(11,1.16) 6=6{^)+- 


2>7zf 




6 




It has been pointed out by Krat [1] that the conditions introduced 
by Chandrasekhar can cause serious errors, since the action of 
masses exterior to the sphere ^ ^ seems to be neglected. More¬ 
over, the question arises, whether both series expansions are con¬ 
vergent on this sphere or not. The remark of Poisson about a 
spherical layer (1.4.15) in which the expressions (1.4.13) given by 
Laplace are not valid could hold in the case under consideration 
too. Attempts were made to close the gap between those two ex¬ 
pressions, but the question does not seem to be clear enough. The 
convergence of the basic series (11.1.12) in this theory of figures of 
compressible masses is just assumed and not proved subsequently 
to be correct. Krat [1] suggested a generalization of Chandrase¬ 
khar's results for the case of a zonal rotation 


(0^ = CqCoI + COo = COod), COi = COi{i) 

In a previous investigation (Krat [1]) he had pointed out that if a 
given law of rotation of a polytropic gaseous mass is incompatible 
with conditions of its equilibrium, this will hold for any gaseous 
mass (generalized polytrop). As an example, the law co^ = 
dz 0 ) 1 $^ is given which does not satisfy the conditions of equilib¬ 
rium of a polytropic mass (despite the fact that the integral of 
(11.1.7') for == 0 does exist). 

For compressible gases in stars, a barocline motion seems to be 
the only possibility. Therefore, it will be very important to find 
precise methods which will enable us to determine all three sets of 
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surfaces, namely, k = const., p = const., and + /?7 = const, 
with any desired approximation. 

11.2. Remarks on The Double-Star Problem 

In its general form, this problem presents such difficulties that 
no method could yet be given to obtain its solution. We have seen 
in the preceding section some simplifications used in the case of a 
slow rotating isolated mass of a compressible gas. The difficulties 
arise first because of the fact that the components of a double star 
are compressible and heterogeneous. Besides the revolution about 
a common mass center, the rotation of both stars with different 
angular velocities is an important factor and the tides may domi¬ 
nate all other phenomena in a close system, Darwin considered 
the case of a rigid body rotation of the whole system. As a first 
approximation for homogeneous incompressible components, the 
ellipsoidal configurations were determined. However, one can 
expect that in a gaseous star the density is essentially increasing 
from the free surface towards the center. Therefore, the larger part 
of the mass is distributed in-the ‘'core,’’ and the actual confiqura- 
tion will differ from an ellipsoidal mass. Bearing this in mind, 
Roche’s model gains in its importance. For a double star, the 
potential due to centrifugal forces has to be combined with gravity 
potentials of two mass points. This model corresponds then 
actually to the restricted three-bodies problem. A set of equi- 
potential surfaces being determined from it, certain indications 
about the paths of particles in stellar atmospheres have been 
obtained, and the discharge of the stream through the “criticar’ 
surfaces have been made probable. But the stars differ very 
much from Roche’s model, and better approximations are required. 

For polytropic stars, two particular problems became important. 
The tidal problem which is reduced to a rectilinear motion of the 
first component in the direction towards the second has been in¬ 
vestigated by Chandrasekhar [1]. For the solution, series expan¬ 
sions similar to those mentioned in the preceding section were used 
and the coefficients determined in terms of a parameter v equal 
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to the ratio of the radius of an undisturbed configuration to the 
distance between the mass centers of star components. It could 
be shown that up to r®, the deformations due to such tides can be 
superposed on those arising because of the adjustment of a double 
star to a shape corresponding to a rigid-body rotation. The figures 
of stars are approximated by ellipsoids if the terms of the order 
V* and are neglected. Further theoretical progress was made by 
Chandrasekhar [1], Krat [1], Russel [1], Cowling [1], and Walter 
[ 1 ]. 

The double-star problem has been treated by a different method 
by Lichtenstein [1]. Upon giving a far more precise formulation of 
the problem, only the general conditions were given in a mathema¬ 
tical form. 

Liapounov’s method can also be applied to the problem of a 
double star, and we shall write the general equations. If we assume 
that the densities of both components are functions of pressure 
>^2 — f 2 {f’)> 'th-6 equations of motion of fluids in a double 

star are 

_ d^r, rd^\ . _ 


because of (1.1.2). Hence on denoting by 
accelerations 

+ = (^+001 


Qi the potential of 


const. 


on a level surface in each component. Since there are only two 
masses involved, the potential has the form 


(11.2.3) U^ = U+U 

i.e., it is the sum of potentials due to two components of a double 
star. In order to make use of the series expansions of Liapounov, 
certain conditions must be satisfied. The basic series expansion 
(4.1.15)-(4.1.16) holds for small (not infinitely small) values of the 
function f which determines the distortion of an ellipsoidal strati¬ 
fication. This can occur either for mean variations of density 
from the center to the free surface of each component, if there are 
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no restrictions about the magnitudes of axes of ellipsoids of 

£i(V^i+ 1, + Qi) £'2 respectively (see Section 

(4.1).) or, if the rotation of each component is very slow. In the 
second case, the ellipsoids and £2 approach spherical shape and 
even for a very large increase in density from the free surface to the 
center, the level surfaces can be spheroids, i.e., C is very small. 
In both cases the distance between the mass centers of the com¬ 
ponents must be large enough in order not to produce an essential 
distortion in the stratification of each. Under these conditions, 
we can write (4.1.15)~(4,1,16) for each component, i.e., 

(11.2.4) 

where the terms are of the order Ci and of the order tg* Cx 
and ta being functions which in (4.1.2) will determine the stratifi¬ 
cations of the first and second star, respectively. The functions U 
will be replaced by V and 0, Section (4.3), if the laws of densities are 
given in the form 

(11.2.5) Hi = 1 + d^(pi (a^) 

To simplify the problem one could assume that = (Jg = d but 
the functions (pi and (p 2 are different and the parameters and aQ 
are to be defined in general by two different conditions. If we 
take the definition of as given in Section (4.1), these conditions 
are the equations (4.1.6) of Liapounov written for each com¬ 
ponent, i.e., 

( 11 . 2 . 6 ) fCidcr^ 

If we consider the stars, equality of masses must replace the 
equality of volumes enclosed by and as mentioned above 
(see Krat [1], p. 148). 

On writing the volume element (4.1.4) in the fonu 

(11.2.7) dV = ^ + mdada 

3 da 

and putting in this expression f = 0 for the ellipsoid E^ we obtain, 
because of the new definition of the parameter a, 



FIGURES OF COMPRESSIBLE MASSES 


169 


{11,2.8) dV = Zl >c da da = 4jrZl M^da 

If we substitute the expression (11.2.6) the term equal to unity 
gives rise obviously to all terms in (4.1.6). The correction due to 
the second term in (11.2.5) can be easily computed, and we have 
two conditions for functions Ci follows: 


(11.2.9) '^ 

^ Jo 

This equation must be used instead of (11.2.6) when the func¬ 
tional equations (11.2.2) will be transformed in integro-differential 
ones and subsequently into integral equations. 

The acceleration potential may present an extreme complex¬ 
ity in more general cases if it does exist. Some approximate ex¬ 
pressions for it can be written only by neglecting certain parts of 
acceleration. By assuming that invariable configurations can 
exist, a rigid body rotation of the whole system yields such a simple 
expression 

( 11 . 2 . 10 ) = ^ 


where denotes the distance of a particle from the axis of rotation 
passing through the mass center of the double star. However, the 
motion of mass centers will not strictly be a Kepler motion. If 
this approximation is made, it may be put that 


( 11 . 2 . 11 ) 




where and are masses of the double-star components and Iq 
the distance between their centers. This distance is the factor 
which affects the approximation in equations (11.2.2) and is used 
for the definition of the parameter v mentioned above. It is 
evident that the ratio of the largest linear dimension in a com¬ 
ponent to the distance 1^ will determine the limit term to be used 
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in an approximation. In a model with homogeneous components, v 
is the only parameter to be used even for close binaries. The rela¬ 
tive order of parameters and v has to be determined for hetero¬ 
geneous components, since to transform the equations ( 11 . 2 . 2 ) we 
must use the series expansions of the form 

and combine the terms of the same order. 

Two equations for functions Ci and Ca 


(11.2.13) 


fU fU QI f + const. 
fU fU Q 2 = f —^ + const. 


for Si 
for Sg 


are of course very simplified in Roche's model in which the po¬ 
tential of another component is reduced to that of a mass point, 
and the fluid is considered as a compressible homogeneous liquid. 
Then, we have 

— Mo 

fXJ f-jj- + =-h const. for Si 

(11.2.14 ^ 

= i) 

f~ + fU + — + const. for 

I ^2 

Usually the figures are approximated by ellipsoids, but a progress 
can be made on substituting in (11.2.14) Liapounov's series ex¬ 
pansion for U in terms of ^i or for U in terms of C 2 . In case of 
(11.2.14) the functions Ci and fg satisfy two separate equations, 
but even in this simplest case the solution of an equation of the 
form (11.2.140) will require very long calculations. Therefore, we 
restrict ourselves now to the last remark which concerns the char¬ 
acter of condensation in a star as given by the law (11.2.5). It 
may be written in a form representing any degree of condensation 
if the parameter 6 is not necessarily a small number. For the law 
X = 1 + and any n > 0 limits for the density are 

determined at the center and at the free surface. They are 
= 1 4 . = 1 + but since 0 < a < 1 , the 
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larger is the value of n; the slower variation of density occurs in 
the ' ‘core” of the star and the more rapid decrease in the “mantle.” 
Maybe this law of densities would better represent the actual 
density distribution in a star and would not increase too much the 
mathematical difficulties of the solution, 
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